Coaxial periodic optical waveguide - OSA Publishing

3 downloads 0 Views 591KB Size Report
Tetsuya Kawanishi and Masayuki Izutsu. Communications Research Laboratory, Ministry of Posts and. Telecommunications, 4-2-1 Nukui-kita, Koganei, Tokyo ...
Coaxial periodic optical waveguide Tetsuya Kawanishi and Masayuki Izutsu Communications Research Laboratory, Ministry of Posts and Telecommunications, 4-2-1 Nukui-kita, Koganei, Tokyo 184-8795, Japan [email protected]

Abstract: Guided modes in a dielectric waveguide structure with a coaxial periodic multi-layer are investigated by using a matrix formula with Bessel functions. We show that guided modes exist in the structure, and that the field is confined in the core which consists of the optically thinner medium. The dispersion curves are discontinuous, so that the modes can exist only in particular wavelength bands corresponding to the stop bands of the periodic structure of the clad. It is possible that the waveguide structure can be applied to filters or optical fibers to reduce nonlinear effects. c 2000 Optical Society of America  OCIS codes: (230.7370) Optical devices, Waveguides

References and links 1. J. C. Knight, T. A. Birks, P. St. J. Russell and D. M. Atkin, “All-silica single-mode optical fiber with photonic crystal cladding,” Opt. Lett. 21, 1547-1549 (1996). 2. J. C. Knight, T. A. Birks, P. St. J. Russell and D. M. Atkin, “All-silica single-mode optical fiber with photonic crystal cladding: errata,” Opt. Lett. 22, 484-485 (1997). 3. T. A. Birks, J. C. Knight and P. St. J. Russell, “Endlessly single-mode photonic crystal fiber,” Opt. Lett. 22, 961-963 (1997). 4. J. C. Knight, T. A. Birks, P. St. J. Russell and J. P. de Sandro, “Properties of photonic crystal fiber and the effective index model,” J. Opt. Soc. Am. 15, 748-752 (1998). 5. J. C. Knight, T. A. Birks, R. F. Cregan, P. St. J. Russell and J. P. de Sandro, “Large mode area photonic crystal fibre,” Electron. Lett. 34, 1347-1348 (1998). 6. J. Broeng, S. E. Barkou, A. Bjarklev, J. C. Knight, T. A. Birks and P. St. J. Russell, “Highly increased photonic band gaps in silica/air structures,” Opt. Commun. 156, 240-244 (1998). 7. J. C. Knight, J. Broeng, T. A. Birks and P. St. J. Russell, “Photonic Band Gap Guidance in Optical Fibers,” Science 282, 1476-1478 (1998). 8. S. E. Barkou, J. Broeng and A. Bjarklev, “Silica-air photonic crystal fiber design that permits waveguiding by a true photonic bandgap effect,” Opt. Lett. 24, 46-48 (1999). 9. A. Ferrando, E. Silverstre, J. J. Miret, P. And´es and M. V. And´es “Full-vector analysis of a realistic photonic crystal fiber,” Opt. Lett. 24, 276-278 (1999). 10. Z. Xue-Heng, “Theory of two-dimensional ‘Fingerprint’ resonators,” Electron.Lett. 25, 1311-1312 (1989). 11. M. Toda, “Single-mode behavior of a circular grating for potential disk-shaped DFB lasers,” IEEE J. Quantum Electron. 26, 473-481 (1990). 12. T. Erdogan and D. G. Hall, “Circularly symmetric distribution feedback semiconductor laser: An analysis,” J.Appl.Phys. 68, 1435-1444 (1990). 13. X. H. Zheng and S. Lacroix, “Mode coupling in circular-cylindrical system and its application to fingerprint resonators,” J.Lightwave Technol. 8, 1509-1516 (1990). 14. C. Wu, T. Makino, J. Glinski, R. Maciejko and S. I. Najafi, “Self-consistent coupled-wave theory for circular gratings on planar dielectric waveguides,” J.Lightwave Technol. 8, 1264-1277 (1991). 15. C. Wu, T. Makino, R. Maciejko, S. I. Najafi and M. Svilans, “Simplified coupled-wave equations for cylindrical waves in circular grating planar waveguides,” J.Lightwave Technol. 10, 1575-1589 (1992). 16. S. Noda, T. Ishikawa, M. Imada and A. Sasaki, “Surface-emitting device with embedded circular grating coupler for possible application on optoelectronic integrated devices,” IEEE J.Photon.Technol.Lett 7, 1397-1399 (1995). 17. A. Yariv and A. Gover, “Equivalence of the coupled-mode and Floquet-Bloch formalisms in periodic optical waveguides,” Appl. Phys. Lett. 26, 537-539 (1975).

#22839 - $15.00 US

(C) 2000 OSA

Received June 01, 2000; Revised June 14, 2000

3 July 2000 / Vol. 7, No. 1 / OPTICS EXPRESS 10

18. A. Y. Cho, A. Yariv and P. Yeh, “Observation of confined propagation in Bragg waveguides,” Appl. Phys. Lett. 30, 471-472 (1977). 19. P. Yeh, A. Yariv and A. Y. Cho, “Optical surface waves in periodic layered media,” Appl. Phys. Lett. 32, 104-105 (1978). 20. W. Ng, P. Yeh, P. C. Chen and A. Yariv, “Optical surface waves in periodic layered medium grown by liquid phase epitaxy,” Appl. Phys. Lett. 32, 370-371 (1978).

1

Introduction

Recently, photonic devices using periodic structures, such as photonic crystals, have received a great deal of attention due to their potential for use to perform new functions. Photonic crystal fibers (PCFs) of which clads are two-dimensional photonic crystal structures have been investigated, experimentally and theoretically [1, 2, 3, 4, 5, 6, 7, 8, 9]. We can categorize PCFs into two types. In the first, refractive index of the core is larger than that of the clad. A guided wave mode, which occurs as a result of the difference of the refractive indices and is similar to that of a conventional optical fiber, is called a conventional mode (C-Mode) in this paper. In the second type, the refractive index of the core is smaller than that of the clad, thus the propagating wave can not be confined in the core, as a result of the difference of the refractive indices. However, when the propagating wave is in the photonic band gap of the photonic crystal, a guided mode may exist, because the wave is reflected into the core due to Bragg reflection. We call this type a Bragg mode (B-Mode). In this type PCF, an air hole can be used as the core to reduce nonlinear effects. In this paper, we investigate a cylindrical waveguide structure having a coaxial periodic structure as its clad. We call it a coaxial periodic optical waveguide (CPOW). Its periodic structure consists of optically denser and thinner media, and it is periodic with respect to the radius. The reflective index of the core is smaller than the average in the clad, therefore, there is not a C-Mode, but B-Modes may exist, and we can expect that the dispersion curves of these modes to be discontinuous. Previously, optical resonators and surface-emitting lasers using circular grating were studied [10, 11, 12, 13, 14, 15, 16]. Their structures are similar to the CPOW and they are also similarly affected by Bragg reflection. But, those structures have slab waveguide structures, and their modes are considered as modified slab waveguide modes by using perturbative methods such as the mode matching theory. Thus, the modes in those structures are basically C-modes, while in this paper the B-modes propagating along the coaxial axis will be discussed. When the EM field in a coaxial structure depends on azimuth angle, the EM field can not be resolved into TE and TM components, thus we should consider them to be hybrid modes as we do in conventional optical fibers. In section 2, we discuss an optical waveguide with one-dimensional photonic bandgap structures to show the basic concept of the B-mode. In section 3, we give an formulation of an EM field in a CPOW by using Bessel functions, and give formal expressions for guided modes. In section 4, we show the numerical results obtained by using the expressions given in section 3, and the finite element method. 2

Waveguide structure with one-dimensional periodic structures

Consider a waveguide structure with one-dimensional periodic structures as shown in Fig. 1. Guided waves are reflected by the periodic structures, and propagate in the gap between the periodic structures [17, 18, 19, 20]. We assume that the periodic structure consists of optically thinner layers (thickness a, refractive index nI = 2I ) and optically denser layers (thickness b, refractive index nII = 2II ), and that n2II is greater than n2I . The wavenumbers of the Bloch wave for the TE- and TM-polarized waves in the periodic

#22839 - $15.00 US

(C) 2000 OSA

Received June 01, 2000; Revised June 14, 2000

3 July 2000 / Vol. 7, No. 1 / OPTICS EXPRESS 11

structure, denoted by kTE and kTM , can be expressed by n21 cos2 θ1 + n22 cos2 θ2 sin(n1 cos θ1 ka) sin(n2 cos θ2 kb) 2n1 n2 cos θ1 cos θ2 + cos(n1 cos θ1 ka) cos(n2 cos θ2 kb) n2 cos2 θ1 + n21 cos2 θ2 sin(n1 cos θ1 ka) sin(n2 cos θ2 kb) cos kTM (a + b) = − 2 2n1 n2 cos θ1 cos θ2 + cos(n1 cos θ1 ka) cos(n2 cos θ2 kb), cos kTE (a + b) = −

(1)

(2)

where k is the wavenumber in free space, and θ1 and θ2 are as defined in Fig. 2. When | cos kTE (a + b)| > 1,

(3)

kTE is a pure imaginary number, and corresponds to an evanescent wave. Similarly, when (4) | cos kTM (a + b)| > 1, kTM is a pure imaginary number. Eqs. (3) and (4) define the stopbands of the periodic structure, in which the wave propagating in the periodic structure is perfectly reflected. In addition, when Eq. (4) is satisfied, | cos kTE (a+b)| is always greater than | cos kTM (a+ b)|. Thus, when a TM-polarized wave is in the stopband, a TE-polarized one is also. The decay constant of a TE-polarized wave is greater than that of a TM-polarized wave. Assuming that the waves in an optical waveguide having a coaxial periodic structure have characteristics similar to a waveguide having one-dimensional periodic structures, Eqs. (3) and (4) can express an approximate necessary condition for guided modes to exist. 3

Electromagnetic field in a coaxial periodic structure

Consider a waveguide having a cylindrical structure like the one shown in Fig. 3. The core of the waveguide is an air hole, whose radius is denoted by r0 , at the center of the cross-section. The clad consists of a coaxial periodic structure, where the radius of the i-th boundary is defined by  n(b + a) + c i = 2n ri = (n = 0, 1, 2, 3, · · ·). (5) n(b + a) + b + c i = 2n + 1 As shown in Fig. 4, the relative permittivity for the ith layer, which corresponds to the region ri−1 < r < ri , is defined by  I i = 2n i = (n = 0, 1, 2, 3, · · ·). (6) II i = 2n + 1 We use the coodinate system in which the z-axis corresponds to the axis of the coaxial structure. The xy components are denoted by a polar form x = r cos φ and y = r sin φ. The z-direction components of the electric field Ez (r) and the magnetic field Hz (r) for the i-th layer are assumed to be expressed by using the first kind Bessel function Jm and the second kind Ym , as follows: Ez (r)

=

{Ai Jm (qi r) + Bi Ym (qi r)} sin(mφ + θm )

(7)

Hz (r)

=

{Ci Jm (qi r) + Di Ym (qi r)} cos(mφ + θm ),

(8)

where qi2 = ω 2 i 0 µ0 −β 2 , m is a non-negative integer number, and θm is an real number. 0 and µ0 denote the permittivity and the permeability of the vacuum.

#22839 - $15.00 US

(C) 2000 OSA

Received June 01, 2000; Revised June 14, 2000

3 July 2000 / Vol. 7, No. 1 / OPTICS EXPRESS 12

The factor exp(jωt − jβz), which expresses the wave propagation along the z-axis, is suppressed for brevity. The field in the core (0 < r < r0 ) is denoted by Eqs. (7) and (8) of i = 0, and should be finite even at r = 0, so that the coefficients of the second kind Bessel functions, B0 and D0 , should be zero. The z- and φ-direction components are expressed in the following matrix form:   Ez (r)  Hz (r)    (9)  Eφ (r)  = U · T (r, i ) · ui Hφ (r)   sin(mφ + θm ) 0 0 0   0 cos(mφ + θm ) 0 0  (10) U ≡   0 0 cos(mφ + θm ) 0 0 0 0 sin(mφ + θm )   Jm (qi r) Ym (qi r) 0 0   0 0 Jm (qi r) Ym (qi r)     jωµ0 Jm (qi r) jωµ0 Ym (qi r)  (11) (qi r) jβmYm (qi r) T (r, i ) ≡  − jβmJm − + + qi qi   q2 r q2 r −

i

 jω i Jm (qi r) qi



t

i

 jω i Ym (qi r) qi

(qi r) + jβmJqm 2r i

(qi r) + jβmYqm 2r

ui = [Ai , Bi , Ci , Di ] .

i

(12)

The boundary conditions for tangential components of the electric and magnetic fields at r = ri are expressed by T (ri , i ) · ui = T (ri , i+1 ) · ui+1 .

(13)

In the case of m = 0, where the field is independent of φ, we can discuss the field in the waveguide by using the concept of TE- and TM-polarized waves, since there is no conection between the coefficient couple (Ai , Bi ) and the couple (Ci , Di ). The TMpolarized wave, in which Ez = 0, can be expressed by Ai and Bi . The TE-polarized wave, in which Hz = 0, can be expressed by Ci and Di . However, when m = 0, TE- or TM-polarized waves do not exist alone, thus the guided mode should be treated as a hybrid mode. By using Eq. (13) iteratively with putting initial values for the coefficients A0 and C0 in the core, we can obtain, the EM field in the waveguide. However, because of numerical calculation error, it is difficult to obtain a solution which decreases exponentially with respect to the radius r. In general, there exists an exponentially increasing solution when an exponentially decreasing solution exists in a periodic structure. Numerical calculation error contains both decreasing and increasing solutions, so that the solution which increases along the forward direction of the calculation steps becomes dominant after several steps, and this solution becomes independent of initial values. Thus, we can obtain decreasing solutions, with respect to the radius r, by using Eq.(13) iteratively, toward the core from a layer whose radius is much greater than the wavelength of the propagating wave. 4

Numerical calculation

By using u0 =

n 

Ri · un+1 ,

Ri ≡ T −1 (ri , i ) · T (ri , i+1 ),

(14)

i=0

we get the coefficients in the core A0 , B0 , C0 , and D0 , where arbitrary values are put on An , Bn , Cn , and Dn (n = 100). Solutions for the guided modes are obtained by searching

#22839 - $15.00 US

(C) 2000 OSA

Received June 01, 2000; Revised June 14, 2000

3 July 2000 / Vol. 7, No. 1 / OPTICS EXPRESS 13

the exponentially decreasing ones, with respect to r, for which coefficients B0 and D0 are nearly equal to zero, with varying parameters k = ω0 µ0 and β/k. The solutions are normalized to be A0 = 1 or C0 = 1. The thickness and the refractive index of the layers in the periodic structure are assumed as follows: k0 a = k0 b = 1.0 and II /I = 2.0, where k0 is an arbitrary unit wavenumber of the structure. We can obtain dispersion curves by plotting the points where guided modes exist. Figs. 5 and 6, respectively, show the dispersion curves for k0 c = 1 and k0 c = 5. Here, we consider the case of m = 0, so that the TE- and TM-polarized modes can be obtained separately. The solution which satisfies B0 = 0 is TM-polarized, and for D0 = 0 it is TE-polarized. Here the x- and y-axes show, respectively, k/k0 and β/k. Fig. 7 shows the points which correspond to guided modes for k0 r = 100. In this case, the radius of the core is much larger than the wavelength, so that the region where the points exist approximately corresponds to the stopband of the coaxial periodic structure. Actually, the points shown in Figs. 5 and 6 are in the stopband shown in Fig. 7. For comparison, the stopband of the waveguide with one-dimensional periodic structures, calculated by using Eqs. (1) and (2), is shown in Fig. 8. Except for the region of β/k ∼ 1, the stopbands shown in Figs. 7 and 8 are analogous. When β/k ∼ 1, the wavelength with respect to r is large, so that the difference between a plane and cylindrical waves becomes significant. In addition, we can see that the stopband of the TM-polarized wave is contained by that of the TE-polarized wave. The fields for the cases of k = 1.2k0 and 3.0k0 are shown in Figs. 9 and 10, where k0 c = 5.0 (see Fig. 6). To clarify the decay due to exponential decay, the y-axis is assigned to electric or magnetic field multiplied by the square root of r, so that we can compensate for the intrinsic decay of the Bessel functions. The hatched region indicates the optically denser (2n + 1-th) layers. The derivatives of the fields with respect to r are discontinuous at the boundaries of the layers. The discontinuity becomes large when β/k is large, as we can see from Eq. (11). In the case of m = 0, we should deal with a couple of exponentially increasing solutions whose increasing factors are different. If the difference is not small, the solution for which factor is larger becomes dominant, so that it is difficult to obtain linear combinations of these two increasing solutions. To discuss wave field distribution and propagation loss of guided modes in a CPOW, EM wave field was calculated by using the finite element method (FEM). We consider a CPOW model structure which has nine layers with parameters k0 a = k0 b = 1.0, k0 c = 5.0, 2π/k0 = 1.5µm and l = 7.5µm, where l denotes the length of the model structure along the z-axis. In the FEM, EM wave can be excited by distinct guided modes of the ports which are defined in model structures. When the mode in the port has rotational symmetry order of m, the excited EM wave field in the CPOW also has the same symmetry. Thus, we can specify the structure of excited EM wave field by choosing the input wave mode. To carry out numerical calculations, the ports should be so simple that the input wave modes can be specified easily. We put circular waveguides at each ends of the CPOW as input and output ports, where the radius of the circular waveguide rf and the length lf are defined as rf = 1µm and lf = 2.5µm. The outer surface of the CPOW is assumed to be radation boundary, and that of the port is perfect conductor. There are guided modes TEmn and TMmn in the ports. Numerical results were obtained by using HP HFSS 5.4. Fig. 11 shows the EM wave intensity distribution in the model excited by the TE01 mode of k = 1.2k0 in the input port, where input wave power is assumed to be 1W. We can see that the EM wave power is confined in the core of the CPOW. The transmittance to the output port |T | and the reflectivity |R| equal 0.986 and 0.007, respectively. Thus, the loss due to radiation from the CPOW is approximately 0.007. Figs. 12, 13 and 14 show the distributions excited by TE21 , TM01 and TM21 , respectively. Guided modes exists in each cases. The modes excited by TE21 and TM21

#22839 - $15.00 US

(C) 2000 OSA

Received June 01, 2000; Revised June 14, 2000

3 July 2000 / Vol. 7, No. 1 / OPTICS EXPRESS 14

are corresponding to the m = 0 cases, that is the hybrid modes. Transmittances are smaller than that of TE01 , because the input and output ports are optimized for the TE01 mode input. 5

Conclusion

We investigated the EM wave in a coaxial periodic optical waveguide. Existence of TEand TM-polarized guided modes and hybrid guided modes was demonstrated by using both Bessel functions expansion and the finite element method. The guided wave can propagate in the air hole with slight loss, and it can exist only in the stopband. By using the waveguide with the coaxial periodic structure, we can obtain fibers which have functions of filter, or which can reduce nonlinear effects.

Fig. 1. A waveguide structure with one-dimensional periodic structures.

#22839 - $15.00 US

(C) 2000 OSA

Received June 01, 2000; Revised June 14, 2000

3 July 2000 / Vol. 7, No. 1 / OPTICS EXPRESS 15

Fig. 2. Definition of propagation direction of a wave in a periodic structure.

Fig. 3. Coaxial periodic optical waveguide.

#22839 - $15.00 US

(C) 2000 OSA

Received June 01, 2000; Revised June 14, 2000

3 July 2000 / Vol. 7, No. 1 / OPTICS EXPRESS 16

Fig. 4. Cross-section.

#22839 - $15.00 US

(C) 2000 OSA

Received June 01, 2000; Revised June 14, 2000

3 July 2000 / Vol. 7, No. 1 / OPTICS EXPRESS 17

Fig. 5. Dispersion curves for k0 c = 1.0

Fig. 6. Dispersion curves for k0 c = 5.0

#22839 - $15.00 US

(C) 2000 OSA

Received June 01, 2000; Revised June 14, 2000

3 July 2000 / Vol. 7, No. 1 / OPTICS EXPRESS 18

Fig. 7. Stopband of a coaxial periodic optical waveguide.

Fig. 8. Stopband of a one-dimensonal photonic band gap structure.

#22839 - $15.00 US

(C) 2000 OSA

Received June 01, 2000; Revised June 14, 2000

3 July 2000 / Vol. 7, No. 1 / OPTICS EXPRESS 19

Fig. 9. z-direction components of magnetic field of TE-polarized mode and electric field of TM-polarized mode, where k = 1.2k0 and β=0.7859[TE], 0.5270[TM].

Fig. 10. z-direction components of magnetic field of TE-polarized mode and electric field of TM-polarized mode, where k = 3.0k0 and β=0.9672[TE], 0.9672[TM].

#22839 - $15.00 US

(C) 2000 OSA

Received June 01, 2000; Revised June 14, 2000

3 July 2000 / Vol. 7, No. 1 / OPTICS EXPRESS 20

Fig. 11. Intensity of EM wave field in CPOW, excited by TE01 mode. |T | = 0.986, |R| = 0.007.

Fig. 12. Intensity of EM wave field in CPOW, excited by TE21 mode. |T | = 0.479, |R| = 0.049.

#22839 - $15.00 US

(C) 2000 OSA

Received June 01, 2000; Revised June 14, 2000

3 July 2000 / Vol. 7, No. 1 / OPTICS EXPRESS 21

Fig. 13. Intensity of EM wave field in CPOW, excited by TM01 mode. |T | = 0.440, |R| = 0.106.

Fig. 14. Intensity of EM wave field in CPOW, excited by TM21 mode. |T | = 0.049, |R| = 0.493.

#22839 - $15.00 US

(C) 2000 OSA

Received June 01, 2000; Revised June 14, 2000

3 July 2000 / Vol. 7, No. 1 / OPTICS EXPRESS 22