OSA's Digital Library

Optics Express

Optics Express

  • Editor: C. Martijin de Sterke
  • Vol. 15, Iss. 9 — Apr. 30, 2007
  • pp: 5342–5359
« Show journal navigation

Determination of the mode reflection coefficient in air-core photonic bandgap fibers

Vinayak Dangui, Michel J. F. Digonnet, and Gordon S. Kino  »View Author Affiliations


Optics Express, Vol. 15, Issue 9, pp. 5342-5359 (2007)
http://dx.doi.org/10.1364/OE.15.005342


View Full Text Article

Acrobat PDF (1193 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

Using an eigenmode decomposition technique, we numerically determine the backreflection coefficient of the modes of air-core photonic bandgap fibers for flat terminations. This coefficient is found to be very small for the fundamental air-guided mode, of the order of 10-5 to 10-6, in contrast with the surface and bulk modes, which exhibit significantly higher reflections, by about three to four orders of magnitude. For the Crystal Fibre HC-1550-2 fiber, we find a reflection coefficient of 1.9×10-6 for an air termination, and approximately 3.3% for a silica termination. We also find that the Fresnel approximation is ill suited for the determination of the modal reflection coefficient, and instead propose a more accurate new formula based on an averaged modal index.

© 2007 Optical Society of America

1. Introduction

Air-core photonic bandgap fibers (PBFs) [1–3

1. J. Broeng, D. Mogilevstev, S. E. Barkou, and A. Bjarklev, “Photonic crystal fibers: A new class of optical waveguides,” Opt. Fiber Technol. 5, 305–330 (1999). [CrossRef]

] offer the enticing theoretical prospect of very low nonlinearities, reduced sensitivity to environmental factors such as temperature [4

4. V. Dangui, H. K. Kim, M. J. F. Digonnet, and G. S. Kino, “Phase sensitivity to temperature of the fundamental mode in air-guiding photonic-bandgap fibers,” Opt. Express 13, 6669–6684 (2005). [CrossRef] [PubMed]

], potentially low backscattering and hence very low propagation losses, and they are finding interesting novel applications in various photonics fields, including fiber sensors and high-power laser delivery. As is the case for conventional solid-core fibers, the performance of a fiber optic system based on PBFs is expected to be adversely affected by the reflection that takes place at the end of an air-core fiber, whether the fiber is terminated in air (free end) or butt-coupled or spliced to another air-core or solid-core fiber. Backreflection from such terminations are notoriously detrimental to the performance of fiber interferometers, communication systems, and lasers and amplifiers, as they result in increased noise, instabilities, and other deleterious effects. [5

5. D. M. Dagenais, K. P. Koo, and F. Bucholtz, “Fiber interferometry limitations due to parasitic optical cavities,” Lasers and Electro-Optics Society Conference, San Jose, CA USA (1993).

] While the reflection at the end of a photonic crystal fiber has been studied, [6

6. J. Corbett, A. Dabirian, T. Butterley, N. A. Mortensen, and J. R. Allington-Smith, “The coupling performance of photonic crystal fibres in fibre stellar interferometry,” Mon. Not. Roy. Astron. Soc. 368, 203–210 (2006).

] evaluating the magnitude of the reflection at the end of an air-core fiber, such as the HC-1550-2, [7

7. Crystal Fibre website, http://www.crystal-fibre.com

] has never been reported before and it is therefore becoming a pressing issue.

The backreflection from the end of a conventional fiber terminated by a semi-infinite medium of refractive index n 0 is characterized by a power reflection coefficient that is usually calculated by the well-known approximate Fresnel formula:

R=(neffn0neff+n0)2,
(1)

where neff is the effective index of the mode considered. The Fresnel approximation is equivalent to replacing the waveguide with a medium of uniform refractive index equal to the mode effective index, and to approximate the mode by a plane wave at normal incidence. As a result, the Fresnel approximation performs best for weakly-guiding waveguides with a small index contrast across the waveguide cross-section, for modes possessing a small numerical aperture (lower order modes), and for a strong refractive index contrast between the waveguide and the termination medium. Conversely, this approximation has proved to be inaccurate to evaluate the reflection coefficient of strongly guiding waveguides, [8

8. D. Khalil, “Reflection at the end of strongly guiding dielectric waveguide,” Opto-Electronics and Communications Conference, Yokohama, Japan , 1, 352–353 (2002).

] and it does not account for the spatial intensity distribution of the mode. Because PBFs are strongly guided waveguides, and because we are mostly concerned with the reflection of a leaky mode (the fundamental air-guided core mode) into air (low refractive index contrast between waveguide and termination medium), we expect the Fresnel approximation to be ineffective in this context, and that a more elaborate calculation will be required.

In this paper, we detail for the first time a computation procedure for the derivation of the reflection coefficient of a PBF mode into a medium of arbitrary refractive index, which we then apply to calculate numerically the reflection coefficient of the fundamental core mode, as well as representative surface and bulk modes, of an air-core fiber, and their dependencies on wavelength and fiber parameters. We then study the influence of PBF design parameters on the reflection coefficient, in order to design PBF components with very low reflections for improved performance in fiber devices. We also derive a simple formula based on the Fresnel approximation that provides accurate predictions of the backreflection from the end of a PBF, especially for fibers with a high air-filling ratio (above 90 %).

2. Theory

The configuration studied is represented in Fig. 1. It consists of a semi-infinite PBF, extending over the z < 0 region (region I). The z > 0 region (region II) is a semi-infinite medium of uniform index n 0, such as air. The PBF parameters are as follows: the cladding crystal periodicity is Λ, the core radius is R, and the cladding air hole radius is ρ. The transverse plane is hereafter defined as the interface plane terminating the PBF. We consider the reflection at this interface of a mode field propagating through the PBF in the +z direction.

Fig. 1. Configuration studied for the computation of the modal reflection coefficient of a PBF terminated in a uniform medium. The boundary between regions I and II is marked by the red dotted line, and defines the transverse plane.

To calculate the reflection coefficient, we use the formalism of eigenmode decomposition.[9

9. P. Gerard, P. Benech, D. Khalil, R. Rimet, and S. Tedjini, “Towards a full vectorial and modal technique for the analysis of integrated optics structures: The Radiation Spectrum Method (RSM),” Opt. Commun. 140, 128–145 (1997). [CrossRef]

,10

10. D. Marcuse, Theory of dielectric optical waveguides, (Academic Press, 1974).

] The incident and reflected fields (in region I) are decomposed into a sum of PBF modes and the transmitted field (in region II) into a sum of plane waves. The transverse incident electrical field is denoted EI(r) and can either be a core mode, a cladding mode, or a surface mode of the PBF, labeled E0(r). The reflected transverse electrical field, labeled ER(r), propagates in the -z direction and is written as a superposition of core and cladding mode fields of the PBF, which are collectively labeled Em(r).

ER(r)=m=0NMamEm(r),
(2)

where the summation is carried over (NM+1) modes, and am is the amplitude of the signal reflected into mode m. Furthermore, only modes of the same symmetry class as the incident mode will be present in the decomposition of the reflected field. For most application, we are concerned with the amount of light backreflected into the incident mode, which is labeled as m = 0, since all other coefficients are much smaller (by many orders of magnitude). The intensity or power reflection coefficient into the fundamental mode is:

R=a02,
(3)

In the same fashion, the transverse transmitted electrical field propagating in region II in the positive z direction is denoted ET(r) and decomposed in a sum of plane waves propagating in free space. The transverse magnetic fields are decomposed in the same fashion. Using the continuity of the transverse electric and magnetic fields, we obtain the following relationships at the PBF interface plane:

EI+ER=ET,
(4)
HI+HR=HT,
(5)

As described in the appendix, we can then use the orthogonality of the transverse fields of the waveguide modes, and the orthogonality of the transverse fields of the plane waves, to transform these equations into a system of linear equations, which are then solved numerically to recover the reflection coefficients am.

The PBFs were modeled using the Stanford Photonic Bandgap Fiber (SPBF) code, a finite-difference mode solver developed for the purpose of modeling the modes of waveguides with complex refractive index profiles such as air-core fibers. [11

11. V. Dangui, M. J. F. Digonnet, and G. S. Kino, “A fast and accurate numerical tool to model the modal properties of photonic-bandgap fibers,” Opt. Express 14, 2979–2993 (2006). [CrossRef] [PubMed]

] In the simulations carried out for this paper, the SPBF code was used with periodic boundary conditions, with a simulation step of Λ/50 and a grid size of 10Λ×10Λ. We used typically 260 PBF modes and between 2000 and 5000 plane waves for the computation of the reflection coefficients.

Fig. 2. Typical matrix A used for computation of the modal backreflection.

To compute the projection of the PBF modes onto plane waves, we used a two-dimensional discrete Fourier transform (2D DFT) of the computed fields. Figure 2 illustrates a typical matrix A. The non-zero coefficients are illustrated as blue dots, and the zero coefficients are left out (shown as white spaces). The upper left corner indicates the projection of PBF modes onto each other; the lower right corner shows the projection of plane waves onto each other. The remaining portions of the matrix show the projection of PBF modes onto plane waves propagating in air. The matrix size is 2082 × 2082 elements, of which nearly 950,000 are non-zero complex elements. This sheer number of non-zero elements makes it impossible to invert such a matrix accurately directly, so we resorted to an iterative method, described in Ref [12

12. H. . van der Vorst, “BI-CGSTAB: A fast and smoothly converging variant of BI-CG for the solution of nonsymmetric linear systems,” SIAM J. Sci. Stat. Comput. 13, 631–644 (1992). [CrossRef]

].

3. Variation of the fundamental mode reflection coefficient across the PBF bandgap

Fig. 3. Refractive index profile of the PBF modeled in this section, for which ρ=0.47Λ and R=0.8Λ. The gray areas indicate air (refractive index equal to 1); the black areas indicate silica (refractive index chosen equal to 1.45).

In the first part of this work, we modeled the PBF whose refractive index profile is illustrated in Fig. 3. The cladding holes' radius is ρ = 0.47Λ, and the core radius R = 0.8Λ. Figure 4 represents the calculated bandgap of this fiber, as well as the dispersion curves of all of its core modes. The bandgap extends from λ = 0.56Λ to λ = 0.64Λ. Due to the choice of core radius, [13

13. H. K. Kim, J. Shin, S. H. Fan, M. J. F. Digonnet, and G. S. Kino, “Designing air-core photonic-bandgap fibers free of surface modes,” IEEE J. Quantum Eectron. 40, 551–556 (2004). [CrossRef]

] this fiber does not support surface modes, and it has only two core modes, the two degenerate, orthogonally polarized fundamental modes. Figure 5 shows the longitudinal Poynting vector profile of one of the fundamental modes at λ = 0.60Λ. The fundamental core modes have a very large fraction of their intensity (over 90 %) concentrated in air, which is the reason why a weak reflection at an air interface is expected.

The Fourier transform intensity of this same mode's transverse electric field, plotted in logarithmic scale in the wavevector domain, is shown in Fig. 6. The black circle shows the boundary for all wavevectors such that kx2+ky2k02=(2πλ)2, i.e., all wavevectors admissible for plane waves propagating in air. Since the fundamental mode profile exhibits little variation across the bandgap, as the wavelength varies across the bandgap, the most significant change occurs in the value of k 0, and thus in the radius of this circle. As the wavelength increases, the radius of the circle decreases, and thus fewer plane waves can be excited past the PBF termination plane, yielding a decrease in transmitted field. We therefore expect the fundamental-mode reflection to increase with wavelength. At the same time, since the circle boundary lies in a region where the Fourier transform intensity of the electric field is very weak (roughly 0.1% of its maximum value), we expect a weak dependence on wavelength.

Fig. 4. Bandgap and dispersion curve of the PBF shown in Fig. 3. The fundamental modes are depicted in red, and the bandgap limits are shown in black.
Fig. 5. Fundamental mode longitudinal Poynting vector profile, at a wavelength of λ = 0.6Λ, for the PBF of Fig. 3.
Fig. 6. Fourier transform intensity, in logarithmic scale, of the transverse electric field of the mode plotted in Fig. 5. The black circle marks all wavevectors such that k 2 x+k 2 yk 2 0, corresponding to plane waves that can be excited in air.

This is indeed found to be the case. Figure 7 shows the reflection coefficient of the x-polarized fundamental mode, defined as |a 0|2, when the fiber is terminated in air, calculated at several wavelengths across the bandgap. The reflection coefficient remains weak across the entire bandgap, its coefficient varying between ~0.6×10-5 at the short-wavelength end of the bandgap and ~1.3×10-5 at the long-wavelength end. We also observe a general trend towards an increase in the backreflection as the wavelength increases, which was predicted earlier on the basis of the smaller wavevectors associated with longer wavelengths: the mode energy is transmitted in air over a smaller range of wavevectors, and fewer plane waves can be excited, thus resulting in an increased reflection coefficient.

Fig. 7. Fundamental mode reflection coefficient for an air termination, for the PBF of Fig. 3, across the bandgap, as a function of the normalized wavelength λ/Λ. The incident field is the mostly x-polarized fundamental mode.

4. Reflection coefficient of various PBF modes

It was also interesting to determine the reflection coefficient for the other kinds of modes supported by air-core fibers, in particular higher order core modes, surface modes, and cladding modes. To do so, we retained the same cladding holes radius (ρ = 0.47Λ) but increased the fiber core radius to R = 2.0Λ, so that the fiber now supports a large number of core modes, as well as some surface modes. As is well known, the surface-mode energy is concentrated mostly in the silica at the core-cladding boundary. The core region of this fiber is shown in Fig. 8, together with the intensity contour lines of an exemplary surface mode.

Fig. 8. Surface mode electric field intensity profile for a PBF with R = 2Λ and ρ = 0.47Λ at a wavelength λ = 0.6Λ.

The reflection coefficient was also calculated for multiple representative bulk modes (also for reflection into the backward-traveling fundamental mode of same polarization). All reflection coefficients were approximately equal to 1%. This is again entirely consistent with the field distributions of these bulk modes, which are strongly localized on the thin silica veins of the photonic-crystal cladding. This feature brings their reflection coefficient up to a value of the same magnitude as a reflection from silica into air (~3.37%).

Fig. 9. Comparison between the computed reflection coefficients and the Fresnel approximation predictions for a few surface and core modes at a wavelength of λ = 0.6Λ, for a PBF with ρ = 0.47Λ and R = 2Λ.

5. Effect of PBF design parameters on the fundamental core mode reflection coefficient

Unlike in a conventional fiber, the magnitude of the backreflection at the end of a PBF is strongly influenced by the fiber parameters. These parameters can therefore be adjusted to minimize backreflection. One such parameter is the core radius. A proper choice of core radius can prevent the PBF from supporting surface modes. [13

13. H. K. Kim, J. Shin, S. H. Fan, M. J. F. Digonnet, and G. S. Kino, “Designing air-core photonic-bandgap fibers free of surface modes,” IEEE J. Quantum Eectron. 40, 551–556 (2004). [CrossRef]

] At the same time, the absence of surface modes has been linked to a reduction in the fraction of fundamental mode power localized in the silica regions of the PBF. [4

4. V. Dangui, H. K. Kim, M. J. F. Digonnet, and G. S. Kino, “Phase sensitivity to temperature of the fundamental mode in air-guiding photonic-bandgap fibers,” Opt. Express 13, 6669–6684 (2005). [CrossRef] [PubMed]

] We thus expect that the reflection of the fundamental mode will be higher for core radii that support surface modes. Since larger cores support fundamental modes that are more strongly confined in the air core (i.e., the fundamental mode intensity at the core-cladding interface is weaker), we also expect that as the core radius increases, the fundamental mode reflection should generally decrease.

To confirm these trends, we plot in Fig. 10 the calculated dependence of the fundamental mode reflection coefficient on the core radius, for the same cladding holes radius as used before (ρ = 0.47Λ) and at about the center of the bandgap (wavelength λ = 0.6Λ). As expected, as the core radius increases, the reflection tends to decrease, except for radii where the PBF supports surface modes, such as for radii between 1.05Λ and 1.3Λ, and between 1.45Λ and 1.7Λ.[13

13. H. K. Kim, J. Shin, S. H. Fan, M. J. F. Digonnet, and G. S. Kino, “Designing air-core photonic-bandgap fibers free of surface modes,” IEEE J. Quantum Eectron. 40, 551–556 (2004). [CrossRef]

] For those radii, the fundamental mode reflection coefficient increases very sharply and significantly, again because of the substantial increase in the fraction of fundamental mode energy in silica.[4

4. V. Dangui, H. K. Kim, M. J. F. Digonnet, and G. S. Kino, “Phase sensitivity to temperature of the fundamental mode in air-guiding photonic-bandgap fibers,” Opt. Express 13, 6669–6684 (2005). [CrossRef] [PubMed]

] This result points to another reason for avoiding surface modes: they increase backreflection of the fundamental (and other) core modes.

The minimum backreflection predicted for a single-mode PBF is determined to occur for a radius of R = 0.95 Λ, and is approximately equal to 2.9×10-6. Further reduction in the backreflection can be achieved by selecting a radius R = 1.9Λ, where the backreflection is predicted to be 4.2×10-7, but at the expense of the PBF no longer being single-mode.

Fig. 10. Variation of the fundamental mode reflection with PBF core radius, for a PBF cladding with ρ = 0.47Λ, at a wavelength of λ = 0.6Λ.

6. Effect of the termination medium refractive index on the reflection coefficient

While the reflections of the core modes of a PBF terminated in air are found to be low, much stronger reflections are expected for a solid or liquid termination medium (such as silica or lithium niobate, for example), because the core modes then propagate into a medium with a much greater index than their own. Since the Fresnel approximation is more accurate for a large index contrast between waveguide and termination medium, it is also expected to give better estimates of the reflection coefficients.

To investigate the influence of the termination medium index, we show in Fig. 11 the reflection of the fundamental mode for a PBF with R = 0.8Λ and ρ = 0.47Λ, and a wavelength λ = 0.6Λ, as a function of the termination refractive index. The reflections predicted by the Fresnel approximation are shown as a red curve. While both methods are in agreement for large termination medium indices, there are as expected significant differences for a termination medium index n 0 close to the refractive index of air. When n 0 is close to 1, the Fresnel approximation ceases to be valid and yields inaccurate answers, as we have seen before, by as much as 3 orders of magnitude. One interesting feature of these curves is the location of their minima. The Fresnel approximation predicts a minimum reflection for a termination medium of refractive index equal to neff < 1, while the calculated reflection exhibits a minimum for a termination medium index of approximately 1.005.

For larger values of the termination medium index, the agreement between both approaches is much improved. Notably, the calculated reflection coefficient for a silica termination is found to be 3.13%, differing by 8% with the value of 3.37% expected for a plane wave at normal incidence to an air-to-silica interface. The Fresnel approximation predicts a reflection of 3.67% for the fundamental mode at the same air-silica interface, within 17% of the computed value. This reflection coefficient is also approximately the reflection we expect from the butt-coupling of this PBF to a conventional index-guiding silica fiber, such as an SMF28 fiber.

Fig. 11. Variation of the reflection coefficient with termination medium refractive index, for a PBF structure with ρ = 0.47Λ, R = 0.8Λ, at a wavelength λ = 0.6Λ.

7. Simplified formulation

nmodal=n(r)ET(r)2dSET(r)2dS,
(6)

R=(nmodaln0nmodal+n0)2,
(7)

Fig. 12. Plot of the averaged modal index of the fundamental mode as a function of normalized core radius for a PBF with ρ = 0.47Λ, at a wavelength of λ = 0.6Λ.

Fig. 13. Variation of the calculated reflection (blue diamonds) and reflection predicted by the averaged modal index approximation (red squares) of the fundamental mode as a function of the core radius for a PBF with ρ = 0.47Λ, at a wavelength λ = 0.6Λ.

It is therefore of interest to consider PBFs where the angular spread of the fundamental mode is much reduced, so as to have a mode behaving nearly like a plane wave under normal incidence at a PBF termination. A larger air-filling ratio than those considered heretofore (around 85 %) is required for that purpose, and we elected to model the HC-1550-2 fiber by Crystal Fibre, [7

7. Crystal Fibre website, http://www.crystal-fibre.com

] which has an air-filling ratio of over 90%. Figure 14 shows a scanning electron microscope picture of the PBF structure, with the entire photonic crystal cladding visible (8 layers of holes). We modeled the PBF using the refractive index structure shown in part in Fig. 15, using 5 layers of holes and periodic boundary conditions.

Fig. 14. Scanning electron microscope picture of the cross-section of the HC-1550-2 PBF by Crystal Fibre. [7]
Fig. 15. Refractive index profile cross-section of the PBF model used to approximate the fiber HC-1550-2. Grey areas mark air (refractive index 1) and black areas mark silica (refractive index chosen equal to 1.45).

The longitudinal Poynting vector component of the fundamental mode at a wavelength of λ = 0.4Λ is shown in Fig. 16. There is very strong confinement of the mode field into the air core, with a field intensity overlap with silica now less than 1%. The Fourier transform intensity of this same mode's transverse electric field, plotted in logarithmic scale in the wavevector domain, is shown in Fig. 17. The black circle shows the boundary for all wavevectors such that kx2+ky2k02=(2πλ)2, i.e. all wavevectors admissible for plane waves propagating in air. The energy of the fundamental mode shows a much smaller angular spread, of the order of 3°. The calculated value for the fundamental mode reflection coefficient is 1.9×10-6 for an air-termination.

Fig. 16. Normalized plot of the longitudinal Poynting vector component for the y-polarized fundamental air-guided mode of the PBF HC-1550-2 model, at a wavelength of λ = 0.4Λ.
Fig. 17. Fourier transform intensity, in logarithmic scale, of the transverse electric field of the fundamental mode plotted in Fig. 16. The black circle marks all wavevectors such that k 2 x + k 2 yk 2 0, corresponding to plane waves that can be excited in air.

Fig. 18. Dependence of the reflection coefficient on termination medium refractive index, calculated with two different models. The PBF refractive index profile is shown in Fig. 15.

8. Conclusion

We have used an eigenmode decomposition method to compute the backreflection coefficient of various modes in PBFs, for terminations in a uniform medium, since the assumptions on which the commonly used Fresnel approximation relies cease to be valid for air-core fibers. The predicted values range from 10-5 to 10-6 for the fundamental mode, a few 10-3 for surface modes, and around 10-2 for bulk modes. Furthermore we observed strong variations of the fundamental mode reflection as a function of core radius, due to the presence of surface modes, which enhance the overlap of the fundamental mode with the silica. For the Crystal Fibre HC-1550-2 PBF we predicted a reflection coefficient of 1.9×10-6 for an air termination, and 3.3% for a silica termination, which would correspond to the splice reflection of a PBF with an index-guiding conventional silica optical fiber.

Appendix: Eigenmode decomposition method

ER(r)=m=0NMamEm(r),
(A1)

where a 0 is the amplitude reflection coefficient of the incident mode considered into itself. The modes present in the decomposition of the reflected field are of the same symmetry class as the fundamental mode. The intensity reflection coefficient is:

R=a02,
(A2)

In the same fashion, the transverse transmitted electrical field propagating in region II in the positive z direction is denoted ET(r) and decomposed in a sum of plane waves propagating in free space:

ET(r)=kΩ(bk1EW,k1+bk2EW,k2)eikrS,
(A3)

where Ω is the set of plane waves used to decompose ET, s is a normalization constant, and E1W,k and E2W,k denote the two amplitudes along orthogonal polarizations of a plane wave propagating in air, with a wavevector k=[kx ky kz] satisfying:

k2=kx2+ky2+kz2k02=(2πλ)2,
(A4)

The transverse magnetic fields are decomposed in the same fashion:

HI(r)=H0(r),
(A5)
HR(r)=amm=0NMHm(r),
(A6)
HT(r)=kΩ(bk1HW,k1+bk2HW,k2)eikrS,
(A7)

Using the continuity of the transverse electric and magnetic fields, we obtain the following relationships:

EI+ER=ET
E0(r)+m=0NMamEm(r)=kΩ(bk1EW,k1+bk2EW,k2)eikrS,
(A8)
HI+HR=HT
H0(r)+m=0NMamHm(r)=kΩ(bk1HW,k1+bk2HW,k2)eikrS,
(A9)

and we define a hermitian product over the function space by:

AB=2A*(r)B(r)dS,
(A10)

where the integral is calculated across the transverse plane.

We can then exploit the orthogonality of the transverse fields of the waveguide modes, and the orthogonality of the transverse fields of the plane waves, to derive the reflection coefficients am and transmission coefficients bk:

E02(1+a0)=k∊Ω(bk1E0Ew,k1eikrs+bk2E0Ew,k2eikrs),
(A11)
Em2am=k∊Ω(bk1EmEw,k1eikrs+bk2EmEw,k2eikrs),m1,
(A12)
Hw,k1eikrs2bk1=(1a0)Hw,k1eikrsH0m=1NMamHw,k1eikrsHm,
(A13)
Hw,k2eikrs2bk2=(1a0)Hw,k2eikrsH0m=1NMamHw,k2eikrsHm,
(A14)

yielding (1+NM+2NW) linear equations to compute the (NM+1) am and 2NW bk coefficients. Equations (A11) and (A12) are derived from Eq. (A8), and Eq. (A13) and (A14) are derived from Eq. (A9).

These linear equations are encoded in a sparse matrix A of size (1+NM+2NW)×(1+NM+2NW), containing 1+NM+2NW+4(1+NM)NW nonzero terms. The reflection coefficients ν are determined by solving an equation of the form:

Aν=w,
(A15)

which is solved by performing an incomplete LU decomposition of A, which is used as a preconditioner to lower the conditioning number of A, and we thus solve:

(LU)1Aν=(LU)1w,
(A16)

which allows for faster and more accurate solving of the reflection coefficient vector through the numerical method of the stabilized biconjugate gradients method BICGSTAB,[11

11. V. Dangui, M. J. F. Digonnet, and G. S. Kino, “A fast and accurate numerical tool to model the modal properties of photonic-bandgap fibers,” Opt. Express 14, 2979–2993 (2006). [CrossRef] [PubMed]

] with an approximate solution ν0 with an error measured by:

Aν0ww=ε1012,
(A17)

thus ensuring the very high accuracy of the solution within a few iterations of BICGSTAB (fewer than 10 usually), in a few seconds of calculation.

Acknowledgments

This work was supported by Litton Systems, Inc., a wholly owned subsidiary of Northrop Grumman Corporation.

References and links

1.

J. Broeng, D. Mogilevstev, S. E. Barkou, and A. Bjarklev, “Photonic crystal fibers: A new class of optical waveguides,” Opt. Fiber Technol. 5, 305–330 (1999). [CrossRef]

2.

J. C. Knight, T. A. Birks, R. F. Cregan, P. St. J. Russell, and J. P. Sandro, “Photonic crystals as optical fibers-physics and applications,” Opt. Mater. 11, 143–151 (1999). [CrossRef]

3.

R. S. Windeler, J. L. Wagener, and D. J. Giovanni, “Silica-air microstructured fibers: Properties and applications,” Optical Fiber Communications Conference, San Diego, CA, USA (1999).

4.

V. Dangui, H. K. Kim, M. J. F. Digonnet, and G. S. Kino, “Phase sensitivity to temperature of the fundamental mode in air-guiding photonic-bandgap fibers,” Opt. Express 13, 6669–6684 (2005). [CrossRef] [PubMed]

5.

D. M. Dagenais, K. P. Koo, and F. Bucholtz, “Fiber interferometry limitations due to parasitic optical cavities,” Lasers and Electro-Optics Society Conference, San Jose, CA USA (1993).

6.

J. Corbett, A. Dabirian, T. Butterley, N. A. Mortensen, and J. R. Allington-Smith, “The coupling performance of photonic crystal fibres in fibre stellar interferometry,” Mon. Not. Roy. Astron. Soc. 368, 203–210 (2006).

7.

Crystal Fibre website, http://www.crystal-fibre.com

8.

D. Khalil, “Reflection at the end of strongly guiding dielectric waveguide,” Opto-Electronics and Communications Conference, Yokohama, Japan , 1, 352–353 (2002).

9.

P. Gerard, P. Benech, D. Khalil, R. Rimet, and S. Tedjini, “Towards a full vectorial and modal technique for the analysis of integrated optics structures: The Radiation Spectrum Method (RSM),” Opt. Commun. 140, 128–145 (1997). [CrossRef]

10.

D. Marcuse, Theory of dielectric optical waveguides, (Academic Press, 1974).

11.

V. Dangui, M. J. F. Digonnet, and G. S. Kino, “A fast and accurate numerical tool to model the modal properties of photonic-bandgap fibers,” Opt. Express 14, 2979–2993 (2006). [CrossRef] [PubMed]

12.

H. . van der Vorst, “BI-CGSTAB: A fast and smoothly converging variant of BI-CG for the solution of nonsymmetric linear systems,” SIAM J. Sci. Stat. Comput. 13, 631–644 (1992). [CrossRef]

13.

H. K. Kim, J. Shin, S. H. Fan, M. J. F. Digonnet, and G. S. Kino, “Designing air-core photonic-bandgap fibers free of surface modes,” IEEE J. Quantum Eectron. 40, 551–556 (2004). [CrossRef]

14.

D. C. Allan, N. F. Borrelli, M. T. Gallagher, D. Müller, C. M. Smith, N. Venkataraman, J. A. West, P. Zhang, and K. W. Koch, “Surface modes and loss in air-core photonic band-gap fibers,” in Photonic Crystal Materials and Devices, A. Adibi, A. Scherer, and S. Yu Lin, eds., Proc. SPIE 5000,161–174 (2003). [CrossRef]

15.

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

16.

J. A. West, C. M. Smith, N. F. Borrelli, D. C. Allan, and K. W. Koch, “Surface modes in air-core photonic band-gap fibers,” Opt. Express 12, 1485–1496 (2004). [CrossRef] [PubMed]

OCIS Codes
(000.4430) General : Numerical approximation and analysis
(060.2310) Fiber optics and optical communications : Fiber optics
(120.5700) Instrumentation, measurement, and metrology : Reflection

ToC Category:
Photonic Crystal Fibers

History
Original Manuscript: August 16, 2006
Revised Manuscript: November 6, 2006
Manuscript Accepted: November 7, 2006
Published: April 18, 2007

Citation
Vinayak Dangui, Michel J. Digonnet, and Gordon S. Kino, "Determination of the mode reflection coefficient in air-core photonic bandgap fibers," Opt. Express 15, 5342-5359 (2007)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-15-9-5342


Sort:  Author  |  Year  |  Journal  |  Reset  

References

  1. J. Broeng, D. Mogilevstev, S. E. Barkou, and A. Bjarklev, "Photonic crystal fibers: A new class of optical waveguides," Opt. Fiber Technol. 5, 305-330 (1999). [CrossRef]
  2. J. C. Knight, T. A. Birks, R. F. Cregan, P. St. J. Russell, and J. P. Sandro, "Photonic crystals as optical fibers-physics and applications," Opt. Mater. 11, 143-151 (1999). [CrossRef]
  3. R. S. Windeler, J. L. Wagener, and D. J. Giovanni, "Silica-air microstructured fibers: Properties and applications," Optical Fiber Communications Conference, San Diego, CA, USA (1999).
  4. V. Dangui, H. K. Kim, M. J. F. Digonnet, and G. S. Kino, "Phase sensitivity to temperature of the fundamental mode in air-guiding photonic-bandgap fibers," Opt. Express 13, 6669-6684 (2005).H [CrossRef] [PubMed]
  5. D. M. Dagenais, K. P. Koo, and F. Bucholtz, "Fiber interferometry limitations due to parasitic optical cavities," Lasers and Electro-Optics Society Conference, San Jose, CA, USA (1993).
  6. J. Corbett, A. Dabirian, T. Butterley, N. A. Mortensen, and J. R. Allington-Smith, "The coupling performance of photonic crystal fibres in fibre stellar interferometry," Mon. Not. Roy. Astron. Soc. 368, 203-210 (2006).
  7. Crystal Fibre website, http://www.crystal-fibre.com
  8. D. Khalil, "Reflection at the end of strongly guiding dielectric waveguide," Opto-Electronics and Communications Conference, Yokohama, Japan, 1, 352-353 (2002).
  9. P. Gerard, P. Benech, D. Khalil, R. Rimet, and S. Tedjini, "Towards a full vectorial and modal technique for the analysis of integrated optics structures: The Radiation Spectrum Method (RSM)," Opt. Commun. 140, 128-145 (1997). [CrossRef]
  10. D. Marcuse, Theory of dielectric optical waveguides, (Academic Press, 1974).
  11. V. Dangui, M. J. F. Digonnet, and G. S. Kino, "A fast and accurate numerical tool to model the modal properties of photonic-bandgap fibers," Opt. Express 14, 2979-2993 (2006).H [CrossRef] [PubMed]
  12. H. . van der Vorst, "BI-CGSTAB: A fast and smoothly converging variant of BI-CG for the solution of nonsymmetric linear systems," SIAM J. Sci. Stat. Comput. 13, 631-644 (1992). [CrossRef]
  13. H. K. Kim, J. Shin, S. H. Fan, M. J. F. Digonnet, and G. S. Kino, "Designing air-core photonic-bandgap fibers free of surface modes," IEEE J. Quantum Eectron. 40, 551-556 (2004). [CrossRef]
  14. D. C. Allan, N. F. Borrelli, M. T. Gallagher, D. Müller, C. M. Smith, N. Venkataraman, J. A. West, P. Zhang, and K. W. Koch, "Surface modes and loss in air-core photonic band-gap fibers," in Photonic Crystal Materials and Devices, A. Adibi, A. Scherer, S. Yu Lin, eds., Proc. SPIE 5000,161-174 (2003). [CrossRef]
  15. 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). [CrossRef] [PubMed]
  16. J. A. West, C. M. Smith, N. F. Borrelli, D. C. Allan, and K. W. Koch, "Surface modes in air-core photonic band-gap fibers," Opt. Express 12, 1485-1496 (2004). [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