Waveguide couplers with new power splitting ratios made possible by cascading of short multimode interference sections

doi:10.1364/OE.15.001588

« Show journal navigation

Waveguide couplers with new power splitting ratios made possible by cascading of short multimode interference sections

David J. Y. Feng, T. S. Lay, and T. Y. Chang »View Author Affiliations

Optics Express, Vol. 15, Issue 4, pp. 1588-1593 (2007)
http://dx.doi.org/10.1364/OE.15.001588

View Full Text Article

Acrobat PDF (247 KB)

Journal Search
Article Lookup

Article Tools

Citations

Alert me when this article is cited
Export Citation/Save

Full Text
Article Info
References (14)
Cited By
Figures (3)
Metrics

Abstract

We show that it is possible to obtain 2×2 waveguide couplers with new power splitting ratios for cross coupling of 7%, 64%, 80% and 93% by cascading two short MMI sections. These couplers have simple geometry and low loss. They offer valuable new possibilities for designing waveguide power taps, high-Q ring resonators, ladder-structure optical filters, and loop-mirror partial reflectors.

1. Introduction

The 3-dB multi-mode interference (MMI) coupler, in either the 1×2 or the 2×2 form, has been extensively used as power splitter/combiners in Mach-Zehnder interferometers, ring resonators, and other optical waveguide circuits owing to their relative compactness [1

L. B. Soldano and E. C. M. Pennings, “Optical multi-mode interference devices based on self-imaging: principles and applications,” J. Lightwave Technol. 13,615–627 (1995). [CrossRef]

]. Couplers for unequal splitting of power are also needed for other photonic-integration applications such as power taps, high-Q ring resonators [2–4], ladder-structure optical filters [5

S. Matsuo, Y. Yoshikuni, T. Segawa, Y. Ohiso, and H. Okamoto, “A widely tunable optical filter using ladder-type structure,” IEEE Photonics Technol. Lett. 15,1114–1116 (2003). [CrossRef]

], and loop-mirror partial reflectors [6–8]. For 2×2 couplers, the conventional constant-width MMI design is known to be capable of providing only seven different cross-coupling ratios, i.e. K’s (defined as the cross output power divided by the sum of bar and cross output powers), of 1, 0.85, 0.72, 0.5, 0.28, 0.15, and 0 [9

M. Bachmann, P. A. Besse, and H. Melchior:, “Overlapping-image multimode interference couplers with a reduced number of self-images for uniform and nonuniform power splitting,” Appl. Opt. 34,6898–6910 (1995). [CrossRef] [PubMed]

, 10

J. Leuthold and C. H. Joyner, “Multimode interference couplers with tunable power splitting ratios,” J. Lightwave Technol. 19,700–707 (2001). [CrossRef]

]. For low-loss power taps and high-Q ring resonators, a K value that is even lower than 0.15 would be very useful. For loop-mirror partial reflectors [6–8], a K value that is even closer to 0.5 than 0.72 would also be highly desirable. Although it is possible to obtain arbitrary splitting ratios by using a bent MMI coupler [11

Q. Lai, M. Bachmann, W. Hunziker, P. A. Besse, and H. Melchior, “Arbitrary ratio power splitters using angled silica on silicon multimode interference couplers,” Electron. Lett. ,32,1576-1577 (1996). [CrossRef]

,12

D. S. Levy, Y. M. Li, R. Scarmozzino, and R. M. Osgood Jr. “A multimode interference-based variable power splitter in GaAs-AlGaAs,” IEEE Photonics Technol. Lett. 9,1373–1375 (1997) [CrossRef]

], for typical mask resolution and process variability, the final K value obtained after fabrication is more susceptible to error due to the very small bend angle involved and the very high angular sensitivity of the K value. We show in this letter that it is possible to realize new K values of 0.07, 0.64, 0.80, and 0.93 by cascading two short MMI sections. The devices have low loss and do not involve angled or curved geometries.

2. Building-block MMI’s

The MMI sections to be discussed in this letter are all 2×2 waveguide couplers. A schematic diagram of a generic 2×2 MMI coupler is shown in Fig. 1(a). For this group of MMI sections, we set the input and output waveguide pairs of all the devices to have the same center-to-center separation s. This distance s will serve as the common unit of length scaling for all MMI sections. In order to minimize the size of the couplers, s is generally chosen to be as short as the lithographic process to be used would permit. The beat length L_π , between the two lowest guided modes of the same polarization in the MMI waveguide is given according to the standard mode propagation analysis (MPA) by [1

L. B. Soldano and E. C. M. Pennings, “Optical multi-mode interference devices based on self-imaging: principles and applications,” J. Lightwave Technol. 13,615–627 (1995). [CrossRef]

, 9

]

L_{π} = \frac{4 n_{s} W_{e}^{2}}{3 λ_{0}} = \frac{4 n_{s} r^{2} s^{2}}{3 λ_{0}} \equiv A r^{2} s^{2},

(1)

where n_s is the effective index of the slab waveguide from which the MMI waveguide of effective width W_e is formed, λ ₀ is the vacuum wavelength, and r≡W_e /s. Equation (1) holds for both TE and TM polarizations [1

L. B. Soldano and E. C. M. Pennings, “Optical multi-mode interference devices based on self-imaging: principles and applications,” J. Lightwave Technol. 13,615–627 (1995). [CrossRef]

Fig. 1. Schematic diagrams of a generic 2×2 MMI coupler: (a) in plan view showing the notations used in this paper, and (b) in cross-section showing the basic layer structure in the waveguide ridge.

In addition to analyses by MPA, we also carry out numerical simulations by the Beam Propagation Method (BPM) to obtain more realistic assessments. In these studies, the epitaxial slab waveguide on InP substrate is assumed to have a 1.79-μm thick (mostly InAlAs) upper cladding, a 0.45-μm thick core layer containing 3 quantum wells with the absorption edge at 1.48-μm, and a 0.23-μm thick lower cladding (also mostly InAlAs). The basic layer structure in the waveguide ridge is shown schematically in Fig. 1(b). The side walls of the MMI section are assumed to be etched down to the substrate to keep the number of confined modes large. The effective index of the TE fundamental slab mode is 3.2147 at λ₀=1.55 μm. In order to match the low-order modes in 2-D simulation to those in 3-D simulation, the equivalent effective index of the lateral cladding is set to 2.3. The effective width, w, of the access waveguides is assumed to be 2.2 μm and their separation is set to be s = 5 μm.

Four relatively short MMI sections with K ≥ 0.5 have been considered as the building blocks of cascaded MMI’s. For MMI-A (K = 0.5), we have We=s+w=7.2μm=1.44s. For given values of s and w, this is the minimum practical effective width for a low loss MMI [1

L. B. Soldano and E. C. M. Pennings, “Optical multi-mode interference devices based on self-imaging: principles and applications,” J. Lightwave Technol. 13,615–627 (1995). [CrossRef]

]. The device length L_A is (3/2)L _πA=3.1104As ²=215 μm. The actual length in μm is obtained not by Eq. (1) but by a 3-D BPM simulation tool [13

BeamPROP, version 5.1, RSoft Inc., NY (2005)

]. For MMI-B (K=0.5, r = 3), we have W_e = 3s = 15 μm and L_B = (1/2)L _πB=4.5As ²=311 μm. The access waveguides for MMI-B are located exactly at ±W_e /6 from the centerline of the multimode waveguide. For MMI-C (K=0.85, r = 2), we have W_e = 2s = 10 μm and L_C = (3/4)L _πC=3As ²=207 μm. The access waveguides are located exactly at ±W_e /4 from the centerline of the MMI-C waveguide. For MMI-D (K=0.72, r = 2.5), we have W_e = 2.5s = 12.5 μm, L_D = (3/5)L _πD=3.75As ²=259 μm. The geometric centers of the access-waveguide pairs at two ends are offset in opposite directions by ±0.25s from the center of the MMI-D waveguide.

The cross- and bar-transfer functions of these basic MMI’s are deduced from Bachmann et al. [9

]. and summarized in Table I. In this table, β_U is the phase constant for the fundamental mode in the MMI section U, and L_U is the length of that MMI section. The ± signs in the subscripts are significant only for MMI-D in which the access waveguides are not symmetrically placed. The asymmetry leads to two cross optical paths of unequal lengths. According to the 3-D BPM simulations, the total insertion losses of these four basic MMI’s lye in the range of 4 to 5%.

Table I Transfer functions of MMI-A, MMI-B, MMI-C, and MMI-D

Transfer Functions
	Cross-state	Bar-state
MMI -A	$H A_{X} = \frac{1}{\sqrt{2}} \exp (\begin{matrix} - j β_{A} L_{A} + j \frac{π}{4} \end{matrix})$	$H A_{B} = \frac{1}{\sqrt{2}} \exp (\begin{matrix} - j β_{A} L_{A} - j \frac{π}{4} \end{matrix})$
MMI -B	$H B_{X} = \frac{1}{\sqrt{2}} \exp (\begin{matrix} - j β_{B} L_{B} - j \frac{π}{4} \end{matrix})$	$H B_{B} = \frac{1}{\sqrt{2}} \exp (\begin{matrix} - j β_{B} L_{B} + j \frac{π}{4} \end{matrix})$
MMI -C	${HC}_{X} = \cos (\begin{matrix} - \frac{π}{8} \end{matrix}) \exp (\begin{matrix} - j β_{C} L_{C} - j \frac{π}{8} \end{matrix})$	${HC}_{B} = \cos (\begin{matrix} - \frac{3 π}{8} \end{matrix}) \exp (\begin{matrix} - j β_{C} L_{C} + j \frac{3 π}{8} \end{matrix})$
MMI -D	${HD}_{X \pm} = \frac{2}{\sqrt{5}} \cos (\begin{matrix} - \frac{π}{10} \end{matrix}) \exp [- j β_{D} L_{D} - j \frac{(2 \pm 1) π}{10}]$	${HD}_{B} = \frac{2}{\sqrt{5}} \cos (\begin{matrix} - \frac{3 π}{10} \end{matrix}) \exp (\begin{matrix} - j β_{D} L_{D} - j \frac{3 π}{10} \end{matrix})$

3. Cascaded 2×2 MMI couplers

Since all four basic MMI sections are 2×2 couplers with the same s value, they can be simply aligned and cascaded. In this paper, the MMI sections are actually butt jointed without any interconnecting waveguides. When MMI-U and MMI-V are cascaded, the transfer matrix of the cascaded device can be obtained by the following matrix multiplication.

[\begin{matrix} {HUV}_{B \mp} & {HUV}_{X \pm} \\ {HUV}_{X \mp} & {HUV}_{B \pm} \end{matrix}] = [\begin{matrix} {HV}_{B} & {HV}_{X \pm} \\ {HV}_{X \mp} & {HV}_{B} \end{matrix}] [\begin{matrix} {HU}_{B} & {HU}_{X \pm} \\ {HU}_{X \mp} & {HU}_{B} \end{matrix}],

(2)

where HUV_X± and HUV_B± are respectively the transfer functions to the cross output port and the bar output port for the cascaded device. The K-value is given by

KUV = \frac{{∣ {HUV}_{X \pm} ∣}^{2}}{{∣ {HUV}_{X \pm} ∣}^{2} + {∣ {HUV}_{B \pm} ∣}^{2}},

(3)

We have examined the results of all possible combinations using Eq. (2) and verified their consistency with the results of 3-D BPM simulations. The resultant K-values are: KAB=0, KAC=0.15 (0.146), KAD=0.07 (0.075), KBC=0.85 (0.854), KBD=0.93 (0.925), KCD=0.64 (0.641), and KDD=0.80 (0.80). The alternative values given between parentheses are the MPA values. Among these, 0.07, 0.93, 0.64, and 0.80 are new K values that are so far not known to be possible by using a constant-width MMI [9

, 10

J. Leuthold and C. H. Joyner, “Multimode interference couplers with tunable power splitting ratios,” J. Lightwave Technol. 19,700–707 (2001). [CrossRef]

]. The optical field maps of these MMI’s obtained by 2-D BPM simulation are shown in Fig. 2. Additional information on these four cases are given below.

KAD=0.07

The total length of this cascaded device is [(3/2)L _πA+(3/5)L _πD]=[3.1104As ²+3.75As ²] =474 μm. The MPA and matrix multiplication give |HAD_X|²=(0.075) and |HAD_B|²=(0.925). The total insertion loss is 5% according to the 3-D simulation.

KBD=0.93

This cascaded device is complementary to the 7% A+D device. The reason behind this complementary relation is that the transfer functions of MMI-B are complementary in phase changes to those of MMI-A as is evident in Table I. The total length is [(1/2)L _πB+(3/5)L _πD]=[4.5As ²+3.75As ²]=570 μm. The total insertion loss is 6%.

KCD=0.64

The total length is [(3/4)L _πC+(3/5)L _πD]=[3As ²+3.75As ²]=466 μm. The MPA and matrix multiplication give predictions of |HCD_X|²=(0.643) and |HCD_B|²=(0.357). The total insertion loss is 5%.

KDD=0.80

In this case, two MMI-D’s are cascaded with the double images shifting successively in the same direction. The total length is (6/5)L _πD=7.5As ²=518 μm. The MPA and matrix multiplication give predictions of |HDD_X|²=(0.80) and |HDD_B|²=(0.20). The total insertion loss is 6%. If the second MMI-D section is flipped over horizontally, the result is a constant-width MMI-2D of length (6/5)L _πD with K2D=0.72 (0.724).

Fig. 2. The simulated 2-D field maps for four different cascaded MMI couplers with K=0.07 (a), 0.93 (b), 0.64 (c), and 0.80 (d). The values for total transmittance and K are obtained by 2-D BPM simulations.

Since the cascaded devices are butt jointed without interconnecting waveguides, the total insertion losses of the cascaded devices are much less than the sum of the total insertion losses of the constituent MMI sections. For instance, the insertion loss of the straight-edged MMI-2D is essentially the same as that of the shorter MMI-D, and the insertion loss of the step-edged MMI-DD is only 1% higher than that of MMI-2D. For all four devices, cascading increases the total insertion loss by less than 2% over those of the single sections alone. This means that, at the junction plane, the optical field is so well focused around the double images that hardly any optical field is intercepted by the abrupt change of width at the junction. At the junction plane, the near absence of optical field away from the focused images also makes the devices relatively insensitive to process induced distortions of the steps, provided that the location of the junction plane does not deviate from the actual image plane by much more than the equivalent Rayleigh length of the images.

4. Wavelength sensitivity of cascaded couplers

When the input wavelength deviates from the design wavelength of 1.55 μm, the total transmittance decreases because L _πU and L _πV are inversely proportional to the wavelength and the device length no longer matches the correct sum of their specific multiples. The results of 3-D simulation using another simulation tool based on modal expansion [14

FimmProp, version 4.3, Photon Design, Oxford, UK (2004)

] are plotted in Fig. 3(a). This simulation tool tracks the wavelength dependence of the mode profiles exactly. The calculated values of the total transmittance are, however, about 2% lower than those obtained by the other simulation tool [13

BeamPROP, version 5.1, RSoft Inc., NY (2005)

]. This is partly because only the part of the output optical field that matches exactly the fundamental mode in the output waveguide is taken as the output. Otherwise, the differences are not significant. The bandwidth at 85% total transmittance (i.e. -1 dB from the maximum) is approximately inversely proportional to the device length. This is a general property of all MMI’s. Each of the MMI images mimics the Gaussian-like input field and has the same effective Rayleigh length. The bandwidth is proportional to the fractional change of the wavelength that moves the image by an effective Rayleigh length. Each of the MMI images moves in proportion to the device length multiplied by the fractional change of the wavelength. Additional simulations show that there is no discernible difference between the 1-dB bandwidths of MMI-DD and MMI-2D. These are two MMI’s of the same width and length, one with stepped edges and the other with perfectly straight edges. The result shows that the steps on the edges do not cause noticeable additional loss when the wavelength deviates from the design value.

The dependences of the K values on the wavelength deviation are plotted in Fig. 3(b). Within the 1-dB bandwidth of each of the four devices, the variation of the K value is found to be less than 0.02.

Fig. 3. (a) The dependence of the total transmittances on the deviation of the input wavelength from the design wavelength. (b) The corresponding variation of the K values.

5. Applications

The K values of 7% and 93% are obviously very useful when one wishes to tap off only a small amount of optical power from the waveguide. The K values of 93% and 80% would also be useful in the construction of ladder-structure optical filters [5

S. Matsuo, Y. Yoshikuni, T. Segawa, Y. Ohiso, and H. Okamoto, “A widely tunable optical filter using ladder-type structure,” IEEE Photonics Technol. Lett. 15,1114–1116 (2003). [CrossRef]

]. The value of 7% is very attractive for constructing a high-Q ring resonator. The filtering property of a symmetric single-ring resonator is very similar to that of a symmetric Fabry-Perot resonator. The coupling factor K is equivalent to 1-R in the Fabry-Perot resonator where R is the mirror reflectance. The finesse of a lossless symmetric ring resonator is, therefore, given by

F = \frac{π \sqrt{1 - K}}{K} .

(4)

For previously available K values of 0.15, and 0.5, the (lossless) finesses are respectively 19.3, and 4.44. By using the available new K value of 0.07, one can obtain a much higher (lossless) finesse of 43.3.

A partially reflecting loop mirror can be constructed by using a 2×2 coupler [6–8

V. M. Menon, W. Tong, C. Li, F. Xia, I. Glesk, P. R. Prucnal, and S. R. Forrest, “All-optical wavelength conversion using a regrowth-free monolithically integrated Sagnac interferometer,” IEEE Photonics Technol. Lett. 15,254–256 (2003). [CrossRef]

]. The access waveguides at the far end would be connected together through a waveguide loop. At the near end, one access waveguide would serve as the input and the other as the output. For a lossless loop mirror, the input reflectance is given by

R = {∣ 2 H_{B} H_{X \pm} ∣}^{2} .

(5)

For K=|H_X±|²=0.5, 0.64, 0.72 (and 0.28), 0.80, 0.85 (and 0.15), and 0.93 (and 0.07), we get respectively R=1, 0.92, 0.81, 0.64, 0.51, and 0.27. The transmittance to the output port is 1-R. The newly available K value of 64% makes it possible to obtain a high reflectance of 92% with 8% output coupling.

6. Conclusions

In summary, we have obtained new power splitting ratios of K=0.07, 0.64, 0.80, and 0.93 by cascading two short MMI sections. These couplers have low insertion loss and simple geometries requiring no angled or curved patterns. They provide new opportunities to improve the performance of waveguide power taps, high-Q ring resonators, ladder-structure optical filters, and loop-mirror partial reflectors.

Acknowledgements

This work was supported by the Ministry of Education, Taiwan under the Aim for the Top University Plan.

References and links

1.	L. B. Soldano and E. C. M. Pennings, “Optical multi-mode interference devices based on self-imaging: principles and applications,” J. Lightwave Technol. 13,615–627 (1995). [CrossRef]
2.	D. G. Rabus, M. Hamacher, U. Troppenz, and H. Heidrich, “Optical filters based on ring resonators with integrated semiconductor optical amplifiers in GaInAsP-InP,” IEEE J. Sel. Topics Quantum Electron. 8,1405–1411 (2002). [CrossRef]
3.	G. Lenz and C. K. Madsen, “General optical all-pass filter structures for dispersion control in WDM systems,” J. Lightwave Technol. 17,1248–1254 (1999). [CrossRef]
4.	S. Suzuki, K. Oda, and Y. Hibino, “Integrated-optic double-ring resonators with a wide free spectral range of 100 GHz,” J. Lightwave Technol. 13,1766–1770 (1995) [CrossRef]
5.	S. Matsuo, Y. Yoshikuni, T. Segawa, Y. Ohiso, and H. Okamoto, “A widely tunable optical filter using ladder-type structure,” IEEE Photonics Technol. Lett. 15,1114–1116 (2003). [CrossRef]
6.	V. M. Menon, W. Tong, C. Li, F. Xia, I. Glesk, P. R. Prucnal, and S. R. Forrest, “All-optical wavelength conversion using a regrowth-free monolithically integrated Sagnac interferometer,” IEEE Photonics Technol. Lett. 15,254–256 (2003). [CrossRef]
7.	M. G. Kane, I. Glesk, J. P. Sokoloff, and P. R. Prucnal, “Asymmetric optical loop mirror: analysis of an alloptical switch,” Appl. Opt. 33,6833–6842 (1994). [CrossRef] [PubMed]
8.	D. B. Mortimore, “Fiber loop reflectors,” J. Lightwave Technol. 6,1217–1224 (1988). [CrossRef]
9.	M. Bachmann, P. A. Besse, and H. Melchior:, “Overlapping-image multimode interference couplers with a reduced number of self-images for uniform and nonuniform power splitting,” Appl. Opt. 34,6898–6910 (1995). [CrossRef] [PubMed]
10.	J. Leuthold and C. H. Joyner, “Multimode interference couplers with tunable power splitting ratios,” J. Lightwave Technol. 19,700–707 (2001). [CrossRef]
11.	Q. Lai, M. Bachmann, W. Hunziker, P. A. Besse, and H. Melchior, “Arbitrary ratio power splitters using angled silica on silicon multimode interference couplers,” Electron. Lett. ,32,1576-1577 (1996). [CrossRef]
12.	D. S. Levy, Y. M. Li, R. Scarmozzino, and R. M. Osgood Jr. “A multimode interference-based variable power splitter in GaAs-AlGaAs,” IEEE Photonics Technol. Lett. 9,1373–1375 (1997) [CrossRef]
13.	BeamPROP, version 5.1, RSoft Inc., NY (2005)
14.	FimmProp, version 4.3, Photon Design, Oxford, UK (2004)

OCIS Codes
(130.2790) Integrated optics : Guided waves
(130.3120) Integrated optics : Integrated optics devices

ToC Category:
Integrated Optics

History
Original Manuscript: January 3, 2007
Revised Manuscript: January 25, 2007
Manuscript Accepted: January 26, 2007
Published: February 19, 2007

Citation
David J. Y. Feng, T. S. Lay, and T. Y. Chang, "Waveguide couplers with new power splitting ratios made possible by cascading of short multimode interference sections," Opt. Express 15, 1588-1593 (2007)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-15-4-1588

Sort: Author | Year | Journal | Reset

References

L. B. Soldano and E. C. M. Pennings, "Optical multi-mode interference devices based on self-imaging: principles and applications," J. Lightwave Technol. 13, 615-627 (1995). [CrossRef]
D. G. Rabus, M. Hamacher, U. Troppenz, and H. Heidrich, "Optical filters based on ring resonators with integrated semiconductor optical amplifiers in GaInAsP-InP," IEEE J. Sel. Topics Quantum Electron. 8, 1405-1411 (2002). [CrossRef]
G. Lenz and C. K. Madsen, "General optical all-pass filter structures for dispersion control in WDM systems," J. Lightwave Technol. 17, 1248-1254 (1999). [CrossRef]
S. Suzuki, K. Oda, and Y. Hibino, "Integrated-optic double-ring resonators with a wide free spectral range of 100 GHz," J. Lightwave Technol. 13, 1766-1770 (1995) [CrossRef]
S. Matsuo, Y. Yoshikuni, T. Segawa, Y. Ohiso, and H. Okamoto, "A widely tunable optical filter using ladder-type structure," IEEE Photonics Technol. Lett. 15, 1114-1116 (2003). [CrossRef]
V. M. Menon, W. Tong, C. Li, F. Xia, I. Glesk, P. R. Prucnal, and S. R. Forrest, "All-optical wavelength conversion using a regrowth-free monolithically integrated Sagnac interferometer," IEEE Photonics Technol. Lett. 15, 254-256 (2003). [CrossRef]
M. G. Kane, I. Glesk, J. P. Sokoloff, and P. R. Prucnal, "Asymmetric optical loop mirror: analysis of an all-optical switch," Appl. Opt. 33, 6833-6842 (1994). [CrossRef] [PubMed]
D. B. Mortimore, "Fiber loop reflectors," J. Lightwave Technol. 6, 1217-1224 (1988). [CrossRef]
M. Bachmann, P. A. Besse, and H. Melchior:, "Overlapping-image multimode interference couplers with a reduced number of self-images for uniform and nonuniform power splitting," Appl. Opt. 34, 6898-6910 (1995). [CrossRef] [PubMed]
J. Leuthold C. H. Joyner, "Multimode interference couplers with tunable power splitting ratios," J. Lightwave Technol. 19, 700-707 (2001). [CrossRef]
Q. Lai, M. Bachmann, W. Hunziker, P. A. Besse, and H. Melchior, "Arbitrary ratio power splitters using angled silica on silicon multimode interference couplers," Electron. Lett., 32, 1576-1577 (1996). [CrossRef]
D. S. Levy, Y. M. Li, R. Scarmozzino, and R. M. OsgoodJr. "A multimode interference-based variable power splitter in GaAs-AlGaAs," IEEE Photonics Technol. Lett. 9, 1373-1375 (1997) [CrossRef]
BeamPROP, version 5.1, RSoft Inc., NY (2005).
FimmProp, version 4.3, Photon Design, Oxford, UK (2004).

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.

Click here to see a list of articles that cite this paper