OSA's Digital Library

Optics Express

Optics Express

  • Editor: C. Martijn de Sterke
  • Vol. 20, Iss. 19 — Sep. 10, 2012
  • pp: 20980–20991
« Show journal navigation

Analysis of light scattering from surface roughness in hollow-core photonic bandgap fibers

Eric Numkam Fokoua, Francesco Poletti, and David J. Richardson  »View Author Affiliations


Optics Express, Vol. 20, Issue 19, pp. 20980-20991 (2012)
http://dx.doi.org/10.1364/OE.20.020980


View Full Text Article

Acrobat PDF (8772 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

We present a theoretical method for analyzing radiation loss from surface roughness scattering in hollow-core photonic bandgap fibers (HC-PBGFs). We treat the scattering process as induced dipole radiation and combine statistical information about surface roughness, mode field distribution and fibre geometry to accurately describe the far-field scattering distribution and loss in fibers with an arbitrary cross-sectional distribution of air holes of any shape. The predicted angular scattering distribution, total scattering loss and the loss wavelength dependence are all shown to agree well with reported experimental data. Our method yields a simpler result than that obtained by more complex approaches and is to the best of our knowledge the first successful attempt to accurately describe roughness scattering in HC-PBGFs.

© 2012 OSA

1. Introduction

However, the lowest loss value reported to date for HC-PBGFs stands at 1.7dB/km [5

5. B. J. Mangan, L. Farr, A. Langford, P. J. Roberts, D. P. Williams, F. Couny, M. Lawman, M. Mason, S. Coupland, R. Flea, H. Sabert, T. A. Birks, J. C. Knight, and P. St. J. Russell, “Low loss (1.7 dB/km) hollow core photonic bandgap fiber,” in Proceedings of Optical Fiber Communication Conference (2004), paper PDP24.

], nearly an order of magnitude higher than in typical single mode fibers. Indeed, HC-PBGFs suffer from different loss mechanisms. Though the finite extent of the periodic cladding implies that the guided modes are always leaky, the leakage loss can be reduced to negligible levels by incorporating an appropriate number of rings of air holes in the cladding. This has led to the identification of scattering from the many air-glass interfaces present within HC-PBGFs as the fundamental mechanism dominating their loss [6

6. P. J. Roberts, F. Couny, H. Sabert, B. J. Mangan, D. P. Williams, L. Farr, M. W. Mason, A. Tomlinson, T. A. Birks, J. C. Knight, and P. S. J. Russell, “Ultimate low loss of hollow-core photonic crystal fibres,” Opt. Express 13236–244 (2005). [CrossRef] [PubMed]

]. The intrinsic roughness on these surfaces is due to thermally excited surface capillary waves (SCWs) present at the surface of molten glass which freeze-in as the glass goes through phase transition during the fibre draw. The thermodynamic origin of this roughness makes the scattering loss mechanism as fundamental as Rayleigh scattering in solid fibers. However, unlike Rayleigh scattering which is solely determined by the material choice, the surface scattering also crucially depends on the specific fiber design. Enlarging the air core or moving the photonic bandgap to longer wavelengths for example are known ways to reduce the roughness scattering loss [4

4. M. N. Petrovich, F. Poletti, A. van Brakel, and D. J. Richardson, “Robustly single mode hollow core photonic bandgap fiber,” Opt. Express 16, 4337–4346 (2008). [CrossRef] [PubMed]

, 6

6. P. J. Roberts, F. Couny, H. Sabert, B. J. Mangan, D. P. Williams, L. Farr, M. W. Mason, A. Tomlinson, T. A. Birks, J. C. Knight, and P. S. J. Russell, “Ultimate low loss of hollow-core photonic crystal fibres,” Opt. Express 13236–244 (2005). [CrossRef] [PubMed]

].

Paramount to obtaining further loss reduction in HC-PBGFs is the development of an accurate mathematical description of the scattering process. The standard theoretical treatment of scattering from waveguide imperfections considers the scattering as coupling between the guided modes and all the radiation modes of the structure [7

7. D. Marcuse, Theory of Dielectric Optical Waveguides, 2nd ed. (Academic Press, 1991).

, 8

8. A. W. Snyder and J. D. Love, Optical Waveguide Theory (Chapman and Hall, 1983).

]. Roberts et al. [9

9. P. J. Roberts, F. Couny, H. Sabert, B. J. Mangan, T. A. Birks, J. C. Knight, and P. S. Russell, “Loss in solid-core photonic crystal fibers due to interface roughness scattering,” Opt. Express 13, 7779–7793 (2005). [CrossRef] [PubMed]

] have applied such an approach to solid-core photonic crystal fibers, using a method involving complicated Green tensor math which in practice limits its applicability to structures with circular air holes. Phan-Huy et al [10

10. M.-C. Phan-Huy, J.-M. Moison, J. A. Levenson, S. Richard, G. Melin, M. Douay, and Y. Quiquempois, “Surface roughness and light scattering in a small effective area microstructured fiber,” J. Lightwave. Technol. 27, 1597–1604 (2009). [CrossRef]

] have treated the same problem in solid-core PCFs by breaking down the process into coupling to higher order core modes and subsequent leakage, but their approach ignores the direct coupling to radiation modes. Radiation modes were also neglected when Dangui et al. [11

11. V. Dangui, M. J. F. Digonnet, and G. S. Kino, “Modeling of the propagation loss and backscattering in air-core photonic-bandgap fibers,” J. Lightwave Technol. 17, 3783–3789 (2009). [CrossRef]

] numerically implemented the coupled-mode theory for HC-PBGFs, with their result requiring the computation of an impractically large number of coupling coefficients for several thousand cladding modes.

Here we extend an approach pioneered by Rawson [12

12. E. G. Rawson, “Theory Of scattering by finite dielectric needles illuminated parallel to their axes,” J. Opt. Soc. Am. 62, 1284–1286 (1972). [CrossRef]

, 13

13. E. G. Rawson, “Analysis of scattering from fiber waveguides with irregular core surfaces,” Appl. Opt. 13, 2370–2377 (1974). [CrossRef] [PubMed]

] and later used by Mazumder et al. [14

14. P. Mazumder, S. L. Logunov, and S. Raghavan, “Analysis of excess scattering in optical fibers,” J. Appl. Phys. 96, 4042–4049 (2004). [CrossRef]

] for weakly guiding fibers and derive a simpler method in which the far-field distribution of the scattered light is calculated using dipole radiation theory. Our method is suitable for describing scattering from boundary imperfections in a broad range of optical waveguides and in particular HC-PBGFs with air holes of arbitrary shape. The method predicts both an angular distribution of scattered light and overall attenuation values that agree well with experiments. We also show that the predicted wavelength dependence of the scattering loss agrees well with reported experimental data. Finally, we highlight the strong dependence of the loss on the roughness power spectral density, by showing that surfaces with the same value of root mean square (rms) roughness but different power density spectra can yield very different attenuation values.

2. Formulation

As shown in Fig. 1 the optical field of air-guided modes in HC-PBGFs always has some unavoidable overlap with the surrounding glass surfaces which generates scattered radiation and hence loss since these surfaces are rough. The intrinsic roughness causes all the air-glass interfaces to depart from those of an ideal fiber. If we consider a single air hole, its boundary imperfections can be described by a random function f (s, z) where s is the curvilinear coordinate along the hole perimeter.

Fig. 1 Overlap between Ex component of the field of the fundamental guided mode and glass surfaces in HC-PBGFs with an illustration of the intrinsic roughness present at air-glass interfaces. The contour line are 2dB apart down to 30dB lower than the maximum value of the electric field at the center of the core.

Our aim in this work is to calculate the far-field distribution of the light scattered by the collection of all the dipoles in a section of length 2L of the perturbed waveguide. This section is chosen so that it contains all the roughness spectral components of interest, but is also short enough for the incident field E0 to be approximately constant throughout. At point P(ϑ, ϕ) (see Fig. 2) on a distant sphere of radius R, the scattered field element from a single dipole located at the origin is [16

16. J. D. Jackson, Classical Electodynamics, 3rd ed. (John Wiley and Sons, 1998)

]:
dEs(ϑ,ϕ)=k024πε0R(r×p)×r
(3)
where k0 is the free-space propagation constant. The total field at P is obtained by adding together all the scattered field elements while taking into account the relative phases. Taking as reference a ray that would be scattered from the origin, the phase difference for a ray scattered from position (s, z) can be calculated with the help of Fig. 2 as:
Φ(ϑ,ϕ,s,z)=βzk0zcosϕk0ysinϕ
(4)

Fig. 2 (a) Illustration of scattering point P on a distant sphere and (b) the relative phase between a ray scattered by the dipole at (si, z) and a reference ray.

In this expression, βz represents the phase build up of the incident mode field between z = 0 and z, while k0zcosϕk0ysinϕ is the phase accumulated by the reference ray in the same time. It has been assumed that all the scattered rays propagate through air, which is a reasonable approximation for the high air-filling fraction periodic cladding of HC-PBGFs. y′ = ycosϑxsinϑ is the new y coordinate of the radiating dipole in the frame obtained by rotating the original xy axes by an angle ϑ around the z-axis. Combining Eq. (3) and (4), the total scattered field at P from a single air hole is
Es(ϑ,ϕ)=z=LLdEs(ϑ,ϕ)exp(iΦ)
(5)
where the closed path integral is along the hole contour. In the case of multiple air holes such as in HC-PBGFs, we assume that the random distortions on each surface are identically distributed but uncorrelated. Using j as an index to label each air hole, the total scattered field is therefore obtained by a summation over all surfaces, which with the help of Eq. (3) and (2) becomes:
Es(ϑ,ϕ)=jz=LLdEs(ϑ,ϕ)exp(iΦ)=jk02α04πε0Rz=LLf(sj,z)(r×E0)×rexp(iΦ)dsjdz
(6)

Equation (6) contains a very general form of the boundary distortion. In simple structures, when the roughness can be separated into axial and azimuthal parts, the solution to Eq. (6) follows along the lines of the work by Mazumder et al. [14

14. P. Mazumder, S. L. Logunov, and S. Raghavan, “Analysis of excess scattering in optical fibers,” J. Appl. Phys. 96, 4042–4049 (2004). [CrossRef]

]. In the case of roughness due to frozen-in capillary waves at glass-air interfaces as in HC-PBGFs that we discuss in section 4.1, the most relevant features of the roughness and those that dictate the angular distribution of scattered light are those along the fiber axis. We therefore restrict f (si, z) to f (z) only. The z–dependent part of the above expression can be rewritten as:
LLf(z)ei(βk0cosϕ)zdzF˜(βk0cosϕ)
(7)
where F̃ is the spatial Fourier transform of the roughness distribution. Expanding the incident field at the scattering interface in terms of its components perpendicular (E0⊥) and parallel(E0||) to the Pz plane and carrying out the vector product, Eq. (6) becomes:
Es(ϑ,ϕ)=k02α04πε0RF˜(βk0cosϕ)jeik0ysinϕ[E0v+(E0zsinϕE0||cosϕ)u]dsj=k02α04πε0RF˜(βk0cosϕ)×[Esuu+Esvv]
(8)
where u and v are the unit vectors shown in Fig. 2. This result shows that the far-field distribution of scattered light is the product of a component depending on the roughness spectrum and a quantity that depends on the incident electric field at the interface, the geometry and optical properties of the fibre. Indeed, each spatial frequency component of the roughness scatters light in a specific direction ϕ, a fact that conforms with the more complicated analysis of coupled-mode theory or volume-current methods [7

7. D. Marcuse, Theory of Dielectric Optical Waveguides, 2nd ed. (Academic Press, 1991).

, 8

8. A. W. Snyder and J. D. Love, Optical Waveguide Theory (Chapman and Hall, 1983).

]. This also implies that only the spatial frequencies satisfying βk0κβ + k0 contribute to radiation loss.

Equation (8) contains all the details of the scattered field but often the main quantity of interest is the average distribution of scattered power such as that typically measured in angularly resolved scattering (ARS) experiments:
P(ϕ)=12πϑ=02π12cnε0|Es(ϑ,ϕ)|2R2sinϕdϑ.
(9)
To describe the statistical properties of the surface roughness, the roughness auto-correlation R(u) and the power spectral density (PSD) which are a Fourier transform pair need to be introduced. The power spectral density of the roughness is defined as [26

26. J. M. Elson and J. M. Bennett, “Relation between the angular dependence of scattering and the statistical properties of optical surfaces,” J. Opt. Soc. Am. 6931–47 (1979). [CrossRef]

]
S(κ)=limL[12L|LLf(z)eiκzdz|2]=12L|F˜(κ)|2,
(10)
where 〈...〉 is the ensemble average over a large set of surfaces with equivalent statistics. By inserting Eq. (8) into Eq. (9), we obtain with the help of Eq. (10):
P(ϕ)=2L2π(k02α04πε0)2S(βk0cosϕ)×ϑ=02π12cnε0[|Esu|2+|Esv|2]sinϕdϑ.
(11)
To account for the fact that some of the scattered light may be recaptured by the waveguide for example if emitted at an angle falling within the total internal reflection or the the photonic bandgap acceptance cones, we denote by L(ϕ) the fraction of light scattered in the polar direction ϕ that is actually lost. If P0 is the total optical power carried by the incident guided mode, the exponential loss coefficient is obtained by integration as:
αsc2π12LP0ϕ=0πP(ϕ)L(ϕ)dϕ.
(12)

Note that in deriving Eq. (11), we have assumed that the scattered light propagates in air, which is a reasonable approximation for within the microstructured area. In all HC-PBGFs however, the microstructure is surrounded by a silica outer cladding, and in ARS experiments an index-matched fluid is often used between the fiber and an external detector. In this case, the scattered light is refracted at the interface between the microstructured cladding and the glass outer jacket with refractive index ngl. The measured scattering angle ϕm is therefore obtained through Snell’s law ngl cos ϕm = cos ϕ.

3. Application to single-mode fibers

For a first verification of the accuracy of our model, we consider the classical problem of radiation loss due to scattering from a sinusoidal diameter variation in a standard single-mode fiber [17

17. A. W. Snyder, “Radiation losses due to variations of radius on dielectric or optical fibers,” IEEE Trans. Microwave Theory Tech. MT18(9) 608–615 (1970). [CrossRef]

, 18

18. D. Marcuse, “Radiation losses of HE11 mode of a fiber with sinusoidally perturbed core boundary,” Appl. Opt. 14, 3021–3025 (1975). [CrossRef] [PubMed]

]. The structure which is shown in the inset of Fig. 3 is a fibre with a core with radius r = a + bsin(2πz/Λ) and refractive index n1, and a cladding with refractive index n2. When the sinusoidal perturbation is small enough to satisfy the assumptions in our treatment (bλ), the scattering loss suffered by the x–polarized fundamental HE11 mode can be calculated using equation (11). Since the cladding is no longer air, the relative phase difference of equation (4) needs to be adjusted by replacing k0 with k0n2. For a pure sinusoidal perturbation over a long fiber section (L ≫ Λ), the power spectral density is
S(κ)=b24δ(κ2πΛ)
(13)
Thus, not surprisingly, scattering only occurs in the direction that satisfies the phase-matching condition:
βk0n2cosϕ2πΛ=0cosϕ=β2π/Λk0n2.
(14)
Rather than showing the radiation pattern which is sharply peaked at this angle due to the delta PSD, it is more interesting to plot the loss coefficient as a function of the escaping angle when the spatial wavelength of the perturbation is changed. Such an example is shown in Fig. 3 where we have used n1 = 1.46, n2 = 1.458, a = 5μm and the wavelength of λ = 1μm is such that the waveguide is just below the cut-off of the first higher-order mode (V = 2.4). To show the compatibility of our method with common numerical tools, the field components of the fundamental mode were calculated using a fully vectorial finite element method, and equation (8) was evaluated to obtain the scattered field. With L(ϕ) = 1 for all angles larger than the critical angle and L(ϕ) = 0 otherwise, the normalized scattering loss coefficient shown in Fig. 3 is identical to that obtained by Marcuse using a far more complex coupled mode theory [7

7. D. Marcuse, Theory of Dielectric Optical Waveguides, 2nd ed. (Academic Press, 1991).

]. This not only validates our simple approach but also shows that it can be implemented in combination with any arbitrary mode-solving numerical tool.

Fig. 3 Normalized exponential scattering loss coefficient of the fundamental LP01 mode due to sinusoidal core diameter variations of a step-index fiber with the parameters shown in the inset.. The loss is plotted as a function of the escaping angle ϕ which is given by ϕ = cos−1((β − 2π/Λ)/n2k0). The result plotted here and obtained by solving Eq. (11) is identical to that obtained by coupled-mode theory (see [7], page 159)

4. Roughness scattering in HC-PBGFs

4.1. Surface roughness due to frozen-in surface capillary waves

As already mentioned, thermally excited surface capillary waves impose an intrinsic roughness on glass surfaces. On a flat two-dimensional surface, they lead to a two-dimensional power spectral density given by the well-known expression [20

20. J. Jäckle and K. Kawasaki, “Intrinsic roughness of glass surfaces,” J. Phys. Condens. Matter 7, 4351–4358 (1995). [CrossRef]

]:
S(κ)=kBTgρg+γ|κ|2,
(15)
where κ is the two dimensional surface wavevector, ρ the glass density and g the gravity constant. kB represents Boltzmann’s constant, Tg the glass transition temperature and γ the surface tension of the glass, which for pure silica are ∼ 1500K and ∼ 1J.m−2 respectively.

In the case of HC-PBGFs and holey fibers in which the glass surfaces are cylindrical rather than flat, the roughness PSD becomes periodic along the azimuthal direction and hence can be expanded into an infinite Fourier series. Neglecting the gravity contribution (since the fibers are made under an imposed vertical flow), integrating over all the azimuthal wavevectors leads to the 1-D form given by [6

6. P. J. Roberts, F. Couny, H. Sabert, B. J. Mangan, D. P. Williams, L. Farr, M. W. Mason, A. Tomlinson, T. A. Birks, J. C. Knight, and P. S. J. Russell, “Ultimate low loss of hollow-core photonic crystal fibres,” Opt. Express 13236–244 (2005). [CrossRef] [PubMed]

]:
Sz(κ)=kBTg4πγκ.
(16)

The 2-D expression of the roughness PSD in Eq. (15) requires an upper wavevector cut-off |κu| for the surface roughness to be physically acceptable, and gravity imposes another lower wavevector cut-off determined by the capillary length, which is of the order of ∼ 4mm [21

21. T. Sarlat, A. Lelarge, E Søndergård, and D. Vandembroucq, “Frozen capillary waves on glass surfaces: an AFM study,” Eur. Phys. J. B 54, 121–126 (2006). [CrossRef]

]. Equation (16) implies a logarithmic roughness autocorrelation along z which also requires the existence of low and high spatial frequency cut-offs to be physically acceptable. Since the high frequency cut-off is estimated to be of the order of a few molecular lengths (∼ 0.5nm), it is of no particular concern here as the highest frequency β + k0 contributing to scattering falls well below this value. In HC-PBGFs, the origins of a low-frequency cut-off are not yet known. Atomic Force microscopy (AFM) measurements on the inner surface of air-holes in HC-PBGFs and holey fibers revealed a roughness PSD consistent with frozen-in surface capillary waves of Eq. (16) between κ = 0.2μm−1 and κ = 30μm−1 but could not measure the low spatial frequency end of the spectrum [6

6. P. J. Roberts, F. Couny, H. Sabert, B. J. Mangan, D. P. Williams, L. Farr, M. W. Mason, A. Tomlinson, T. A. Birks, J. C. Knight, and P. S. J. Russell, “Ultimate low loss of hollow-core photonic crystal fibres,” Opt. Express 13236–244 (2005). [CrossRef] [PubMed]

,10

10. M.-C. Phan-Huy, J.-M. Moison, J. A. Levenson, S. Richard, G. Melin, M. Douay, and Y. Quiquempois, “Surface roughness and light scattering in a small effective area microstructured fiber,” J. Lightwave. Technol. 27, 1597–1604 (2009). [CrossRef]

]. The results presented in this paper are therefore obtained by imposing that spatial frequencies below the arbitrary cut-off point κc = 0.1μm−1 have the same PSD value S(κc).

4.2. Angular distribution of scattered power and loss

Using expression (16) for the PSD of surface capillary waves and the cut-off previously mentioned, we calculate the scattered field distribution according to Eq. (11). The fibre used in our example is an ideal HC-PBGF with a structure closely matching that of a low-loss fibre we recently reported [22

22. N. V. Wheeler, M. N. Petrovich, R. Slavík, N. Baddela, E. Numkam Fokoua, J. R. Hayes, D. R. Gray, F. Poletti, and D. J. Richardson “Wide-bandwidth, low-loss, 19-cell hollow core photonic band gap fiber and its potential for low latency data transmission,” in Proceedings of Optical Fiber Communication Conference (2012), paper PDP5A.2.

], a 19c fiber with air holes similar to hexagons with rounded corners (such as shown in Fig. 1) and in which the pitch is Λ ∼ 4.4μm and d/Λ ∼ 0.975. Figure 4(a) shows plots of P(ϕm)/(2LP0) as a function of scattering angle ϕm for a number of guided modes within the fiber, calculated with the assumption that all scattered light is lost, or L(ϕ) = 1 for all ϕ. The curve for the fundamental mode shows good qualitative agreement with the measured ARS data from Roberts et al, referring to a not too dissimilar fiber design [6

6. P. J. Roberts, F. Couny, H. Sabert, B. J. Mangan, D. P. Williams, L. Farr, M. W. Mason, A. Tomlinson, T. A. Birks, J. C. Knight, and P. S. J. Russell, “Ultimate low loss of hollow-core photonic crystal fibres,” Opt. Express 13236–244 (2005). [CrossRef] [PubMed]

]. Because the roughness described by Eq. (16) has high PSD at small spatial frequencies, scattering occurs primarily in the forward direction, as can be seen. Our model goes one step further in allowing us to fully evaluate the differential scattering loss between the guided modes of the fiber, which can be seen to be non-negligible for these fibers.

Fig. 4 (a)Angular scattering distribution from frozen-in surface capillary waves for a few guided modes of the HC-PBGF in [22]. (b) Measured and calculated loss contributions for the fundamental guided mode. A cut-off at κc = 0.1μm−1 has been imposed on the roughness PSD. In performing the cutback loss measurement, care was taken to predominantly excite the fundamental mode only. s2 and time of flight measurements showed an extinction higher than 22dB for all the higher order modes at the output (see [22].

Figure 4(b) shows a plot of the calculated roughness scattering loss and confinement loss across the photonic bandgap for the fundamental mode of the same fibre, along with the measured attenuation spectrum. As can be seen, our model suggests that the total loss is dominated by the roughness scattering contribution near the center of the bandgap, while the bandgap edges are determined by confinement loss. The calculated scattering loss values are within a factor of 2 of the measured ones, which is an excellent agreement considering the simplicity of our method.

4.3. Wavelength dependence of scattering loss

A crucial factor in determining the ultimate loss performance of HC-PBGFs is its wavelength dependence. If the effects of material dispersion are rigidly ignored, the wavelength λc at which the lowest loss value occurs is proportionally scaled when an HC-PBGF with a given cross-section is drawn to different dimensions [23

23. J. D. Joannopoulos, S. G. Johnson, J. N. Winn, and R. D. Meade, Photonic Crystal: Molding the Flow of Light, 2nd ed. (Princeton University Press, 2008).

]. It has been shown experimentally that when a given fiber cross-section is drawn to different diameters, its minimum loss decreases as λc3 until a wavelength around 2μm where infrared absorption mechanisms begin to dominate. However, a theoretical understanding of this dependence has not hitherto been provided [6

6. P. J. Roberts, F. Couny, H. Sabert, B. J. Mangan, D. P. Williams, L. Farr, M. W. Mason, A. Tomlinson, T. A. Birks, J. C. Knight, and P. S. J. Russell, “Ultimate low loss of hollow-core photonic crystal fibres,” Opt. Express 13236–244 (2005). [CrossRef] [PubMed]

,24

24. F. Benabid and P. J. Roberts, “Linear and nonlinear optical properties of hollow core photonic crystal fiber,” J. Mod. Opt. 2, 87–124 (2011). [CrossRef]

]. Despite the explicit k04 in Eq. (11), our model confirms the experiments by predicting a scattering loss that indeed follows a λca trend, with a = 3.06 in the case of the low-loss HC-PBGF, as shown in Fig. 5. In general, a depends on both the specific design of the fibre cross-section and the roughness PSD, and it has been found to assume values between 2.5 and 3.5 for a range of fibers we studied. This confirms that the scattering loss and its wavelength dependence are very much dependent on the design parameters.

Fig. 5 Predicted wavelength dependence of the minimum scattering loss coefficient when the fibre in [22] is rigidly scaled to larger transverse cross-sections.

On the same figure, we also show how the normalized interface field intensity
F=(ε0μ0)12holeperimeters|E|2dscrosssectionE×H*dA
(17)
which is routinely used to compare fiber designs scales as a function of wavelength [6

6. P. J. Roberts, F. Couny, H. Sabert, B. J. Mangan, D. P. Williams, L. Farr, M. W. Mason, A. Tomlinson, T. A. Birks, J. C. Knight, and P. S. J. Russell, “Ultimate low loss of hollow-core photonic crystal fibres,” Opt. Express 13236–244 (2005). [CrossRef] [PubMed]

, 25

25. R. Amezcua-Correa, N. G. R. Broderick, M. N. Petrovich, F. Poletti, and D. J. Richardson, “Design of 7 and 19 cells core air-guiding photonic crystal fibers for low-loss, wide bandwidth and dispersion controlled operation,” Opt. Express 15, 17577–17586 (2007). [CrossRef] [PubMed]

]. Here, F was calibrated using measurements and simulations for the fibre in [22

22. N. V. Wheeler, M. N. Petrovich, R. Slavík, N. Baddela, E. Numkam Fokoua, J. R. Hayes, D. R. Gray, F. Poletti, and D. J. Richardson “Wide-bandwidth, low-loss, 19-cell hollow core photonic band gap fiber and its potential for low latency data transmission,” in Proceedings of Optical Fiber Communication Conference (2012), paper PDP5A.2.

] at a wavelength of 1.5μm so that F = 0.0174μm−1 corresponds to 3.5dB/km. As can be seen, this scalar quantity does not explain the observed wavelength dependence of the loss and is thus inaccurate in quantifying it. This can be understood by realizing that in addition to the fibre cross-section determining the field at the interfaces and how it changes with wavelength, when the operational wavelength is increased, the range of roughness spatial frequencies that contribute to the scattering also changes and becomes narrower. For example, at a wavelength of 1μm the range of spatial frequencies required to calculate the scattering in air at angles ϕ from 0 to 180° ranges from 0 to 12.55μm−1 while at a wavelength of 2μm, it goes from 0 to 6.275μm−1. Therefore only models that fully account for the roughness PSD like our method (and do not just use its rms value like the F–parameter) can predict the wavelength dependence of the loss.

5. The impact of roughness PSD

Fig. 6 (a)Scattering loss as a function of the spatial frequency κc of the SCW roughness PSD at a wavelength of 1.55μm. (b) Computed total loss across the photonic bandgap for the fibre in [22] for three distinct values of the roughness frequency cut-off.

Figure 6(b) shows the total loss across the photonic bandgap for three values of κc. A decrease in κc corresponds to a higher PSD at the short spatial frequencies responsible for stronger scattering at longer wavelengths. As a result, the loss at longer wavelengths becomes higher than at shorter wavelengths, as can be readily seen in the figure. Additionally, we note that the scattering loss is directly proportional to the ratio Tg/γ which we have assumed is 1500K/(J.m−2), an increase in this ratio will therefore shift the scattering loss curve upwards.

SCWs may not be the only mechanism generating surface roughness in a HC-PBGF. While they have been found to dominate the roughness spectrum in HC-PBGFs between the spatial frequencies of 0.2 and 30μm−1, roughness from other sources may be present at frequencies below those practically measurable. It is well known that most realistic surfaces have roughness statistics that present one or more of three basic components: a long-range waviness, short-range random roughness and periodicity [26

26. J. M. Elson and J. M. Bennett, “Relation between the angular dependence of scattering and the statistical properties of optical surfaces,” J. Opt. Soc. Am. 6931–47 (1979). [CrossRef]

]. These additional roughness components, corresponding to long range waviness and possibly originating from drawing parameter fluctuations during fiber fabrication, may ultimately have an effect on the overall loss [11

11. V. Dangui, M. J. F. Digonnet, and G. S. Kino, “Modeling of the propagation loss and backscattering in air-core photonic-bandgap fibers,” J. Lightwave Technol. 17, 3783–3789 (2009). [CrossRef]

].

To investigate this, we study the loss that would be produced in our HC-PBGF by roughness components with an exponential or Gaussian auto-correlation. The analytical expressions for these are R(u) = σ2exp(−|u|/Lc) and R(u)=σ2exp(u2/Lc2) respectively, with Lc being the correlation length and σ the corresponding rms roughness. The PSD arising from these auto-correlation functions are a Lorentzian S(κ)=2σ2Lc/(1+Lc2κ2) and a Gaussian S(κ)=σ2Lc/πexp(κ2Lc2/4) respectively.

Figure 7(a) shows the scattering loss contribution arising from such roughness components as a function of the correlation length Lc at wavelength λ = 1.55μm, assuming σ of 0.1nm. In both cases, the loss increases with Lc, peaks near Lc 4mm and then decreases quickly, as typically observed with exponential or Gaussian roughness [27

27. D. Marcuse, “Mode conversion caused by surface imperfection of a dielectric slab waveguide,” Bell Syst. Tech. J. 48, 3187–3215 (1969).

]. In a broad range of correlation lengths (70μm − 3.7cm) the scattering loss from such roughness components alone is higher than the measured value, which would suggest discarding this range. If however a long-range exponential roughness component with Lc of the order of a few centimeters is present in our HC-PBGF in addition to SCWs, this leads to a predicted scattering loss value very close to the measured 3.5dB/km. The contribution of this additional roughness component would be many orders of magnitude below that of SCWs over AFM-measurable spatial frequencies and would therefore be compatible with reported measurements, as can be seen in Fig. 7(b). Although clearly speculative, this confirms the importance of accurate roughness information in the low spatial frequency region to allow us making even more accurate loss predictions.

Fig. 7 (a)Scattering loss from exponential and Gaussian roughness as function of correlation length.(b) Roughness PSD for frozen-in SCWs, Gaussian and exponential roughness. If exponential or gaussian roughness with the correlation lengths shown are present in the fiber, their contribution is negligible in the region where SCWs dominate.

6. Conclusion

We have presented a simple theoretical method based on dipole radiation to describe light scattering from surface roughness in HC-PBGFs. Our method shows excellent agreement with the more standard but considerably more complex approach of coupled-mode theory. Our model not only returns an accurate value of loss, but also gives detailed information about the far-field distribution of scattered light, including the angular distribution of scattered power typically measured in angularly resolved scattering experiments. Using our model and a number of reasonable assumptions about the surface roughness, the predicted loss from frozen-in SCWs alone is within a factor of 2 for a state-of-the-art HC-PBGF. Realistic considerations on other possible origins of roughness such as a long-range waviness, especially in the spectral region not accessible by current roughness measurement methods, can lead to more accurate loss predictions. Our simple approach allows calculation of scattering losses from any microstructured fiber of arbitrary geometry for which the power spectral density of the roughness is known and it can be used in combination with any numerical method to calculate the field intensity of the guided mode at the interfaces.

Acknowledgments

This work was supported by the EU 7th Framework Programme under grant agreement 228033 ( MODE-GAP).

References and links

1.

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

2.

R. F. Cregan, B. J. Mangan, J. C. Knight, T. A. Birks, P. S. Russell, P. J. Roberts, and D. C. Allan, “Single-mode photonic band gap guidance of light in air,” Nature (London) 285, 1537–1539 (1999).

3.

C. M. Smith, N. Venkataraman, M. T. Gallagher, D. Muller, J. A. West, N. F. Borrelli, D. C. Allan, and K. W. Koch, “Low-loss hollow-core silica/air photonic bandgap fibre,” Nature (London) 424, 657–659 (2003). [CrossRef]

4.

M. N. Petrovich, F. Poletti, A. van Brakel, and D. J. Richardson, “Robustly single mode hollow core photonic bandgap fiber,” Opt. Express 16, 4337–4346 (2008). [CrossRef] [PubMed]

5.

B. J. Mangan, L. Farr, A. Langford, P. J. Roberts, D. P. Williams, F. Couny, M. Lawman, M. Mason, S. Coupland, R. Flea, H. Sabert, T. A. Birks, J. C. Knight, and P. St. J. Russell, “Low loss (1.7 dB/km) hollow core photonic bandgap fiber,” in Proceedings of Optical Fiber Communication Conference (2004), paper PDP24.

6.

P. J. Roberts, F. Couny, H. Sabert, B. J. Mangan, D. P. Williams, L. Farr, M. W. Mason, A. Tomlinson, T. A. Birks, J. C. Knight, and P. S. J. Russell, “Ultimate low loss of hollow-core photonic crystal fibres,” Opt. Express 13236–244 (2005). [CrossRef] [PubMed]

7.

D. Marcuse, Theory of Dielectric Optical Waveguides, 2nd ed. (Academic Press, 1991).

8.

A. W. Snyder and J. D. Love, Optical Waveguide Theory (Chapman and Hall, 1983).

9.

P. J. Roberts, F. Couny, H. Sabert, B. J. Mangan, T. A. Birks, J. C. Knight, and P. S. Russell, “Loss in solid-core photonic crystal fibers due to interface roughness scattering,” Opt. Express 13, 7779–7793 (2005). [CrossRef] [PubMed]

10.

M.-C. Phan-Huy, J.-M. Moison, J. A. Levenson, S. Richard, G. Melin, M. Douay, and Y. Quiquempois, “Surface roughness and light scattering in a small effective area microstructured fiber,” J. Lightwave. Technol. 27, 1597–1604 (2009). [CrossRef]

11.

V. Dangui, M. J. F. Digonnet, and G. S. Kino, “Modeling of the propagation loss and backscattering in air-core photonic-bandgap fibers,” J. Lightwave Technol. 17, 3783–3789 (2009). [CrossRef]

12.

E. G. Rawson, “Theory Of scattering by finite dielectric needles illuminated parallel to their axes,” J. Opt. Soc. Am. 62, 1284–1286 (1972). [CrossRef]

13.

E. G. Rawson, “Analysis of scattering from fiber waveguides with irregular core surfaces,” Appl. Opt. 13, 2370–2377 (1974). [CrossRef] [PubMed]

14.

P. Mazumder, S. L. Logunov, and S. Raghavan, “Analysis of excess scattering in optical fibers,” J. Appl. Phys. 96, 4042–4049 (2004). [CrossRef]

15.

S. G. Johnson, M. L. Povinelli, M. Soljacic, A. Karalis, S. Jacobs, and J. D. Joannopoulos, “Roughness losses and volume-current methods in photonic-crystal waveguides,” Appl. Phys. B 81, 283–293 (2005). [CrossRef]

16.

J. D. Jackson, Classical Electodynamics, 3rd ed. (John Wiley and Sons, 1998)

17.

A. W. Snyder, “Radiation losses due to variations of radius on dielectric or optical fibers,” IEEE Trans. Microwave Theory Tech. MT18(9) 608–615 (1970). [CrossRef]

18.

D. Marcuse, “Radiation losses of HE11 mode of a fiber with sinusoidally perturbed core boundary,” Appl. Opt. 14, 3021–3025 (1975). [CrossRef] [PubMed]

19.

F. Poletti, N. G. R. Broderick, D. J. Richardson, and T. M. Monro, “The effect of core asymmetries on the polarization properties of hollow core photonic bandgap fibers,” Opt. Express 139115–9124 (2005). [CrossRef] [PubMed]

20.

J. Jäckle and K. Kawasaki, “Intrinsic roughness of glass surfaces,” J. Phys. Condens. Matter 7, 4351–4358 (1995). [CrossRef]

21.

T. Sarlat, A. Lelarge, E Søndergård, and D. Vandembroucq, “Frozen capillary waves on glass surfaces: an AFM study,” Eur. Phys. J. B 54, 121–126 (2006). [CrossRef]

22.

N. V. Wheeler, M. N. Petrovich, R. Slavík, N. Baddela, E. Numkam Fokoua, J. R. Hayes, D. R. Gray, F. Poletti, and D. J. Richardson “Wide-bandwidth, low-loss, 19-cell hollow core photonic band gap fiber and its potential for low latency data transmission,” in Proceedings of Optical Fiber Communication Conference (2012), paper PDP5A.2.

23.

J. D. Joannopoulos, S. G. Johnson, J. N. Winn, and R. D. Meade, Photonic Crystal: Molding the Flow of Light, 2nd ed. (Princeton University Press, 2008).

24.

F. Benabid and P. J. Roberts, “Linear and nonlinear optical properties of hollow core photonic crystal fiber,” J. Mod. Opt. 2, 87–124 (2011). [CrossRef]

25.

R. Amezcua-Correa, N. G. R. Broderick, M. N. Petrovich, F. Poletti, and D. J. Richardson, “Design of 7 and 19 cells core air-guiding photonic crystal fibers for low-loss, wide bandwidth and dispersion controlled operation,” Opt. Express 15, 17577–17586 (2007). [CrossRef] [PubMed]

26.

J. M. Elson and J. M. Bennett, “Relation between the angular dependence of scattering and the statistical properties of optical surfaces,” J. Opt. Soc. Am. 6931–47 (1979). [CrossRef]

27.

D. Marcuse, “Mode conversion caused by surface imperfection of a dielectric slab waveguide,” Bell Syst. Tech. J. 48, 3187–3215 (1969).

OCIS Codes
(060.2310) Fiber optics and optical communications : Fiber optics
(060.2400) Fiber optics and optical communications : Fiber properties
(290.5880) Scattering : Scattering, rough surfaces

ToC Category:
Fiber Optics and Optical Communications

History
Original Manuscript: June 15, 2012
Revised Manuscript: August 13, 2012
Manuscript Accepted: August 14, 2012
Published: August 29, 2012

Citation
Eric Numkam Fokoua, Francesco Poletti, and David J. Richardson, "Analysis of light scattering from surface roughness in hollow-core photonic bandgap fibers," Opt. Express 20, 20980-20991 (2012)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-20-19-20980


Sort:  Author  |  Year  |  Journal  |  Reset  

References

  1. P. St. J. Russell, “Photonic crystal fibers,” Science299, 358–362 (2003). [CrossRef] [PubMed]
  2. R. F. Cregan, B. J. Mangan, J. C. Knight, T. A. Birks, P. S. Russell, P. J. Roberts, and D. C. Allan, “Single-mode photonic band gap guidance of light in air,” Nature (London)285, 1537–1539 (1999).
  3. C. M. Smith, N. Venkataraman, M. T. Gallagher, D. Muller, J. A. West, N. F. Borrelli, D. C. Allan, and K. W. Koch, “Low-loss hollow-core silica/air photonic bandgap fibre,” Nature (London)424, 657–659 (2003). [CrossRef]
  4. M. N. Petrovich, F. Poletti, A. van Brakel, and D. J. Richardson, “Robustly single mode hollow core photonic bandgap fiber,” Opt. Express16, 4337–4346 (2008). [CrossRef] [PubMed]
  5. B. J. Mangan, L. Farr, A. Langford, P. J. Roberts, D. P. Williams, F. Couny, M. Lawman, M. Mason, S. Coupland, R. Flea, H. Sabert, T. A. Birks, J. C. Knight, and P. St. J. Russell, “Low loss (1.7 dB/km) hollow core photonic bandgap fiber,” in Proceedings of Optical Fiber Communication Conference (2004), paper PDP24.
  6. P. J. Roberts, F. Couny, H. Sabert, B. J. Mangan, D. P. Williams, L. Farr, M. W. Mason, A. Tomlinson, T. A. Birks, J. C. Knight, and P. S. J. Russell, “Ultimate low loss of hollow-core photonic crystal fibres,” Opt. Express13236–244 (2005). [CrossRef] [PubMed]
  7. D. Marcuse, Theory of Dielectric Optical Waveguides, 2nd ed. (Academic Press, 1991).
  8. A. W. Snyder and J. D. Love, Optical Waveguide Theory (Chapman and Hall, 1983).
  9. P. J. Roberts, F. Couny, H. Sabert, B. J. Mangan, T. A. Birks, J. C. Knight, and P. S. Russell, “Loss in solid-core photonic crystal fibers due to interface roughness scattering,” Opt. Express13, 7779–7793 (2005). [CrossRef] [PubMed]
  10. M.-C. Phan-Huy, J.-M. Moison, J. A. Levenson, S. Richard, G. Melin, M. Douay, and Y. Quiquempois, “Surface roughness and light scattering in a small effective area microstructured fiber,” J. Lightwave. Technol.27, 1597–1604 (2009). [CrossRef]
  11. V. Dangui, M. J. F. Digonnet, and G. S. Kino, “Modeling of the propagation loss and backscattering in air-core photonic-bandgap fibers,” J. Lightwave Technol.17, 3783–3789 (2009). [CrossRef]
  12. E. G. Rawson, “Theory Of scattering by finite dielectric needles illuminated parallel to their axes,” J. Opt. Soc. Am.62, 1284–1286 (1972). [CrossRef]
  13. E. G. Rawson, “Analysis of scattering from fiber waveguides with irregular core surfaces,” Appl. Opt.13, 2370–2377 (1974). [CrossRef] [PubMed]
  14. P. Mazumder, S. L. Logunov, and S. Raghavan, “Analysis of excess scattering in optical fibers,” J. Appl. Phys.96, 4042–4049 (2004). [CrossRef]
  15. S. G. Johnson, M. L. Povinelli, M. Soljacic, A. Karalis, S. Jacobs, and J. D. Joannopoulos, “Roughness losses and volume-current methods in photonic-crystal waveguides,” Appl. Phys. B81, 283–293 (2005). [CrossRef]
  16. J. D. Jackson, Classical Electodynamics, 3rd ed. (John Wiley and Sons, 1998)
  17. A. W. Snyder, “Radiation losses due to variations of radius on dielectric or optical fibers,” IEEE Trans. Microwave Theory Tech.MT18(9) 608–615 (1970). [CrossRef]
  18. D. Marcuse, “Radiation losses of HE11 mode of a fiber with sinusoidally perturbed core boundary,” Appl. Opt.14, 3021–3025 (1975). [CrossRef] [PubMed]
  19. F. Poletti, N. G. R. Broderick, D. J. Richardson, and T. M. Monro, “The effect of core asymmetries on the polarization properties of hollow core photonic bandgap fibers,” Opt. Express139115–9124 (2005). [CrossRef] [PubMed]
  20. J. Jäckle and K. Kawasaki, “Intrinsic roughness of glass surfaces,” J. Phys. Condens. Matter7, 4351–4358 (1995). [CrossRef]
  21. T. Sarlat, A. Lelarge, E Søndergård, and D. Vandembroucq, “Frozen capillary waves on glass surfaces: an AFM study,” Eur. Phys. J. B54, 121–126 (2006). [CrossRef]
  22. N. V. Wheeler, M. N. Petrovich, R. Slavík, N. Baddela, E. Numkam Fokoua, J. R. Hayes, D. R. Gray, F. Poletti, and D. J. Richardson “Wide-bandwidth, low-loss, 19-cell hollow core photonic band gap fiber and its potential for low latency data transmission,” in Proceedings of Optical Fiber Communication Conference (2012), paper PDP5A.2.
  23. J. D. Joannopoulos, S. G. Johnson, J. N. Winn, and R. D. Meade, Photonic Crystal: Molding the Flow of Light, 2nd ed. (Princeton University Press, 2008).
  24. F. Benabid and P. J. Roberts, “Linear and nonlinear optical properties of hollow core photonic crystal fiber,” J. Mod. Opt.2, 87–124 (2011). [CrossRef]
  25. R. Amezcua-Correa, N. G. R. Broderick, M. N. Petrovich, F. Poletti, and D. J. Richardson, “Design of 7 and 19 cells core air-guiding photonic crystal fibers for low-loss, wide bandwidth and dispersion controlled operation,” Opt. Express15, 17577–17586 (2007). [CrossRef] [PubMed]
  26. J. M. Elson and J. M. Bennett, “Relation between the angular dependence of scattering and the statistical properties of optical surfaces,” J. Opt. Soc. Am.6931–47 (1979). [CrossRef]
  27. D. Marcuse, “Mode conversion caused by surface imperfection of a dielectric slab waveguide,” Bell Syst. Tech. J.48, 3187–3215 (1969).

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