OSA's Digital Library

Optics Express

Optics Express

  • Editor: Michael Duncan
  • Vol. 13, Iss. 1 — Jan. 10, 2005
  • pp: 267–274
« Show journal navigation

Empirical relations for simple design of photonic crystal fibers

Kunimasa Saitoh and Masanori Koshiba  »View Author Affiliations


Optics Express, Vol. 13, Issue 1, pp. 267-274 (2005)
http://dx.doi.org/10.1364/OPEX.13.000267


View Full Text Article

Acrobat PDF (257 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

In order to simply design a photonic crystal fiber (PCF), we provide numerically based empirical relations for V parameter and W parameter of PCFs only dependent on the two structural parameters — the air hole diameter and the hole pitch. We demonstrate the accuracy of these expressions by comparing the proposed empirical relations with the results of full-vector finite element method. Through the empirical relations we can easily evaluate the fundamental properties of PCFs without the need for numerical computations.

© 2005 Optical Society of America

1. Introduction

Theoretical descriptions of PCFs have traditionally been based on numerical approaches, such as the plane wave expansion method [2

2. S.G. Johnson and J.D. Joannopoulos, “Block-iterative frequency-domain methods for Maxwell’s equations in a planewave basis,” Opt. Express8, 173–190 (2001), http://www.opticsexpress.org/abstract.cfm?URI=OPEX-8-3-173. [CrossRef] [PubMed]

], the multipole method (MM) [3

3. 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). [CrossRef]

], the finite element method (FEM) [4

4. M. Koshiba, “Full-vector analysis of photonic crystal fibers using the finite element method,” IEICE Trans. Electron. E85-C, 881–888 (2002).

, 5

5. K. Saitoh and M. Koshiba, “Full-vectorial imaginary-distance beam propagation method based on finite element scheme: Application to photonic crystal fibers,” IEEE J. Quantum Electron. 38, 927–933 (2002). [CrossRef]

], and so on, because of the relatively complex cross section of a PCF for which rotational symmetry is absent. However, numerical simulations are, in general, time-consuming and costly. Recently, an analytical approach based on the V parameter (normalized frequency) frequently used in the design of conventional optical fibers has been developed for index-guiding PCFs [6

6. M. Koshiba and K. Saitoh, “Applicability of classical optical fiber theories to holey fibers,” Opt. Lett. 29, 1739–1741 (2004). [CrossRef] [PubMed]

]. By appropriately defining the V parameter, various unique properties of PCFs can be qualitatively understood within the framework of well-established classical fiber theories without heavy numerical computations [6

6. M. Koshiba and K. Saitoh, “Applicability of classical optical fiber theories to holey fibers,” Opt. Lett. 29, 1739–1741 (2004). [CrossRef] [PubMed]

]. Although the V parameter offers a simple way to design a PCF, a limiting factor is that a numerical method is still required for obtaining the accurate effective cladding index. If we can get empirical relations for not only V parameter but also W parameter (normalized transverse attenuation constant) only dependent on the wavelength and the structural parameters, they would be very useful for simple design of PCFs.

The aim of this work is to provide the empirical relations for V parameter and W parameter of PCFs based on the fundamental geometrical parameters — the air hole diameter and the hole pitch. We demonstrate the accuracy of these expressions by comparing the proposed empirical relations with the results of full-vector FEM. Through the empirical relations we can easily evaluate the fundamental properties of PCFs without the need for numerical computations.

Fig. 1. Index-guiding photonic crystal fiber.

2. The V parameter expression

We consider a PCF with a triangular lattice of holes as shown in Fig. 1, where d is the hole diameter, Λ is the hole pitch, and the refractive index of silica is 1.45. In the center an air hole is omitted creating a central high index defect serving as the fiber core.

Recently, we have claimed that the triangular PCFs can be well parameterized in terms of the V parameter [6

6. M. Koshiba and K. Saitoh, “Applicability of classical optical fiber theories to holey fibers,” Opt. Lett. 29, 1739–1741 (2004). [CrossRef] [PubMed]

] that is given by

V=2πλaeffnco2nFSM2=U2+W2
(1)

with

U=2πλaeffnco2neff2
(2)
W=2πλaeffneff2nFSM2
(3)

where λ is the operating wavelength, nco is the core index, nFSM is the cladding index, defined as the effective index of the so-called fundamental space-filling mode in the triangular air hole lattice [7

7. T.A. Birks, J.C. Knight, and P.St.J. Russell, “Endlessy single-mode photonic crystal fiber,” Opt. Lett. 22, 961–963 (1997). [CrossRef] [PubMed]

], neff is the effective index of the fundamental guided mode, and aeff is the effective core radius that here is assumed to be Λ/3 [4

4. M. Koshiba, “Full-vector analysis of photonic crystal fibers using the finite element method,” IEICE Trans. Electron. E85-C, 881–888 (2002).

, 6

6. M. Koshiba and K. Saitoh, “Applicability of classical optical fiber theories to holey fibers,” Opt. Lett. 29, 1739–1741 (2004). [CrossRef] [PubMed]

]. The parameters U and W are called, respectively, the normalized transverse phase and attenuation constants. Mortensen et al. proposed the following effective V parameter [8

8. N.A. Mortensen, J.R. Folkenberg, M.D. Nielsen, and K.P. Hansen, “Modal cutoff and the V parameter in photonic crystal fibers,” Opt. Lett. 28, 1879–1881 (2003). [CrossRef] [PubMed]

] for triangular PCFs:

Veff=2πλΛneff2nFSM2
(4)

and reported the empirical relation for Veff of Eq. (4) [9

9. M.D. Nielsen and N.A. Mortensen, “Photonic crystal fiber design based on the V-parameter,” Opt. Express 11, 2762–2768 (2003), http://www.opticsexpress.org/abstract.cfm?URI=OPEX-11-21-2762. [CrossRef] [PubMed]

]. However, this definition is intrinsically different from the original V parameter definition in step-index fiber (SIF) theory and corresponds to the W parameter. Therefore it seems to be difficult to apply the design principle of SIFs straightforwardly to PCFs. So we adopt the V parameter definition of Eq. (1). Although we can estimate the fundamental properties of PCFs using the V parameter in Eq. (1) [6

6. M. Koshiba and K. Saitoh, “Applicability of classical optical fiber theories to holey fibers,” Opt. Lett. 29, 1739–1741 (2004). [CrossRef] [PubMed]

], a limiting factor for using Eq. (1) is that a numerical method is required for obtaining the accurate effective cladding index nFSM.

Figure 2 shows V values calculated through vector FEM [4

4. M. Koshiba, “Full-vector analysis of photonic crystal fibers using the finite element method,” IEICE Trans. Electron. E85-C, 881–888 (2002).

] as a function of λ/Λ for d/Λ ranging from 0.20 to 0.80 in steps of 0.05, where data are shown by open circles. By trial and error, we find that each data set in Fig. 2 can be fitted to a function of the form

V(λΛ,dΛ)=A1+A21+A3exp(A4λΛ)
(5)

and the results are indicated by the solid curves. For accurate fitting, the data sets are truncated at V=0.85. In Eq. (5) the fitting parameters Ai (i=1 to 4) depend on d/Λ only. The data are well described by the following expression

Ai=ai0+ai1(dΛ)bi1+ai2(dΛ)bi2+ai3(dΛ)bi3
(6)

and the coefficients a i0 to a i3 and b i1 to b i3 are given in Table 1. For λ/Λ<2 and V>0.85 the expression of Eq. (5) gives values of V which deviates less than 1.3% from the corrected values obtained from Eq. (1).

Using the effective V parameter in Eq. (5), the effective cladding index nFSM can be obtained without the need for numerical computations. Figure 3 shows nFSM as a function of λ/Λ for d/Λ ranging from 0.2 to 0.8 in steps of 0.1, where the open circles show numerical results from a fully vectorial FEM [4

4. M. Koshiba, “Full-vector analysis of photonic crystal fibers using the finite element method,” IEICE Trans. Electron. E85-C, 881–888 (2002).

] and the solid curves show the results from Eqs. (1) and (5) with aeff=Λ/√3. Mortensen et al. proposed the empirical expression for the value of nFSM [10

10. N.A. Mortensen, M.D. Nielsen, J.R. Folkenberg, A. Petersson, and H.R. Simonsen, “Improved large-mode-area endlessly single-mode photonic crystal fibers,” Opt. Lett. 28, 393–395 (2003). [CrossRef] [PubMed]

] to directly fit the effective cladding index, however, the results are not so accurate. On the other hand, the expression of Eq. (5) gives values of nFSM which deviates less than 0.25% from the values obtained through vector FEM for λ/Λ<1.5 and V>0.85.

As in Ref. [6

6. M. Koshiba and K. Saitoh, “Applicability of classical optical fiber theories to holey fibers,” Opt. Lett. 29, 1739–1741 (2004). [CrossRef] [PubMed]

], from Eq. (5) the cutoff condition is given by V=2.405, as in conventional SIFs. Using the empirical relation of Eq. (5) and various formulas in terms of the V parameter found for SIFs, we can easily estimate the fundamental properties of PCFs, such as mode field diameter, beam divergence, splice loss, and so on [6

6. M. Koshiba and K. Saitoh, “Applicability of classical optical fiber theories to holey fibers,” Opt. Lett. 29, 1739–1741 (2004). [CrossRef] [PubMed]

].

Table 1. Fitting coefficients in Eq. (6).

table-icon
View This Table
| View All Tables
Fig. 2. Effective V parameter as a function of λ/Λ.
Fig. 3. Effective cladding index as a function of λ/Λ.

3. The W parameter expression

Figure 4 shows W values calculated through vector FEM [5

5. K. Saitoh and M. Koshiba, “Full-vectorial imaginary-distance beam propagation method based on finite element scheme: Application to photonic crystal fibers,” IEEE J. Quantum Electron. 38, 927–933 (2002). [CrossRef]

] as a function of λ/Λ for d/Λ ranging from 0.20 to 0.80 in steps of 0.05, where data are shown by open circles. Again, by trial and error, we find that each data set in Fig. 4 can be fitted to the same function in Eq. (5) as

W(λΛ,dΛ)=B1+B21+B3exp(B4λΛ)
(7)

and the results are indicated by the solid curves. For accurate fitting, the data sets are truncated at W=0.1. In Eq. (7) the fitting parameters Bi (i=1 to 4) depend on d/Λ only. The data are well described by the following expression

Bi=ci0+ci1(dΛ)di1+ci2(dΛ)di2+ci3(dΛ)di3
(8)

and the coefficients c i0 to c i3 and d i1 to d i3 are given in Table 2. For λ/Λ<2 and W>0.1 the expression of Eq. (7) gives values of W which deviates less than 0.015 from the corrected values obtained from Eq. (3).

Using the V parameter in Eq. (5) and the W parameter in Eq. (7), the effective index of the fundamental mode neff can be obtained without the need for numerical computations. Figure 5 shows neff as a function of λ/Λ for d/Λ ranging from 0.2 to 0.8 in steps of 0.1, where the open circles show numerical results from a fully vectorial FEM [5

5. K. Saitoh and M. Koshiba, “Full-vectorial imaginary-distance beam propagation method based on finite element scheme: Application to photonic crystal fibers,” IEEE J. Quantum Electron. 38, 927–933 (2002). [CrossRef]

] and the solid curves show the results from Eqs. (1), (3), (5), and (7) with=Λ/√3 eff a. For λ/Λ<1.5 and W>0.1 the expressions of Eqs. (5) and (7) give values of neff which deviates less than 0.15% from the values obtained through vector FEM.

Table 2. Fitting coefficients in Eq. (8).

table-icon
View This Table
| View All Tables
Fig. 4. Effective W parameter as a function of λ/Λ.
Fig. 5. Effective index of the fundamental mode neff as a function of λ/Λ.

Next, using the empirical relations of Eqs. (5) and (7) we calculate the chromatic dispersion in PCFs. In order to use universal data for the effective index of the fundamental mode neff, we assume that the waveguide contribution to the dispersion parameter D is independent of material dispersion Dm, namely

D=λcd2neffdλ2+Dm
(9)

where c is the light velocity in a vacuum and Dm is given by the Sellmeier relation. Figure 6 shows the dispersion parameter D as a function of wavelength for d/Λ ranging from 0.2 to 0.8 in steps of 0.1, where the background index of silica is assumed to be 1.45, namely, nco=1.45. The results based on the empirical relations of Eqs. (5) and (7) agree well with the numerical results obtained by FEM [5

5. K. Saitoh and M. Koshiba, “Full-vectorial imaginary-distance beam propagation method based on finite element scheme: Application to photonic crystal fibers,” IEEE J. Quantum Electron. 38, 927–933 (2002). [CrossRef]

]. It is worth noting that, in Ref. [6

6. M. Koshiba and K. Saitoh, “Applicability of classical optical fiber theories to holey fibers,” Opt. Lett. 29, 1739–1741 (2004). [CrossRef] [PubMed]

], the chromatic dispersion of PCFs was calculated by using the Gloge formula [11

11. D. Gloge, “Weakly guiding fibers,” Appl. Opt. 10, 2252–2258 (1971). [CrossRef] [PubMed]

] and V values, while, here, the expressions of Eqs. (5) and (7) are used for direct calculation of the chromatic dispersion.

Fig. 6. Chromatic dispersion as a function of wavelength for (a) Λ=2.0 µm, (b) Λ=2.5 µm, and (c) Λ=3.0 µm. Solid curves, results of empirical relations; dashed curves, results of vector FEM.

4. Conclusions

In order to simply design a PCF, we provided the empirical relations for both V parameter and W parameter of PCFs only dependent on the air hole diameter and the hole pitch. We demonstrated the accuracy of these expressions by comparing the proposed empirical relations with the results of full-vector FEM. Through the empirical relations the fundamental properties of PCFs could be easily estimated without the need for numerical computations.

References and links

1.

P.St.J. Russell, “Photonic crystal fibers,” Science 299, 358–362 (2003). [CrossRef] [PubMed]

2.

S.G. Johnson and J.D. Joannopoulos, “Block-iterative frequency-domain methods for Maxwell’s equations in a planewave basis,” Opt. Express8, 173–190 (2001), http://www.opticsexpress.org/abstract.cfm?URI=OPEX-8-3-173. [CrossRef] [PubMed]

3.

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). [CrossRef]

4.

M. Koshiba, “Full-vector analysis of photonic crystal fibers using the finite element method,” IEICE Trans. Electron. E85-C, 881–888 (2002).

5.

K. Saitoh and M. Koshiba, “Full-vectorial imaginary-distance beam propagation method based on finite element scheme: Application to photonic crystal fibers,” IEEE J. Quantum Electron. 38, 927–933 (2002). [CrossRef]

6.

M. Koshiba and K. Saitoh, “Applicability of classical optical fiber theories to holey fibers,” Opt. Lett. 29, 1739–1741 (2004). [CrossRef] [PubMed]

7.

T.A. Birks, J.C. Knight, and P.St.J. Russell, “Endlessy single-mode photonic crystal fiber,” Opt. Lett. 22, 961–963 (1997). [CrossRef] [PubMed]

8.

N.A. Mortensen, J.R. Folkenberg, M.D. Nielsen, and K.P. Hansen, “Modal cutoff and the V parameter in photonic crystal fibers,” Opt. Lett. 28, 1879–1881 (2003). [CrossRef] [PubMed]

9.

M.D. Nielsen and N.A. Mortensen, “Photonic crystal fiber design based on the V-parameter,” Opt. Express 11, 2762–2768 (2003), http://www.opticsexpress.org/abstract.cfm?URI=OPEX-11-21-2762. [CrossRef] [PubMed]

10.

N.A. Mortensen, M.D. Nielsen, J.R. Folkenberg, A. Petersson, and H.R. Simonsen, “Improved large-mode-area endlessly single-mode photonic crystal fibers,” Opt. Lett. 28, 393–395 (2003). [CrossRef] [PubMed]

11.

D. Gloge, “Weakly guiding fibers,” Appl. Opt. 10, 2252–2258 (1971). [CrossRef] [PubMed]

OCIS Codes
(060.2280) Fiber optics and optical communications : Fiber design and fabrication
(060.2400) Fiber optics and optical communications : Fiber properties

ToC Category:
Research Papers

History
Original Manuscript: December 8, 2004
Revised Manuscript: December 23, 2004
Manuscript Accepted: December 30, 2004
Published: January 10, 2005

Citation
Kunimasa Saitoh and Masanori Koshiba, "Empirical relations for simple design of photonic crystal fibers," Opt. Express 13, 267-274 (2005)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-13-1-267


Sort:  Author  |  Year  |  Journal  |  Reset  

References

  1. P.St.J.  Russell, “Photonic crystal fibers,” Science 299, 358–362 (2003). [CrossRef] [PubMed]
  2. S.G.  Johnson, J.D.  Joannopoulos, “Block-iterative frequency-domain methods for Maxwell’s equations in a planewave basis,” Opt. Express 8, 173–190 (2001), http://www.opticsexpress.org/abstract.cfm?URI=OPEX-8-3-173. [CrossRef] [PubMed]
  3. T.P.  White, B.T.  Kuhlmey, R.C.  McPhedran, D.  Maystre, G.  Renversez, C.M.  de Sterke, L.C.  Botten, “Multipole method for microstructured optical fibers. I. Formulation,” J. Opt. Soc. Am. B 19, 2322–2330 (2002). [CrossRef]
  4. M.  Koshiba, “Full-vector analysis of photonic crystal fibers using the finite element method,” IEICE Trans. Electron. E85-C, 881–888 (2002).
  5. K.  Saitoh, M.  Koshiba, “Full-vectorial imaginary-distance beam propagation method based on finite element scheme: Application to photonic crystal fibers,” IEEE J. Quantum Electron. 38, 927–933 (2002). [CrossRef]
  6. M.  Koshiba, K.  Saitoh, “Applicability of classical optical fiber theories to holey fibers,” Opt. Lett. 29, 1739–1741 (2004). [CrossRef] [PubMed]
  7. T.A.  Birks, J.C.  Knight, P.St.J.  Russell, “Endlessy single-mode photonic crystal fiber,” Opt. Lett. 22, 961–963 (1997). [CrossRef] [PubMed]
  8. N.A.  Mortensen, J.R.  Folkenberg, M.D.  Nielsen, K.P.  Hansen, “Modal cutoff and the V parameter in photonic crystal fibers,” Opt. Lett. 28, 1879–1881 (2003). [CrossRef] [PubMed]
  9. M.D.  Nielsen, N.A.  Mortensen, “Photonic crystal fiber design based on the V-parameter,” Opt. Express 11, 2762–2768 (2003), http://www.opticsexpress.org/abstract.cfm?URI=OPEX-11-21-2762. [CrossRef] [PubMed]
  10. N.A.  Mortensen, M.D.  Nielsen, J.R.  Folkenberg, A.  Petersson, H.R.  Simonsen, “Improved large-mode-area endlessly single-mode photonic crystal fibers,” Opt. Lett. 28, 393–395 (2003). [CrossRef] [PubMed]
  11. D.  Gloge, “Weakly guiding fibers,” Appl. Opt. 10, 2252–2258 (1971). [CrossRef] [PubMed]

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