OSA's Digital Library

Optics Express

Optics Express

  • Editor: Michael Duncan
  • Vol. 12, Iss. 1 — Jan. 12, 2004
  • pp: 90–103
« Show journal navigation

Matrix analysis of microring coupled-resonator optical waveguides

Joyce K. S. Poon, Jacob Scheuer, Shayan Mookherjea, George T. Paloczi, Yanyi Huang, and Amnon Yariv  »View Author Affiliations


Optics Express, Vol. 12, Issue 1, pp. 90-103 (2004)
http://dx.doi.org/10.1364/OPEX.12.000090


View Full Text Article

Acrobat PDF (271 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

We use the coupling matrix formalism to investigate continuous-wave and pulse propagation through microring coupled-resonator optical waveguides (CROWs). The dispersion relation agrees with that derived using the tight-binding model in the limit of weak inter-resonator coupling. We obtain an analytical expression for pulse propagation through a semi-infinite CROW in the case of weak coupling which fully accounts for the nonlinear dispersive characteristics. We also show that intensity of a pulse in a CROW is enhanced by a factor inversely proportional to the inter-resonator coupling. In finite CROWs, anomalous dispersions allows for a pulse to propagate with a negative group velocity such that the output pulse appears to emerge before the input as in “superluminal” propagation. The matrix formalism is a powerful approach for microring CROWs since it can be applied to structures and geometries for which analyses with the commonly used tight-binding approach are not applicable.

© 2004 Optical Society of America

1. Introduction

Coupled optical resonators are becoming important in nonlinear optics research as well as in telecommunication applications [1

1. A. Yariv, Y. Xu, R. K. Lee, and A. Scherer, “Coupled-resonator optical waveguide: a proposal and analysis,” Opt. Lett. 24, 711–713 (1999). [CrossRef]

, 2

2. Y. Xu, R. K. Lee, and A. Yariv, “Propagation and second-harmonic generation of electromagnetic waves in a coupled-resonator optical waveguide,” J. Opt. Soc. Am. B 77, 387–400 (2000). [CrossRef]

, 3

3. N. Stefanou and A. Modinos, “Impurity bands in photonic insulators,” Phys. Rev. B 57, 12 127–12 133 (1998). [CrossRef]

, 4

4. D. N. Christodoulides and N. K. Efremidis, “Discrete temporal solitons along a chain of nonlinear coupled microcavities embedded in photonic crystals,” Opt. Lett. 27, 568–570 (2002). [CrossRef]

, 5

5. S. Mookherjea and A. Yariv, “Kerr-stabilized super-resonant modes in coupled-resonator optical waveguides,” Phys. Rev. E 66, 046 610 (2002). [CrossRef]

, 6

6. J. E. Heebner and R. W. Boyd, “‘Slow’ and ‘fast’ light in resonator-coupled waveguides,” J. Mod. Opt. 49, 2629–2636 (2002). [CrossRef]

]. Systems consisting of a few coupled resonators, say 1<N<5, have been proposed for optical filtering and modulation [7

7. C. K. Madsen, “General IIR optical filter design for WDM applications using all-pass filters,” IEEE J. Lightwave Technol. 18, 860–868 (2000). [CrossRef]

, 8

8. G. Lenz, B. J. Eggleton, C. K. Madsen, and R. E. Slusher, “Optical delay lines based on optical filters,” IEEE J. Quantum Electron. 37, 525–532 (2001). [CrossRef]

, 9

9. B. E. Little, S. T. Chu, W. Pan, D. Ripin, T. Kaneko, Y. Kokubun, and E. Ippen, “Vertically coupled glass microring resonator channel dropping filters,” IEEE Photon. Technol. Lett. 11, 215–217 (1999). [CrossRef]

]. On the other extreme, “large” systems, say N>10, can be regarded as a new type of waveguide termed Coupled-Resonator Optical Waveguide (CROW) with unique and controllable dispersion properties [3

3. N. Stefanou and A. Modinos, “Impurity bands in photonic insulators,” Phys. Rev. B 57, 12 127–12 133 (1998). [CrossRef]

, 1

1. A. Yariv, Y. Xu, R. K. Lee, and A. Scherer, “Coupled-resonator optical waveguide: a proposal and analysis,” Opt. Lett. 24, 711–713 (1999). [CrossRef]

, 2

2. Y. Xu, R. K. Lee, and A. Yariv, “Propagation and second-harmonic generation of electromagnetic waves in a coupled-resonator optical waveguide,” J. Opt. Soc. Am. B 77, 387–400 (2000). [CrossRef]

, 10

10. M. Bayindir, B. Temelkuran, and E. Ozbay, “Tight-binding description of the coupled defect modes in three-dimensional photonic crystals,” Phys. Rev. Lett. 84, 2140–2143 (2000). [CrossRef] [PubMed]

].

The “large” chains (CROWs) have been previously analyzed using a tight-binding formalism [1

1. A. Yariv, Y. Xu, R. K. Lee, and A. Scherer, “Coupled-resonator optical waveguide: a proposal and analysis,” Opt. Lett. 24, 711–713 (1999). [CrossRef]

]. In the tight-binding method, we approximate the electric field of an eigenmode E K of the CROW as a Bloch wave superposition of the individual resonator modes E Ω [1

1. A. Yariv, Y. Xu, R. K. Lee, and A. Scherer, “Coupled-resonator optical waveguide: a proposal and analysis,” Opt. Lett. 24, 711–713 (1999). [CrossRef]

],

EK(r,t)=E0exp(iωKt)nexp(inKΛ)EΩ(rnΛẑ),
(1)

where the nth resonator in the chain is centered at z=nΛ. Under the assumption of symmetric nearest neighbor coupling, the dispersion relation of the CROW is [1

1. A. Yariv, Y. Xu, R. K. Lee, and A. Scherer, “Coupled-resonator optical waveguide: a proposal and analysis,” Opt. Lett. 24, 711–713 (1999). [CrossRef]

]

ωK=Ω[1Δα2+κ1cos(KΛ)],
(2)

where Ω is the resonant frequency of an individual resonator and Δα and κ 1 are defined as

Δα=d3r[ε(r)ε0(r)]EΩ(r)·EΩ(r)
(3a)
κ1=d3r[ε0(rΛẑ)ε(rΛẑ)]×EΩ(r)·EΩ(rΛẑ).
(3b)

ε (r) is the dielectric coefficient of the CROW and ε0(r) is the dielectric coefficient of an individual resonator. Therefore, the coupling parameter κ 1 represents the overlap of the modes of two neighboring resonators and Δα/2 gives the fractional self frequency shift of ωK.

2. Transfer matrix formalism

We first consider an infinite chain of coupled ring resonators in order to obtain its dispersion relation. Both forward and backward propagating waves exist in an individual resonator, as shown in Fig. 1. We assume the coupling region is sufficiently long compared to λ, so the light circulating in one direction in a resonator is phase-matched to only one of the two degenerate counter-propagating modes of the adjacent resonator. Using the notation of Fig. 1, the coupling between two adjacent rings can thus be described by [16

16. A. Yariv, “Universal relations for coupling of optical power between microresonators and dielectric waveguides,” Electron. Lett. 36, 321–322 (2000). [CrossRef]

]

Fig. 1. An infinitely long chain of coupled ring resonators, with the forward and backward propagating field components labelled.
[bnbn+1]=[tκκ*t*][anan+1],[dndn+1]=[tκκ*t*][cncn+1]
(4)

where t and κ are respectively the dimensionless transmission and coupling coefficients over the coupling length. The matrix is unitary and unimodular so that |t|2+|κ|2=1. Defining a vector with the different field components,

xn=[abcd]n,
(5)

Equation (4) can be rewritten as

xn+1=[P00P]xnxn
(6a)
P=1κ[t11t*]
(6b)

As the field propagates around the ring, it accumulates a phase shift and may be attenuated, so

xn=[0QQ0]xnxn
(7a)
Q=[0eiβRπeiβRπ0]
(7b)

In the above definition, R is the ring radius and β=n(ω)ω/c+, where n(ω) is the frequency dependent effective index and α is the loss (or gain) per unit length in the ring. Combining (6) and (7), we have

xn+1=xn
(8)

Equation (8) is completely general. The matrices ℙ and ℚ can be specified at each frequency to account for any frequency dependence of the effective index, loss, and transmission and coupling coefficients.

3. CROW dispersion relation

The field in one resonator of the CROW as specified by xn is

𝓔(ρ,ϕ)=E(ρ)×{anexp[iβR(πϕ)]+dnexp[iβR(πϕ)]0<ϕ<πbnexp[iβR(π+ϕ)]+cnexp[iβR(π+ϕ)]π<ϕ<0
(9)

where ϕ is the azimuthal angle relative to the propagation direction in the counter-clockwise sense, and ρ is the radial co-ordinate. For a mode of an infinite chain of ring resonators, the fields are periodic at the lattice constant, Λ. So applying Bloch’s theorem,

xn+1=exp(iKΛ)xn,
(10)

where K is the CROW propagation constant. Combining this requirement with (8) leads to

Detexp(iKΛ)U=Det(PQ)2exp(i2KΛ)U=0,
(11)

where U is the identity matrix.

We assume lossless propagation and Im(κ)≫Re(κ) for phase-matched coupling. We recall that at the resonant frequency of an individual resonator, Ω, Ωn(Ω)R/c=m, where m is an integer, and n(Ω) is the effective index at Ω. Therefore, approximating n(Ω)≈n(ωK ), we solve Eq. (11) to obtain

sin(ωKΩmπ)=±Im(κ)cos(KΛ),
(12)

which is the desired dispersion relation for a ring CROW. This relation is exact in the sense that it involves no assumption about the coupling strength.

If we expand Eq. (12) in the parameter Δωmπ/Ω, ΔωωK -Ω, we obtain to first order

ωKΩ=1±κ2cos(KΛ),
(13)

where κ 2≡Im(κ)/(). The two dispersion relations corresponding to the ‘±’ coexist for an infinite structure to allow for both forward and backward wave propagation (i.e. positive and negative group velocities). Physically, for a finite structure without reflection and a unidirectional input as in Fig. 3, only the dispersion relation with the matching group and phase velocities as the input wave will be of significance.

Equation (13) is of a form identical to the tight-binding result in Eq. (2). The correction Δα/2 term does not explicitly appear in Eq. (13) since it is accounted for by Re(κ). From Eq. (13), it follows that for Δωmπ/Ω≪1, it is necessary that |κ|≪1. This condition and the absence of all but the nearest neighbor coupling are thus the validity conditions for the tight-binding approximate result (13).

Fig. 2. The exact and cosine-approximate (i.e. tight-binding-approximate) dispersion relations for m=100 and κ=-0.8i.

Figure 2 shows the dispersion relations for κ=-0.8i and m=100 as calculated using the “exact” form in Eq. (12) and the approximated form in Eq. (13). As ωK /Ω increases, the exact dispersion relation deviates more significantly from the cosine form. For smaller values of κ, the deviation from the cosine dispersion relation is reduced.

The Bloch modes of the CROWare given by the eigenvectors of ℙℚ. At each frequency, there are 4 Bloch modes corresponding to the 4 eigenvalues (i.e. values of K). The eigenvalues are exp(-iK 1Λ)≡ξ1, exp[-i(K 1Λ+π)]≡-ξ 1, exp(-iK 2Λ)≡ξ 2, and exp[-i(K 2Λ+π)]≡-ξ 2. The corresponding (un-normalized) eigenvectors are

q̂ξ1=[ζ+γ1ζ+γ1],q̂ξ1=[(ζ+γ)1ζ+γ1],q̂ξ2=[ζγ1ζγ1],q̂ξ2=[(ζγ)1ζγ1],
(14)

where

γ=12t1+exp(4imπωΩ)+2exp(2imπωΩ)(12t2),
(15a)
ζ=12t[1+exp(2imπωΩ)].
(15b)

The 4 eigenvectors are orthogonal to each other and they represent standing waves in each resonator. In the limit of weak coupling |κ|≪1 and ω≈Ω, such that γ≈|κ|≈0, ζ≈1, and ξ 1=ξ 2, the 4 eigenvectors reduce to 2 degenerate eigenvectors, representing the two different superpositions of the clockwise and counter-clockwise propagating waves in a single resonator:

q̂ξ1=q̂ξ2=[1111],q̂ξ1=q̂ξ2=[1111].
(16)

In the limit of strong coupling and γ+ζ1+κ1κ2t, and ζγ1κ1+κ0. The eigenvectors become

q̂ξ1=[2t12t1],q̂ξ1=[2t12t1],q̂ξ2=[0101],q̂ξ2=[0101].
(17)

We observe that 2 field components are significantly stronger than the other. This corresponds to a wave that “zig-zags” through the resonators without making complete round-trips in each resonator. The asymptotic behavior of the eigenvectors confirms the physical picture that as |κ|→0, the modes of the CROW are essentially the modes of the independent resonators, and as |κ|→1, the microrings no longer act as resonators and the CROW modes are essentially conventional waveguide modes.

4. Finite CROWs and a travelling wave picture

For physical realizations of CROWs, we are interested in finite structures with input and output coupling. These properties cannot be readily dealt with using the tight-binding method, but can be easily incorporated into the transfer matrices. Figure 3 shows a typical implementation of a microring CROW: light is coupled into and out of a set of coupled ring resonators via the input and output waveguides. Assuming that the coupling length between waveguides and the CROW is long compared to λ, then only the travelling wave phase-matched to the input can be excited.

Fig. 3. A CROW consisting of N ring resonators with input and output waveguides.

[ab]n+1=PQ[ab]n,
(18)

where P and Q are defined in Eq. (6) and Eq. (7).

By cascading the transfer matrices, PQ, we obtain an expression for the field components at the output of the CROW after N identical rings:

[aN+1bN+1]=PoutQ(PQ)N1Pin[a0b0][ABCD][a0b0],
(19)

where Pin and Pout describe the coupling between the CROW and the input/output waveguides. For a single input to the waveguide, we set a N+1=0. Therefore, the transfer functions at the “through” and “output” ports as shown in Fig. 3 are

b0a0=ABTthr(ω),
(20a)
bN+1a0=CADBTout(ω).
(20b)

As in microring filter design [17

17. B. E. Little, S. T. Chu, H. A. Haus, J. Foresi, and J.-P. Laine, “Microring resonator channel dropping filter,” IEEE J. Lightwave Technol. 15, 998–1005 (1997). [CrossRef]

], the coupling between the waveguides and the CROW can be selected to maximize the flatness of the transmission response. Therefore, a finite CROW can be designed to mimic an infinite CROW over a bandwidth with a sufficiently flat transmission response.

An advantage of the matrix formalism is that it is valid for chains of any length N, which is essential in analyzing any physical realization of a CROW. From the phase response of the transmission function given by Eq. (20), we can deduce the dispersion relation of the structure. However, we note that the travelling wave is not an eigenmode of the CROW, since the Bloch modes as given by the eigenvectors of ℙℚ are standing waves as in Eq. (14). A travelling wave solution is formed by a superposition, either the sum or difference, of the two Bloch modes with equal group velocities (i.e. ξ1 and -ξ1, or ξ2 and -ξ2). The travelling wave is an eigenvector of (PQ)2, and it is verified that the sense of propagation in the rings alternates between clockwise and counter-clockwise with each operation of ℙℚ, as depicted in Fig. 3. Therefore, taking the phase difference accumulated over two rings to be -2KΛ, where K is the Bloch wave vector, such that the phase difference between the output and the input is approximately -(N-1)KΛ, we can determine the CROW dispersion from the finite structure.

As an example, we compute the dispersion relation of a finite CROW consisting of 20 coupled rings with inter-resonator and waveguide-resonator coupling constants of -0.5i. The rings are lossless, and their radius is 16.4µm. n eff is taken to be constant and equal to 1.5. Figure 4 compares the dispersion relation extrapolated from the finite CROW with the dispersion relation of an infinite CROW as given by Eq. (12). The small amplitude ripples are manifested at the resonance frequencies of the finite structure. In the limit of an infinite number of resonators, the resonance peaks will be infinitesimally close to each other and the ripples will be smoothed out.

5. Pulse propagation

5.1. Semi-infinite case

A semi-infinite microring CROWconsists of an infinitely long CROWcoupled to a single input waveguide as in Fig. (5). The input waveguide ensures that only a pulse of positive (or negative) group velocity propagates through the structure. Assuming that the bandwidth of the input pulse is within the bandwidth of the CROW band such that all of the input light is coupled into the waveguide, the field amplitude b 1 in the first resonator is b 1(ω)=-1/κina 0(ω), where κin is the coupling coefficient between the input waveguide and the first resonator. Since |κin |<1, the intensity of the field inside the CROW is higher than that of the input pulse by 1/|κin |2. This does not violate energy conservation, as the increased intensity is a consequence of the reduced group velocity and hence the spatial compression of the pulse inside the CROW. Using the dispersion relation in Eq. (13), the maximum group velocity in the CROW is

vg,max=κΛΩmπ.
(21)

Defining the “slowing” factor as in Ref. [18

18. A. Melloni and F. Morichetti, “Linear and nonlinear pulse propagation in coupled resonator slow-wave optical structures,” Opt. Quantum Electron. 35, 365–379 (2003). [CrossRef]

] to be

S=cneffvg,max,
(22)

and approximating Λ∝2R, S can be expressed as

S=π2κ.
(23)

Therefore, for κin =κ, the intensity inside the rings is roughly enhanced by (23S)2. This result makes intuitive sense since the only loss mechanism for the otherwise lossless resonators is the inter-resonator coupling. Interestingly, even though the energy velocity of the Bloch modes at Ω corresponds to the group velocity v g,max[11

11. A. Yariv and P. Yeh, Optical waves in crystals: Propagation and control of laser radiation (Wiley, New York, 1984).

], the energy velocity of a wave that is fully coupled into the semi-infinite CROW is proportional to |κ|2. Hence, the intensity enhancement is proportional to the energy velocity reduction rather than the group velocity reduction.

Fig. 4. The exact dispersion relation for an infinite CROW and the dispersion relation as extracted from 20 coupled resonators. The rings have a radius of 16.4µm and the inter-resonator coupling is -0.5i.

Also, in contrast to CROWs, for other coupled resonator structures where there is no feedback between the resonators, such as the SCISSOR [19

19. J. E. Heeber, R. W. Boyd, and Q.-H. Park, “SCISSOR solitons and other novel propagation effects in microresonator-modified waveguides,” J. Opt. Soc. Am. B 19, 722–731 (2002). [CrossRef]

], the slowing factor is approximately proportional to 1/|κ|2 in the case of weak coupling. However, a CROW has the advantage that, even in the presence of loss, it is most transmitting for the frequencies of the CROW band, while a side-coupled resonator is most attenuating near the resonant frequency of the resonator.

To analyze the temporal dynamics a pulse launched into the semi-infinite CROW, we adopt a method of analysis that is analogous to pulse propagation in conventional waveguides such as optical fibers. We shall find in the limit of weak coupling, such that Eq. (13) is a good approximation, there exists a closed-form solution to the evolution of any arbitrary input pulse.

Fig. 5. A semi-infinite CROW.

The electric field at where b 1 is taken, 𝓔(t, z=0), can be expressed as the Fourier integral

𝓔(t,z=0)=banddωb1(ω)exp(iωt),
(24a)
b1(ω)=dt2π𝓔(t,z=0)exp(iωt).
(24b)

At z=NΛ, each frequency component, b 1(ω), acquires a phase shift of NKΛ, so the field is

𝓔(t,z=NΛ)=banddωexp(iωt)dt2π𝓔(t,z=0)exp[i(ωt+K(ω)NΛ)].
(25)

However, K(ω) is given by the dispersion relation of the CROW, Eq. (13). Therefore, instead of integrating over frequency in (25), if we integrate over the half of the Brillouin zone that gives the appropriate group velocity (for example, the right half), we obtain

𝓔(t,z=NΛ)=ΛΩκ2eiΩt0πΛdKsin(KΛ)eiKNΛeiΩκ2cos(KΛ)(tt)
dt2π𝓔(t,z=0)eiΩt.
(26)

Equation (26) can be further simplified by letting x=KΛ, and invoking the Jacobi-Anger expansion [20

20. S. Mookherjea and A. Yariv, “Pulse propagation in a coupled-resonator optical waveguide to all orders of dispersion,” Phys. Rev. E 65, 056 601 (2002). [CrossRef]

],

eiΩκcos(x)(tt)=mcmJm[Ωκ2(tt)]cos(mx)
(27a)
cm={1ifm=02imifm>0,
(27b)

to arrive at

𝓔(t,z=NΛ)=Ωκ2eiΩtmcm0πdxsin(x)cos(mx)eixN
dt2πJm[Ωκ2(tt)]𝓔(t,z=0)eiΩt.
(28)

However, αm,N=0π dx sin(x)cos(mx)e -ixN≠0 only for certain values of m and N:

αm,N={iπ4forN=m1andN=m1iπ4forN=m+1andN=m+12(m2+N21)(m2+N21)24m2N2forN+m=even
(29)

So the equation for the pulse envelope E(t, z), such that 𝓔(t, z)=E(t, z)e iΩt, is given by the convolution integral

E(t,z=NΛ)=κ2Ω2πmcmαm,NdtJm[Ωκ2(tt)]E(t,z=0).
(30)

The Fourier transform of a Bessel function Jn (t) is only defined within |2πf|≤1 [21

21. A. D. Poularikas, The handbook af formulas and tables for signal processing (IEEE Press, New York, 1998). [CrossRef]

], which accounts for the finite bandwidth of the CROW, Ω(1-|κ 2|)≤ω≤Ω(1+|κ 2|).

Equation (30) holds for an arbitrary input pulse and its sole assumption is the cosine-approximate dispersion relation, which is valid for small κ. The nonlinear dispersive nature of the CROW is embodied in the summation over the Bessel functions. Figure 6(a) shows the evolution of a Gaussian input pulse E(t, z=0)=exp(-t 2/T 2) as calculated using Eq. (30). Figure 6(b) shows the numerical results obtained from the transfer matrices. The analytical solution is in excellent agreement with the fully numerical approach. As the pulse propagates, even though the main peak travels at the group velocity, the ripples develop only at the tail end of the pulse.

Fig. 6. Evolution of a 2.4ps (FWHM) Gaussian pulse centered about 1.5µm in a semi-infinite CROW with κ2=0.0016. The fields are normalized to the maximum field amplitude in the first resonator. (a) Theoretical results computed using Eq. (30). (b) Results computed numerically with the transfer matrices using a chain of 100 ring resonators (n eff=1.5, R=16µm).

5.2. Finite case

Pulse propagation through finite CROWs can be easily analyzed using the transfer matrices results of Eq. (20). Since Eq. (20) is specified at each frequency, we simply have to find the product between the transfer functions and the spectral components of the input pulse. The temporal behavior follows naturally from the Fourier transform.

Distortionless propagation through an arbitrary finite CROW can always be achieved if the input pulse is sufficiently narrow-band such that the transmission function of the drop port, as defined in Fig. 3, over the pulse bandwidth is near unity. However, short pulses which become distorted as they propagate in the CROW are also of fundamental interest. For this purpose, we take an example consisting of 10 coupled ring resonators of radius 164.5µm and n eff=1.5. The inter-resonator coupling constant is -0.3i and the coupling between the waveguides and CROW is -0.5i. The transfer characteristics of this structure are shown in Fig. (7). We launch a 30.5ps (FWHM) long pulse centered at 1.55µm into the CROW.

Fig. 7. The transmission characteristics of a 10 ring long CROW. The ring radius is 164.5µm and n eff=1.5. Inter-resonator coupling is -0.3i and the waveguide-CROW coupling is -0.5i. A 30.5 ps (FWHM) long pulse centered at 1.55µm is input into the CROW. (a) Transmittance of at the drop port. The dashed line shows the spectrum of the input pulse. (b) Phase response at the drop port. (c) Transmittance at the through port. (d) Phase response at the through port.

Using the transfer matrices, we can examine how a pulse evolves in the CROWby finding the transfer functions associated with an or bn . Figure 8 shows the evolution of the pulse through the CROW. Even though the output pulse is attenuated compared to the input, the field intensity inside the rings can be greater than the input, as in the case of the semi-infinite CROW. The intensity build-up is verified by a FDTD simulation discussed in Sect. 5.3. The significant increase in the intensity of the input pulse inside the CROWcan be used to enhance the strengths of nonlinear optical interactions. As noted earlier, we can account for loss (or gain) in our model by including an imaginary part to the propagation factor β. We have found the transfer matrices give excellent agreement with experimental results [22

22. G. T. Paloczi, Y. Huang, A. Yariv, and S. Mookherjea, “Polymeric Mach-Zehnder interferometer using serially coupled microresonators,” Opt. Express 11, 2666–2671 (2003). [CrossRef] [PubMed]

].

Another interesting effect is the small amplitude ripple that follows the main peak in each resonator. The ripple is travelling from the end of the CROW back to the start at approximately the group velocity of the forward moving pulse. This is analogous to a reflection from the end of a waveguide, though in the microring CROWdescribed here there are no reflection mechanisms as the coupling is assumed to be perfectly phase-matched. Indeed, Fig. 9 shows that the ripple at the through port is delayed from the the drop port pulse by the travelling time between the input and the drop. Therefore, although the microring CROWis composed of “microscopic,” discrete elements, it possesses certain “macroscopic” properties that mimic conventional waveguides.

Fig. 8. The pulse transmission through the CROW described in Fig. (7). The 0th resonator is the input pulse and the 11th resonator is the output pulse at the drop port.
Fig. 9. The input pulse and the output pulses at the drop and through ports of the the CROW described in Fig. (7). The solid vertical line marks the maximum of the input pulse, and the dashed vertical line marks the maximum of the output pulse at the through port. The peak of the through port pulse occurs about 5ps sooner than the peak of the input.

5.3. FDTD simulations

As a test for the transfer matrix method and a confirmation of the intensity build-up and anomalous dispersion, we use a finite difference time domain (FDTD) simulation to study the pulse propagation through two coupled ring resonators. The waveguides and rings are 0.2µm wide. They are set in air and have an index of refraction of 3.5. The rings have a radius of 5µm, and the wavelength dependent effective index, as extrapolated from a separate FDTD simulation of the waveguides, is n eff=3.617-0.5539λ. The coupling between the rings is -0.32i, and the coupling between the rings and the waveguides is -0.4i. A 2.4ps (FWHM) Gaussian pulse is launched into the system, and the fields at the through port, drop port, and inside the rings are monitored. We compare the transfer matrix method with the FDTD simulation in Fig. 10(a), showing that the approaches are in excellent agreement. The anomalous dispersion at the through port and the increase in intensity in the coupled rings are confirmed by the FDTD simulation and are evident in Fig. 10(b).

Fig. 10. FDTD simulation of 2 coupled ring resonators with input and output waveguides. The radius of the rings is 5µm, and the effective index is n eff=3.617-0.5539λ. The inter-resonator coupling is -0.32i and the waveguide-resonator coupling is -0.4i. The input pulse is a 2.4ps (FWHM) Gaussian centered at 1.55µm. (a) Comparison between the FDTD simulation and the transfer matrix method. Output refers to the drop port. (b) Intensity built-up and anomalous dispersion as confirmed by the FDTD simulation.

6. Comparison with Fabry-Perot resonators

The transfer matrix analysis discussed here is general and can be applied to Bragg stacks and Fabry-Perot resonators as well. A Bragg stack can be regarded as a chain of coupled Fabry-Perot etalons. Bayindir et al. have recently demonstrated a CROW operating near 600nm using coupled Fabry-Perot resonators fabricated in a multi-layer fashion [24

24. M. Bayindir, S. Tanriseven, and E. Ozbay, “Propagation of light through localized coupled-cavity modes in one-dimensional photonic band-gap structures,” Appl. Phys. A 72, 117–119 (2001). [CrossRef]

]. There have also been extensive studies in the linear and nonlinear optical propagation in Bragg stacks [11

11. A. Yariv and P. Yeh, Optical waves in crystals: Propagation and control of laser radiation (Wiley, New York, 1984).

, 25

25. W. Chen and D. L. Mills, “Gap solitons and the nonlinear optical-response of superlattices,” Phys. Rev. Lett. 58, 160–163 (1987). [CrossRef] [PubMed]

, 26

26. C. M. de Sterke and J. E. Sipe, “Envelope-function approach for the electrodynamics of nonlinear periodic structures,” Phys. Rev. A 38, 5149–5165 (1988). [CrossRef] [PubMed]

, 27

27. C. M. de Sterke, D. G. Salinas, and J. E. Sipe, “Coupled-mode theory for light propagation through deep nonlinear gratings,” Phys. Rev. E 54, 1969–1989 (1996). [CrossRef]

, 28

28. B. J. Eggleton, R. E. Slusher, C. M. de Sterke, P. A. Krug, and J. E. Sipe, “Bragg grating solitons,” Phys. Rev. Lett. 76, 1627–1630 (1996). [CrossRef] [PubMed]

, 29

29. D. N. Christodoulides and R. I. Joseph, “Slow Bragg solitons in nonlinear periodic structures,” Phys. Rev. Lett. 62, 1746–1749 (1989). [CrossRef] [PubMed]

]. In Fabry-Perot resonators, the coupling parameter κ corresponds to the transmission coefficient, while t corresponds to the reflection coefficient. They can be calculated from the Fresnel coefficients or from an analysis of a Bragg stack in the case of a Fabry-Perot with Bragg end mirrors [11

11. A. Yariv and P. Yeh, Optical waves in crystals: Propagation and control of laser radiation (Wiley, New York, 1984).

]. The two major differences between the Fabry-Perot resonator and the ring resonator are 1) the former can be fully described by a 2×2 transfer matrix while the latter requires a 4×4 matrix, and 2) the reflection and transmission coefficients may contain a real part. Therefore, in contrast to the ±cos(KΛ) dependence in Eq. (13), there is only one dispersion curve for a Fabry-Perot CROW at a given frequency range and there may also be an additional phase shift in the cosine dependence.

In contrast to ring resonators, the coupling between Fabry-Perot resonators with Bragg end mirrors is controlled by Bragg reflection. Even though the coupling coefficient may be more stringently controlled in these structures, ring resonators remain an attractive option for CROWs in planar integrated optical circuits because they can be fabricated in a single lithographic step. Moreover, recent developments in coupled ring resonators in polymer and semiconductor materials illustrate the potential of using ring resonators as constituent elements in CROWs [22

22. G. T. Paloczi, Y. Huang, A. Yariv, and S. Mookherjea, “Polymeric Mach-Zehnder interferometer using serially coupled microresonators,” Opt. Express 11, 2666–2671 (2003). [CrossRef] [PubMed]

, 30

30. Little Optics press release, “Higher order optical filters using microring resonators” (Little Optics, 2003), http://www.littleoptics.com/hofilter.pdf.

].

7. Conclusion

In summary, the transfer matrix method is used to analyze microring coupled resonators. The transfer matrix and tight-binding approaches yield equivalent dispersion relations in the limit of weak coupling. We also study pulse propagation through semi-infinite and finite CROWs to find intensity enhancement as well as anomalous dispersion. The matrix method can account for finite chains, holds for any coupling strength, applies to travelling waves, and can treat heterogeneous chains consisting of an arbitrary mix of resonators and coupling constants. These features make the transfer matrices versatile for device design and for analyzing experimental results of microring CROWs [22

22. G. T. Paloczi, Y. Huang, A. Yariv, and S. Mookherjea, “Polymeric Mach-Zehnder interferometer using serially coupled microresonators,” Opt. Express 11, 2666–2671 (2003). [CrossRef] [PubMed]

].

Acknowledgments

References and links

1.

A. Yariv, Y. Xu, R. K. Lee, and A. Scherer, “Coupled-resonator optical waveguide: a proposal and analysis,” Opt. Lett. 24, 711–713 (1999). [CrossRef]

2.

Y. Xu, R. K. Lee, and A. Yariv, “Propagation and second-harmonic generation of electromagnetic waves in a coupled-resonator optical waveguide,” J. Opt. Soc. Am. B 77, 387–400 (2000). [CrossRef]

3.

N. Stefanou and A. Modinos, “Impurity bands in photonic insulators,” Phys. Rev. B 57, 12 127–12 133 (1998). [CrossRef]

4.

D. N. Christodoulides and N. K. Efremidis, “Discrete temporal solitons along a chain of nonlinear coupled microcavities embedded in photonic crystals,” Opt. Lett. 27, 568–570 (2002). [CrossRef]

5.

S. Mookherjea and A. Yariv, “Kerr-stabilized super-resonant modes in coupled-resonator optical waveguides,” Phys. Rev. E 66, 046 610 (2002). [CrossRef]

6.

J. E. Heebner and R. W. Boyd, “‘Slow’ and ‘fast’ light in resonator-coupled waveguides,” J. Mod. Opt. 49, 2629–2636 (2002). [CrossRef]

7.

C. K. Madsen, “General IIR optical filter design for WDM applications using all-pass filters,” IEEE J. Lightwave Technol. 18, 860–868 (2000). [CrossRef]

8.

G. Lenz, B. J. Eggleton, C. K. Madsen, and R. E. Slusher, “Optical delay lines based on optical filters,” IEEE J. Quantum Electron. 37, 525–532 (2001). [CrossRef]

9.

B. E. Little, S. T. Chu, W. Pan, D. Ripin, T. Kaneko, Y. Kokubun, and E. Ippen, “Vertically coupled glass microring resonator channel dropping filters,” IEEE Photon. Technol. Lett. 11, 215–217 (1999). [CrossRef]

10.

M. Bayindir, B. Temelkuran, and E. Ozbay, “Tight-binding description of the coupled defect modes in three-dimensional photonic crystals,” Phys. Rev. Lett. 84, 2140–2143 (2000). [CrossRef] [PubMed]

11.

A. Yariv and P. Yeh, Optical waves in crystals: Propagation and control of laser radiation (Wiley, New York, 1984).

12.

K. Oda, N. Takato, and H. Toba, “A wide-FSR waveguide double-ring resonator for optical FDM transmission systems,” IEEE J. of Lightwave Technol. 9, 728–736 (1991). [CrossRef]

13.

R. Orta, P. Savi, R. Tascone, and D. Trinchero, “Synthesis of multiple-ring resonator filters for optical systems,” IEEE Photon. Technol. Lett. 7, 1447–1449 (1995). [CrossRef]

14.

J. V. Hryniewicz, P. P. Absil, B. E. Little, R. A. Wilson, and P.-T. Ho, “Higher order filter response in coupled microring resonators,” IEEE Photon. Technol. Lett. 12, 320–322 (2000). [CrossRef]

15.

A. Melloni, R. Costa, P. Monguzzi, and M. Martinelli, “Ring-resonator filters in silicon oxynitride technology for dense wavelength-division multiplexing systems,” Opt. Lett. 28, 1567–1569 (2003). [CrossRef] [PubMed]

16.

A. Yariv, “Universal relations for coupling of optical power between microresonators and dielectric waveguides,” Electron. Lett. 36, 321–322 (2000). [CrossRef]

17.

B. E. Little, S. T. Chu, H. A. Haus, J. Foresi, and J.-P. Laine, “Microring resonator channel dropping filter,” IEEE J. Lightwave Technol. 15, 998–1005 (1997). [CrossRef]

18.

A. Melloni and F. Morichetti, “Linear and nonlinear pulse propagation in coupled resonator slow-wave optical structures,” Opt. Quantum Electron. 35, 365–379 (2003). [CrossRef]

19.

J. E. Heeber, R. W. Boyd, and Q.-H. Park, “SCISSOR solitons and other novel propagation effects in microresonator-modified waveguides,” J. Opt. Soc. Am. B 19, 722–731 (2002). [CrossRef]

20.

S. Mookherjea and A. Yariv, “Pulse propagation in a coupled-resonator optical waveguide to all orders of dispersion,” Phys. Rev. E 65, 056 601 (2002). [CrossRef]

21.

A. D. Poularikas, The handbook af formulas and tables for signal processing (IEEE Press, New York, 1998). [CrossRef]

22.

G. T. Paloczi, Y. Huang, A. Yariv, and S. Mookherjea, “Polymeric Mach-Zehnder interferometer using serially coupled microresonators,” Opt. Express 11, 2666–2671 (2003). [CrossRef] [PubMed]

23.

S. Longhi, M. Marano, M. Belmonte, and P. Laporta, “Superluminal pulse propagation in linear and nonlinear photonic grating structures,” IEEE. J. Sel. Top. Quantum Electron. 9, 4–16 (2003). [CrossRef]

24.

M. Bayindir, S. Tanriseven, and E. Ozbay, “Propagation of light through localized coupled-cavity modes in one-dimensional photonic band-gap structures,” Appl. Phys. A 72, 117–119 (2001). [CrossRef]

25.

W. Chen and D. L. Mills, “Gap solitons and the nonlinear optical-response of superlattices,” Phys. Rev. Lett. 58, 160–163 (1987). [CrossRef] [PubMed]

26.

C. M. de Sterke and J. E. Sipe, “Envelope-function approach for the electrodynamics of nonlinear periodic structures,” Phys. Rev. A 38, 5149–5165 (1988). [CrossRef] [PubMed]

27.

C. M. de Sterke, D. G. Salinas, and J. E. Sipe, “Coupled-mode theory for light propagation through deep nonlinear gratings,” Phys. Rev. E 54, 1969–1989 (1996). [CrossRef]

28.

B. J. Eggleton, R. E. Slusher, C. M. de Sterke, P. A. Krug, and J. E. Sipe, “Bragg grating solitons,” Phys. Rev. Lett. 76, 1627–1630 (1996). [CrossRef] [PubMed]

29.

D. N. Christodoulides and R. I. Joseph, “Slow Bragg solitons in nonlinear periodic structures,” Phys. Rev. Lett. 62, 1746–1749 (1989). [CrossRef] [PubMed]

30.

Little Optics press release, “Higher order optical filters using microring resonators” (Little Optics, 2003), http://www.littleoptics.com/hofilter.pdf.

OCIS Codes
(230.3120) Optical devices : Integrated optics devices
(230.5750) Optical devices : Resonators
(230.7370) Optical devices : Waveguides

ToC Category:
Research Papers

History
Original Manuscript: October 14, 2003
Revised Manuscript: December 19, 2003
Published: January 12, 2004

Citation
Joyce Poon, Jacob Scheuer, Shayan Mookherjea, George Paloczi, Yanyi Huang, and Amnon Yariv, "Matrix analysis of microring coupled-resonator optical waveguides," Opt. Express 12, 90-103 (2004)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-12-1-90


Sort:  Journal  |  Reset  

References

  1. A. Yariv, Y. Xu, R. K. Lee, and A. Scherer, �??Coupled-resonator optical waveguide: a proposal and analysis,�?? Opt. Lett. 24, 711�??713 (1999). [CrossRef]
  2. Y. Xu, R. K. Lee, and A. Yariv, �??Propagation and second-harmonic generation of electromagnetic waves in a coupled-resonator optical waveguide,�?? J. Opt. Soc. Am. B 77, 387�??400 (2000). [CrossRef]
  3. N. Stefanou and A. Modinos, �??Impurity bands in photonic insulators,�?? Phys. Rev. B 57, 12 127�??12 133 (1998). [CrossRef]
  4. D. N. Christodoulides and N. K. Efremidis, �??Discrete temporal solitons along a chain of nonlinear coupled microcavities embedded in photonic crystals,�?? Opt. Lett. 27, 568�??570 (2002). [CrossRef]
  5. S. Mookherjea and A. Yariv, �??Kerr-stabilized super-resonant modes in coupled-resonator optical waveguides,�?? Phys. Rev. E 66, 046 610 (2002). [CrossRef]
  6. J. E. Heebner and R. W. Boyd, �??�??Slow�?? and �??fast�?? light in resonator-coupled waveguides,�?? J. Mod. Opt. 49, 2629�??2636 (2002). [CrossRef]
  7. C. K. Madsen, �??General IIR optical filter design for WDM applications using all-pass filters,�?? IEEE J. Lightwave Technol. 18, 860�??868 (2000). [CrossRef]
  8. G. Lenz, B. J. Eggleton, C. K. Madsen, and R. E. Slusher, �??Optical delay lines based on optical filters,�?? IEEE J. Quantum Electron. 37, 525�??532 (2001). [CrossRef]
  9. B. E. Little, S. T. Chu, W. Pan, D. Ripin, T. Kaneko, Y. Kokubun, and E. Ippen, �??Vertically coupled glass microring resonator channel dropping filters,�?? IEEE Photon. Technol. Lett. 11, 215�??217 (1999). [CrossRef]
  10. M. Bayindir, B. Temelkuran, and E. Ozbay, �??Tight-binding description of the coupled defect modes in three-dimensional photonic crystals,�?? Phys. Rev. Lett. 84, 2140�??2143 (2000). [CrossRef] [PubMed]
  11. A. Yariv and P. Yeh, Optical waves in crystals: Propagation and control of laser radiation (Wiley, New York, 1984).
  12. K. Oda, N. Takato, and H. Toba, �??A wide-FSR waveguide double-ring resonator for optical FDM transmission systems,�?? IEEE J. of Lightwave Technol. 9, 728�??736 (1991). [CrossRef]
  13. R. Orta, P. Savi, R. Tascone, and D. Trinchero, �??Synthesis of multiple-ring resonator filters for optical systems,�?? IEEE Photon. Technol. Lett. 7, 1447�??1449 (1995). [CrossRef]
  14. J. V. Hryniewicz, P. P. Absil, B. E. Little, R. A. Wilson, and P.-T. Ho, �??Higher order filter response in coupled microring resonators,�?? IEEE Photon. Technol. Lett. 12, 320�??322 (2000). [CrossRef]
  15. A. Melloni, R. Costa, P. Monguzzi, and M. Martinelli, �??Ring-resonator filters in silicon oxynitride technology for dense wavelength-division multiplexing systems,�?? Opt. Lett. 28, 1567�??1569 (2003). [CrossRef] [PubMed]
  16. A. Yariv, �??Universal relations for coupling of optical power between microresonators and dielectric waveguides,�?? Electron. Lett. 36, 321�??322 (2000). [CrossRef]
  17. B. E. Little, S. T. Chu, H. A. Haus, J. Foresi, and J.-P. Laine, �??Microring resonator channel dropping filter,�?? IEEE J. Lightwave Technol. 15, 998�??1005 (1997). [CrossRef]
  18. A. Melloni and F. Morichetti, �??Linear and nonlinear pulse propagation in coupled resonator slow-wave optical structures,�?? Opt. Quantum Electron. 35, 365�??379 (2003). [CrossRef]
  19. J. E. Heeber, R. W. Boyd, and Q.-H. Park, �??SCISSOR solitons and other novel propagation effects in microresonator-modified waveguides,�?? J. Opt. Soc. Am. B 19, 722�??731 (2002). [CrossRef]
  20. S. Mookherjea and A. Yariv, �??Pulse propagation in a coupled-resonator optical waveguide to all orders of dispersion,�?? Phys. Rev. E 65, 056 601 (2002). [CrossRef]
  21. A. D. Poularikas, The handbook af formulas and tables for signal processing (IEEE Press, New York, 1998). [CrossRef]
  22. G. T. Paloczi, Y. Huang, A. Yariv, and S. Mookherjea, �??Polymeric Mach-Zehnder interferometer using serially coupled microresonators,�?? Opt. Express 11, 2666�??2671 (2003). [CrossRef] [PubMed]
  23. S. Longhi, M. Marano, M. Belmonte, and P. Laporta, �??Superluminal pulse propagation in linear and nonlinear photonic grating structures,�?? IEEE. J. Sel. Top. Quantum Electron. 9, 4�??16 (2003). [CrossRef]
  24. M. Bayindir, S. Tanriseven, and E. Ozbay, �??Propagation of light through localized coupled-cavity modes in one-dimensional photonic band-gap structures,�?? Appl. Phys. A 72, 117�??119 (2001). [CrossRef]
  25. W. Chen and D. L. Mills, �??Gap solitons and the nonlinear optical-response of superlattices,�?? Phys. Rev. Lett. 58, 160�??163 (1987). [CrossRef] [PubMed]
  26. C. M. de Sterke and J. E. Sipe, �??Envelope-function approach for the electrodynamics of nonlinear periodic structures,�?? Phys. Rev. A 38, 5149�??5165 (1988). [CrossRef] [PubMed]
  27. C. M. de Sterke, D. G. Salinas, and J. E. Sipe, �??Coupled-mode theory for light propagation through deep nonlinear gratings,�?? Phys. Rev. E 54, 1969�??1989 (1996). [CrossRef]
  28. B. J. Eggleton, R. E. Slusher, C. M. de Sterke, P. A. Krug, and J. E. Sipe, �??Bragg grating solitons,�?? Phys. Rev. Lett. 76, 1627�??1630 (1996). [CrossRef] [PubMed]
  29. D. N. Christodoulides and R. I. Joseph, �??Slow Bragg solitons in nonlinear periodic structures,�?? Phys. Rev. Lett. 62, 1746�??1749 (1989). [CrossRef] [PubMed]
  30. Little Optics press release, �??Higher order optical filters using microring resonators�?? (Little Optics, 2003), <a href="http://www.littleoptics.com/hofilter.pdf.">http://www.littleoptics.com/hofilter.pdf.</a>

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.


« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited