Andrew M. Weiner, Editor-in-Chief
Chin-Ping Yu and Hung-Chun Chang
Chin-Ping Yu and Hung-Chun Chang†
Graduate Institute of Electro-Optical Engineering National Taiwan University, Taipei, Taiwan 106-17, R.O.C.
†also with the Department of Electrical Engineering and the Graduate Institute of Communication Engineering, National Taiwan University
Eigenvalue equations for solving full-vector modes of optical waveguides are formulated using Yee-mesh-based finite difference algorithms and incorporated with perfectly matched layer absorbing boundary conditions. The established method is thus able to calculate the complex propagation constants and the confinement losses of leaky waveguide modes. Proper matching of dielectric interface conditions through the Taylor series expansion of the fields is adopted in the formulation to achieve high numerical accuracy. The method is applied to the study of the holey fiber with triangular lattice, the two-core holey fiber, and the air-guiding photonic crystal fiber.
© 2004 Optical Society of America
Po-jui Chiang and Hung-chun Chang
Opt. Express 19(2) 1594-1608 (2011)
Ming-yun Chen, Sen-ming Hsu, and Hung-chun Chang
Opt. Express 17(8) 5965-5979 (2009)
Zhaoming Zhu and Thomas G. Brown
Opt. Express 10(17) 853-864 (2002)
Jinbiao Xiao and Xiaohan Sun
Opt. Express 20(19) 21583-21597 (2012)
Chih-Hsien Lai and Hung-chun Chang
Opt. Express 19(2) 562-569 (2011)
S. Guo, F. Wu, and S. Albin, “Loss and dispersion analysis of microstructured fibers by finite-difference method,” Opt. Express 12, 3341–3352 (2004), http://www.opticsexpress.org/abstract.cfm?URI=OPEX-12-15-3341
K. Saitoh, N. A. Mortensen, and M. Koshiba, “Air-core photonic band-gap fibers: the impact of surface modes,” Opt. Express 12, 394–400 (2004), http://www.opticsexpress.org/abstract.cfm?URI=OPEX-12-3-394
H. K. Kim, J. Shin, S. Fan, M. J. F. Digonnet, and G. S. Kino, “Designing air-core photonic-bandgap fibers free of surface modes,” IEEE J. Quantum Electron. 40, 551–556 (2004).
Y. C. Chiang, Y. P. Chiou, and H. C. Chang, “Improved full-vectorial finite-difference mode solver for optical waveguides with step-index profiles,” J. Lightwave Technol. 20, 1609–1618 (2002).
G. R. Hadley, “High-accuracy finite-difference equations for dielectric waveguide analysis I: Uniform regions and dielectric interfaces,” J. Lightwave Technol. 20, 1210–1218 (2002).
G. R. Hadley, “High-accuracy finite-difference equations for dielectric waveguide analysis II: Dielectric corners,” J. Lightwave Technol. 20, 1219–1231 (2002).
K. Saitoh and M. Koshiba, “Full-vectorial imaginary-distance beam propagation method based on a finite element scheme: Application to photonic crystal fibers,” IEEE J. Quantum Electron. 38, 927–933 (2002).
T. P. White, B. T. Kuhlmey, R. C. Mcphedran, D. Maystre, G. Renversez, C. M. de Sterke, and L. C. Botten, “Multipole method for microstructured optical fibers. I. Formulation,” J. Opt. Soc. Am. B 19, 2322–2330 (2002).
T. Ando, H. Nakayama, S. Numata, J. Yamauchi, and H. Nakano, “Eigenmode analysis of optical waveguides by a Yee-mesh-based imaginary-distance propagation method for an arbitrary dielectric interface,” J. Lightwave Technol. 20, 1627–1634 (2002).
Z. Zhu and T. G. Brown, “Full-vectorial finite-difference analysis of microstructured optical fibers,” Opt. Express 10, 853–864 (2002), http://www.opticsexpress.org/abstract.cfm?URI=OPEX-10-17-853
M. Koshiba and K. Saitoh, “Numerical verification of degeneracy in hexagonal photonic crystal fibers,” IEEE Photon. Technol. Lett. 13, 1313–1315 (2001).
Y. P. Chiou, Y. C. Chiang, and H. C. Chang, “Improved three-point formulas considering the interface conditions in the finite-difference analysis of step-index optical devices,” J. Lightwave Technol. 18, 243–251 (2000).
R. F. Cregan, B. J. Mangan, J. C. Knight, T. A. Birks, P. St. J. Russell, P. J. Roberts, and D. C. Allan, “Single-mode photonic band gap guidance of light in air,” Science 285, 1537–1539 (1999).
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).
R. Mittra and U. Pekel, “A new look at the perfectly matched layer (PML) concept for the reflectionless absorption of electromagnetic waves,” IEEE Microwave Guided Wave Lett. 5, 84–86 (1995).
C. M. Rappaport, “Perfectly matched absorbing boundary conditions based on anisotropic lossy mapping of space,” IEEE Microwave Guided Wave Lett. 5, 90–92 (1995).
C. L. Xu, W. P. Huang, M. S. Stern, and S. K. Chaudhuri, “Full-vectorial mode calculations by finite difference method,” Inst. Elec. Eng. Proc.-J 141, 281–286 (1994).
P. Lüsse, P. Stuwe, J. Schüle, and H.-G. Unger, “Analysis of vectorial mode fields in optical waveguides by a new finite difference method,” J. Lightwave Technol. 12, 487–494 (1994).
J. P. Berenger, “A perfectly matched layer for the absorption of electromagnetic waves,” J. Computat. Phys. 114, 185–200 (1994).
W. C. Chew and W. H. Weedon, “A 3D perfectly matched medium from modified Maxwell’s equations with stretched coordinates,” Microwave Opt. Technol. Lett. 7, 599–604 (1994).
K. S. Yee, “Numerical solution of initial boundary value problems involving Maxwell’s equations in isotropic media,” IEEE Trans. Antennas Propagat. AP-14, 302–307 (1966).
H. C. Chang and C. P. Yu, “Yee-mesh-based finite difference eigenmode analysis algorithms for optical waveguides and photonic crystals,” in OSA Integr. Photon. Res. Dig., 2004, paper IFE4.
H. C. Chang and C. P. Yu, “Yee-mesh-based finite difference optical waveguide eigenmode solver with perfectly matched layer absorbing boundaries conditions,” in Proceedings of 5th Asia-Pacific Engineering Research Forum on Microwaves and Electromagnetic Theory, pp. 158–164, Kyushu University, Fukuoka, Japan, July 29–30, 2004.
Ditkowski, J. S. Hesthaven, and C. H. Teng, “Modeling dielectric interfaces in the FDTD-method: A comparative study,” in 2000 Progress in Electromagnetics Research (PIERS 2000) Proceedings, Cambridge, Massachusetts, July 5–14, 2000.
OSA participates in Crossref's Cited-By Linking service. Citing articles from OSA journals and other participating publishers are listed here.
Alert me when this article is cited.
Click here to see a list of articles that cite this paper
The cross-section of an arbitrary waveguide problem with the PMLs placed at the edges of the computing domain.
Download Full Size | PPT Slide | PDF
Yee’s 2-D mesh for the FDFD method.
The situation of a dielectric interface lying between two sampled points.
Relative errors in the modal index of the fundamental TE-like mode for the square channel waveguide using the present method with three different schemes in dealing with the dielectric interface.
Relative errors in the modal index of the fundamental mode and the corresponding computation time for the strongly guiding optical fiber using the present method with the proper BC matching scheme and the stair-case approximation.
The computing window with PMLs for the holey fiber with three-ring air holes.
(a) The effective indices and (b) the losses of the x-polarized fundamental guided modes in the three-ring holey fiber with a=2.3 µm.
Modal birefringence of the one-ring HF using the present method with the proper BC matching scheme and the stair-case approximation.
The computing window with PMLs for the two-core holey fiber.
(a) The effective indices and (b) the losses of the even modes in the two-core holey fiber with r/a=0.25, 0.30, and 0.35.
(a) The effective indices and (b) the losses of the odd modes in the two-core holey fiber with r/a=0.25, 0.30, and 0.35.
The computing window with PMLs for a hollow PCF with six rings of air holes in the cladding.
(a) The effective indices of the surface modes and fundamental core modes of the hollow PCF of Fig. 11 with air filling fraction f=0.7. (b) Losses of the core modes and the surface modes.
The field distributions of the surface modes and fundamental core modes at points A, B, C, and D indicated in Fig. 13(a).
Table 1. The values of sx
for the PML region I, II, and III.
Table 2. The calculated modal index as a function of the number of grid points along the x-axis and the reference result is n
Equations on this page are rendered with MathJax. Learn more.
The values of sx
for the PML region I, II, and III.
The calculated modal index as a function of the number of grid points along the x-axis and the reference result is n