Compact Supercell Method Based on Opposite Parity for Bragg Fibers

The supercellbased orthonormal basis method is proposed to investigate the modal properties of the Bragg fibers. A square lattice is constructed by the whole Bragg fiber which is considered as a supercell, and the periodical dielectric structure of the square lattice is decomposed using periodic functions (cosine). The modal electric field is expanded as the sum of the orthonormal set of Hermite-Gaussian basis functions based on the opposite parity of the transverse electric field. The propagation characteristics of Bragg fibers can be obtained after recasting the wave equation into an eigenvalue system. This method is implemented with very high efficiency and accuracy. ©2003 Optical Society of America OCIS codes: (060.2280) Fiber design and fabrication, (060.23 I 0) Fiber optics, (230. 7370) Waveguides References and links I. P.Yeh, A. Yariv, E. Marom, "Theory ofBragg fibers," J. Opt. Soc. Am., 68, 1196-1201(1978). 2. Y.Fink, J.N. Winn, S.H. Fan, C.Chen, J. Michel, J.D. Joannopoulos, E.L. Thomas, "A dielectric omnidirectional reflector," Science 282, 1679-1682(1998). 3. Y.Fink, D.J. Ripin, S.Fan, C. Chen, J.D. Joannopoulos, E.L. Thomas, "Guiding optical light in air using an all-dielcetric structure," J. Lightwave Technol. 17, 2039-2041 (1999). 4. M.lbanescu, Y.Fink, S.Fan, E.L. Thomas, J.D. Joannopoulos, "An all-dielectric coaxial waveguide," Science 289, 415-419(2000). 5. S.D. Hart, G.R.Maskaly, B.Temelkuran, P.H. Prideaux, J.D. Joannopoulos, Y.Fink, "External reflection from omnidirectional dielectric mirror fibers," Science 296, 510-513(2002). 6. P. Russell, "Photonic crystal fibers," Science 299, 358-362 (2003). 7. J.C. Knight, P. St. Russell, "New ways to guide light," Science 296, 276-277 (2002). 8. Jes Broeng, Stig E. B. Libori, T. Sondergaard, and Anders Bjarklev, "Analysis of air-guiding photonic bandgap fibers," Opt. Lett. 2S, 96-98 (2000). 9. M. Ibanescu, S.G. Johnson, M. Soljacic, J.D. Joannopoulos, Y.Fink, "Analysis of mode structure in hollow dielectric waveguide fibers," Phys. Rev. E 67, 046608-1-8(2003). 10. J.M. Bassett, A. Argyros, "Elimination of polarization degeneracy in round waveguides," Opt. Express 10, 1342-1346 (2002), http: //www.opticsexpress.org/ abstract. cfm? URI =OPEX-1023-1342. 11. S.G. Johnson, M. Ibanescu, M. Skorobogatiy, 0. Weisberg, T.D. Engeness, M. Soljacic, S.A. Jacobs, J.D. Joannopoulos, Y.Fink, "Low-loss asymptotically single-mode propagation in large-core omni guide fibers," Opt. Express 9, 748-779 (2001), http://www.opticsexpress.org/ abstract. cfm? URI =OPEX-913-748. 12. G. Ourang, Yong Xu, A. Yariv, "Theoretical study on dispersion compensation in air-core Bragg fibers," Opt. Express 10, 899-908 (2002), http: //www.opticsexpress.org/ abstract. cfm? URI =OPEX-1017-899. 13. T.D. Engeness, M. Ibanescu, S.G. Johnson, 0. Weisberg, M.Skorobogatiy, S.Jacobs, Y.Fink, "Dispersion tailoring and compensation by modal interactions in omni guide fibers," Opt. Express 11, 1175-1196 (2003), http: //www.opticsexpress.org/ abstract. cfm? URI =OPEX-11I0-1175. 14. J.A.Mosoriu, E. Silvestre, A. Ferrando, P.Andres, J.J . Miret, "High-index-core Bragg fibers: dispersion properties," Opt. Express 11 , 1400-1405 (2003), http://www.opticsexpress.org/ abstract. cfm? URI =OPEX11 12-1400. #3248 $15.00 us (C)2003 OSA Received October 30, 2003; Revised December I l, 2003 29 December 2003 / Vol. 11, No. 26 / OPTICS EXPRESS 3542 15. E. Silvestre, M.V. Andres, and P. Andres, "Biorthonormal-basis method for the vector description of optical-fiber mode," J. Lightwave Technol. 16, 923-928 (1998). 16. R. D. Meade, A.M. Rappe, K.D. Brommer, J.D. Joannopoulos, and 0.L. Alerhand, "Accurate theoretical analysis ofphotonic band-gap materials," Phys. Rev. B 48, 8434-8437 (1993). 17. Shangping Guo, Feng Wu, Khalid lkram, and Sacharia Albin, "Analysis of circular fibers with arbitrary index profile by Galerkin method," to be published by Opt. Lett .. 18. J.D. Joannopoulos, R.D. Meade, J.N. Winn, "Photonic crystals: molding the flow oflight," (New York, Princeton university press, 1995). 19. Shangping Guo, Sacharia Albin, "Simple plane wave implementation for photonic crystal calculations," Opt. Express 11, 167-175 (2003), http://www.opticsexpress.org/ abstract. cfm? URl=OPEX-l l-2-167. 20. A.W. Snyder, J.D. Love, Optical waveguide theory, (New York: Chapman and Hall, 1983). 21. T.M. Monro, DJ. Richardson, N.G.R. Broderick, P.J. Bennett, "Holey optical fibers: an efficient modal model," J. Ligbtwave Technol. 17, 1093-1102 (1999). 22. Yu-Li Hsueh, E. S. T. Hu, M. E. Marhic, . G. Kazovsky, "Opposite-parity orthonormal function expansion for efficient full-vectorial modeling of holey optical fibers," Opt. Lett. 28, 1188-1190(2003). 23. W. Zhi, R.G. Bin, L.S. Qin, and J.S.Sheng, "Supercell lattice method for photonic crystal fibers ," Opt. Express 11, 980-991 (2003), http:// www.opticsexpress.org/ abstract. cfm? UR!=OPEX-l 1-9-980. 24. R.G. Bin, W. Zhi, L.S. Qin, J.S. Sheng, "Mode Classification and Degeneracy in Photonic Crystal Fiber," Opt. Express 11, 1310-1321 (2003), http: //www.opticsexpress.org/abstract.cfm?URl=OPEX-1 l1 1I 3 l 0.


Introduction
At the beginning of the optical fiber progress, the ring fiber had been investigated by the transfer matrix method [ l].There was not much more attention paid to it because the index difference between every layer was very small.The ring fibers have attracted much interest because of their extraordinary properties when the great index difference dielectric materials were used at the end of last century (2)(3)(4)(5).The ring fiber has some new names called "Coaxial fiber", "Ornniguide fiber", or "Bragg fiber", in which the guiding of the light in the low-index core is due to the photonic band gap (PBG) produced by the periodic cladding, instead of the total internal reflection (TIR).
In the conventional metallic coaxial cable, the entire electromagnetic field is confined between two coaxial metal cylinders.The fundamental electromagnetic mode of a coaxial cable is the transverse electromagnetic (TEM) mode, which has radial symmetry in the electric field distribution and a linear relationship between frequency and wave vector.
Optical waveguide is restricted to the use of dielectric materials at optical wavelengths because of heavy absorption losses in the metal.However, because of the differences in boundary conditions of the electromagnetic fields at metal and dielectric surfaces, it has not previously been possible to recreate a TEM-like mode with all-dielectric materials.Consequently, optical waveguiding is done with the traditional total internal reflection (TIR) mechanism, by which the waveguides can achieve very low losses.
Recently, all-dielectric waveguides have been introduced that confine optical light by means of ID or 2D photonic band gap [6][7][8].These new designs have the potential advantage that light propagates mainly through the empty core of a hollow waveguide, thus minimizing effects associated with material nonlinearities, absorption losses and sharp bending.
The coaxial omniguide Bragg fiber, which combines some of the best features of the metallic coaxial cable and the dielectric waveguides, is an all-dielectric coaxial waveguide and supports a truly single mode in a low-index air core, which is very similar to the TEM mode of the metallic coaxial cable [2][3][4][5].It has a radially symmetric electric field distribution so that the polarization is maintained throughout propagation [4,9-11].It can be designed to be single-mode over a wide range of frequencies [11].In addition, the mode has a point of intrinsic zero dispersion around which a pulse can retain its shape during propagation [12][13][14].Finally, the coaxial omniguide can be used to guide light around sharp bends whose radius of curvature can be as small as the wavelength of the light [4,5].
There are a few different approaches to calculate the modes of the Bragg fibers, such as the semi-analytic approach based on the transfer matrix method [l], the asymptotic matrix method [11,12], the bi-orthonormal-basis method [14,15], and the approach involving a numerical solution of Maxwell's equations in the frequency domain with the use of the conjugate gradient method within the supercell approximation [16].The most recent Galerkin method [ 17] can be used to analyze the circular fibers with arbitrary index profile.
The supercell-based orthonormal basis method is proposed in this article.In section 2, the square lattice of the dielectric structure is constructed by arranging the whole Bragg fiber along x-and y-direction.In Section 3, the transverse electric field is decomposed into Hermite-Gaussian basis based on the opposite parity of the transverse electric field, and an eigenvalue system is deduced from the full-vectorial coupling wave equations.In Section 4, a Bragg fiber, which has the same structure as in Ref [ 11 ], is considered as an example, and the numerical results agree with it.Section 5 is the conclusion.

Supercell of the dielectric structure
The air-guiding Bragg fiber is a cylindrical multilayered fiber, which is a cylindrically symmetric microstructured fiber having a hollow core surrounded by a multilayered cladding made of alternating layers of a higher and a lower refractive-index dielectric.The multiayered structure geometry is characterized by the radial multiplayer period, A, and the low-indexlayer thichness, a, as illustrated in Fig. l(a).The dielectric constants are alternately e 1 and e 2 , where e 2 <e 1 • The dielectric constant of the hollow core is e 3 , and its radius is R .In order to analyze this fiber, a square lattice is constructed by the whole transverse profile of the Bragg fiber that is considered as a supercell, the lattice constant of the square lattice is D, as illustrated in Fig. l(b) .
., ---------- The periodic dielectric constant structure should be transferred into its Fourier transform when one is investigating photonic crystals, photonic crystal fibers or Bragg fibers [18,19], then the wave equation will be solved with the Fourier transform..e{r)={e,, where £; and r 1 are the dielectric constant and the outer radius of the ith layer, and r 0 is set zero, mis the total number of the layers, Eb is the background dielectric constant.Eq. (1) will be analytically expressed as [18]: where k=lkl,f;=m-/lA is the filling ratio.The limitation will be used when k=O in Eq. ( 3).
Because the dielectric constant structure of the Bragg fiber has the xand yaxial symmetric, i.e., £(-xJ1)=~x,-y)=£(xJ1), it can be expressed as a sum of the cosine functions as where (P+l) is the number of the expansion items, Pab, Pab1n are the expansion coefficients which can be analytically evaluated from the Fourier transform Ep(k) in Eq. ( 3) and expressed as Eq. ( 5).
,. Figure 2 is the simulation result of the dielectric constant structure of the Bragg fiber, which has the same structure parameters as the Fig. 4 in Ref. [11].The central area is the aircore with radius R=30A, and A=0.434µm.The dielectric constant of the periodical cladding are £ 1 =4.6 2 and £ 2 =1.6 2 , the background is Ez too.There are 17 layers, and the lower-index layer thichness is a=0.78A.The supercell lattice constant D is set as 1.2(2R+18A), the expansion parameter P=1200.

opposite parity of the mode field
When the confinement loss and the absorption loss are ignored in a normal waveguide which is uniform along the longitudinal direction (z-axis), both components e..(x.y) and e,,(x.y) of the transverse electric field must satisfy the full-vectorial coupling wave equations [20,21]: where pis the propagation constant corresponding to the mode field distribution (e,,, ey), ir-0=21TI). is the wave number of the vacuum.When f(x, y) is an even function about both the x and y directions, it can be proved that ex(x.y) and e,,(x.y)always have opposite parities in the x and y directions for each guided eigenmode (22].It is also to say that, when one transverse component is an even function about one axis, the other will be an odd function about the same axis.For example, if ex(x.y) is an even function about x, e,,(x.y)will be an odd function about x, i.e., e..(-x.y)=ex(x.y),e,,(x.y)=-e,,(x.y).

subscripts
For the case of compactness, two subscripts m and n are introduced to express the opposite parity of the mode electric field, which have the logical value O or 1, and are used to describe the symmetry of the x-component e..(x.y) as ex(-x.y)=(-lre..(x.y) and ex(x,-y)=(-l)"e..(x.y).All the compositions of 'mn' are [00, 01, 10, 11], which can completely express the symmetry of the mode electric field about both axes.For example, the subscript '10' means that ex(x.y) is an odd function about x and an even function about y, e,,(x.y) is an even function about x and an odd function about y.
In order to obtain the characteristics of the modes in the Bragg fibers, the transverse electric field can be expanded using the localized orthonormal Hermite-Gaussian basis functions as follows based on the opposite parity of the electric field.

.O a,b=O
where the bar over the subscript means the logical operator 'NOT', and the subscript 'mn' means that there are four sets of (e,,, ey) with different parity.In Eq. ( 7), Fis the number of the expansion terms, Eat,' (s=x,y) are the expansion coefficients, vr,{s) is the ith order orthonormal Hermite-Gaussian function which is defined as: where H,{s/m,) is the ith Hermite function which has the parity of H,{-slm,)=(-lYH,{slm,), m, is the character width [21,22].
the coupled terms f 4 ' JX and f 4 lY; the fiber with complex dielectric structure can be investigated while it is made of the dielectric constant with negative or positive imaginary part (gain or absorption); the degeneracy property can be obtained from the full vectorial coupling wave equation.The mode characteristics of the Bragg fiber are obtained by numerical computations in which the parameter Fis set as 15.The transverse electric field of the lowest-order mode, which is a degenerated pair of modes HE 11 x and HE 11 y, and the second order modes are quivered in Fig. 3 at wavelength 1550nm.It can be found that light is strongly confined in the low index air core by the PBG due to the outside annular Bragg reflector.The critical property of TEom (especially TEo 1 ) modes is that they have a node in their electric field (E,,) near r=R.The fractional IEl 2 in the cladding, hence the radiation loss, scales as 1/R 3 and 1/R for TEom and other modes (TM, HE, EH), respectively.Because of the substantial discrimination between the single lowest-loss mode TEo 1 and other higher-order guided modes, it allows even a highly multimode omniguide fiber to operate in an effectively single-mode fashion [9,11].

Numerical results
The propagation constant of mode TEo 1 is shown in Fig. 4(a), in which the solid line is obtained by the supercell method, and the cross is computed by the transfer matrix method (TMM).The difference between both approaches is demonstrated in Fig. 4(b ).In our computation for the case with A=0.434J.Un, J=l550nm, i.e. w=0.28x21te/A, which is same as Ref. [ 11 ], the propagation constant is 0.2792093x21t/ A, which has very good agreement with the value /J=0.27926x21t/A in Ref. [l l].It shows the reliability ofour model.

Conclusion
A compact supercell method based on the opposite parity of the transverse electric field for the Bragg fibers is proposed to investigate the modal properties in this paper.The modal electric field is expanded as the sum of the orthonormal set of Hermite-Gaussian basis functions, a square lattice is constructed by the whole transverse profile of the Bragg fiber which is considered as a supercell, and the periodical dielectric structure of the square lattice is decomposed using periodic functions (cosine).Considering the function relations and the orthonormality of Hermite-Gaussian functions, the propagation characteristics of Bragg fibers are obtained after recasting the full vectorial coupling wave equation into an eigenvalue system.Since the expansion of the transverse electric field is complete and compact due to the opposite parity in the even-parity dielectric waveguide, and all the decomposition coefficients and the overlap integrals are evaluated analytically, this model is efficient and accurate.This method is implemented with very high efficiency and agreement with some literatures [ 11,17].

Fig
Fig. I. scheme of the construction of the supercell square lattice of the Bragg fiber, (a) is the radial distribution of the dielectric constant and (b) is the supercell square lattice.