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 Optical Society of America
In their most common form, hollow-core photonic bandgap fibers (HC-PBGFs) consist of a two-dimensional triangular lattice of air holes embedded in a host material with the core formed by removing a number of adjacent capillaries, usually 3, 7 or 19 [1–4]. Light guidance in the air core is therefore achieved by virtue of the out-of plane photonic bandgap of the periodic cladding. When designed properly, HC-PBGFs can support fundamental modes with more than 99% of the guided power propagating in air at wavelengths within the bandgap. As a consequence, Rayleigh scattering and phonon or impurity absorption that dominate the attenuation in more conventional step or graded index fibres have a lower contribution to the overall loss in HC-PBGFs. This raises the possibility of producing HC-PBGFs with very low attenuation values which, together with the nearly three orders of magnitude decrease in effective nonlinearity and the much reduced latency make HC-PBGFs an interesting candidate for next generation long-haul data transmission systems.
However, the lowest loss value reported to date for HC-PBGFs stands at 1.7dB/km , 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 . 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, 6].
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, 8]. Roberts et al.  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  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.  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, 13] and later used by Mazumder et al.  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.
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.
Any local departure from the ideal boundary can be treated as a dipole of volume dV = |f (s, z)|dsdz which is excited by the incident modal field. This approximation is justified by the fact that the rms roughness resulting from frozen-in SCWs is of the order of 0.1nm, much smaller than the wavelength of light. Since the normal component of the electric displacement D and the parallel component of the electric field E are continuous at all interfaces, the induced dipole moment is generally given as:15]. In general, these quantities will depend on the shape of the dipole and the refractive index contrast. Johnson et al.  have computed the polarizabilities α and γ for different dipole shapes and index contrasts. From their results, the effects of neglecting the exact shape of the dipole results in an error of less than 1% for silica-air interfaces at infrared wavelengths (index contrasts of ∼ 0.5 and ∼ 2). Further neglecting the effects of the discontinuity of the normal component of the electric field results in an additional error which is well below 5%. The polarizability additionally depends on the sign of the roughness defect. For example, when the perturbation is such that glass protrudes into air we speak of positive roughness, otherwise, the roughness is said to be negative. We take the polarizability α+ of a positive roughness “bump” as that of a glass sphere suspended in air and that of a negative one α− as the polarizability of an air sphere in a glass background. The two differ in both sign and magnitude [15, 16], but since f is random with zero average, we will neglect all the effects above and write the induced dipole moment in terms of the incident electric field at the interface as:
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 :Fig. 2 as:
In this expression, βz represents the phase build up of the incident mode field between z = 0 and z, while k0zcosϕ − k0y′sinϕ 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 isEq. (3) and (2) becomes:
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. . 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:Eq. (6) becomes: 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, 8]. 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:26] Eq. (8) into Eq. (9), we obtain with the help of Eq. (10):
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, 18]. 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 isFig. 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 . This not only validates our simple approach but also shows that it can be implemented in combination with any arbitrary mode-solving numerical tool.
4. Roughness scattering in HC-PBGFs
Since no analytical solution for wave propagation in HC-PBGFs can be obtained, our general approach is to solve for the field components at the interfaces of the guided modes using a fully vectorial finite element method . A realistic model of the surface roughness inside the fiber is then incorporated and the far-field distribution of the scattered light obtained from Eq. (8). As no assumption has been made on the nature of the incident modal field, Eq. (8) can also be used to analyze the scattering loss of surface modes.
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 :
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 :
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 . 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,10]. 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 , 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 . 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.
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 . It has been shown experimentally that when a given fiber cross-section is drawn to different diameters, its minimum loss decreases as 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,24]. Despite the explicit in Eq. (11), our model confirms the experiments by predicting a scattering loss that indeed follows a 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.
On the same figure, we also show how the normalized interface field intensity6, 25]. Here, F was calibrated using measurements and simulations for the fibre in  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
In order to accurately predict roughness scattering loss in HC-PBGFs, accurate knowledge of the roughness statistics and its PSD in particular is of utmost importance. However, the practical limitations of AFM measurements make it impossible to obtain information about the very important low spatial frequency region of the spectrum and alternative techniques are yet to be implemented. In the previous sections, we have assumed that the PSD for surface capillary waves is cut-off at an arbitrary spatial frequency κc = 0.1μm−1. In Fig. 6(a), we plot the computed scattering loss for our HC-PBGF as a function of the cut-off frequency at the wavelength of λc = 1.55μm. As κc increases, the PSD S(κc) is lower for a broad range of frequencies, leading to the observed decrease in the scattering loss.
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 . 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 .
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 respectively, with Lc being the correlation length and σ the corresponding rms roughness. The PSD arising from these auto-correlation functions are a Lorentzian and a Gaussian 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 . 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.
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.
This work was supported by the EU 7th Framework Programme under grant agreement 228033 ( MODE-GAP).
References and links
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]
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]
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]
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).