Application of continuous wavelet transform for determination of fiber birefringence

Measurements of the group-birefringence dispersion of a microstructured optical fiber using the wavelength-scanning technique are reported. Resulting interferograms are processed using the continuous wavelet transform. We discuss application of this approach for determination of birefringence of few-mode fibers and show that with careful analysis it is possible to determine birefringence dispersion of the higher-order modes in optical fibers. © 2012 Optical Society of America OCIS codes: (060.2300) Fiber measurements; (060.2420) Fibers, polarization-maintaining; (060.4005) Microstructured fibers; (100.7410) Wavelets. References and links 1. A. Ferrando, E. Silvestre, P. Andres, J. Miret, and M. Andres, “Designing the properties of dispersion-flattened photonic crystal fibers,” Opt. Express 9(13), 687–697 (2001). 2. L. Fu, B. K. Thomas, and L. Dong, “Efficient supercontinuum generations in silica suspended core fibers,” Opt. Express 16(24), 19629–19642 (2008). 3. J. R. Hayes, J. C. Flanagan, T. M. Monro, D. J. Richardson, P. Grunewald, and R. Allott, “Square core jacketed air-clad fiber,” Opt. Express 14(22), 10345–10350 (2006). 4. A. Ortigosa-Blanch, J. C. Knight, W. J. Wadsworth, J. Arriaga, B. J. Mangan, T. A. Birks, and P. St. J. Russell, “Highly birefringent photonic crystal fibers,” Opt. Lett. 25(18), 1325–1327 (2000). 5. K. Saitoh, Y. Tsuchida, M. Koshiba, and N. A. Mortensen, “Endlessly single-mode holey fibers: the influence of core design,” Opt. Express 13(26), 10833–10839 (2005). 6. W. J. Wadsworth, R. Percival, G. Bouwmans, J. C. Knight, and P. St. J. Russell, “High power air-clad photonic crystal fibre laser,” Opt. Express 11(1), 48–53 (2003). 7. T. M. Monro, S. Warren-Smith, E. P. Schartner, A. Francois, S. Heng, H. Ebendorff-Heidepriem, and V. S. Afshar, “Sensing with suspended-core optical fibers,” Opt. Fiber Technol. 16(6), 343–356 (2010). 8. Y. Ruan, E. P. Schartner, H. Ebendorff-Heidepriem, P. Hoffmann, and T. M. Monro, “Detection of quantum-dot labelled proteins using soft glass microstructured optical fibers,” Opt. Express 15(26), 17819–17826 (2007). 9. M. Grabka, B. Wajnchold, S. Pustelny, W. Gawlik, K. Skorupski, and P. Mergo, “Experimental and theoretical study of light propagation in suspended-core optical fiber,” Acta Phys. Pol. A 118, 1127–1132 (2010). 10. M. J. Steel, T. P. White, C. Martijn de Sterke, R. C. McPhedran, and L. C. Botten, “Symmetry and degeneracy in microstructured optical fibers,” Opt. Lett. 26(8), 488–490 (2001). 11. A. Kumar and A. Ghatak, Polarization of Light with Applications in Optical Fibers (SPIE Press, 2011) 12. S. R. Norman, D. N. Payne, M. J. Adams, and A. M. Smith, “Fabrication of single-mode fibers exhibiting extremely low polarization birefringence,” Electron. Lett. 15(11), 309–311 (1979). 13. T. Martynkien, M. Szpulak, G. Statkiewicz, G. Golojuch, J. Olszewski, W. Urbanczyk, J. Wojcik, P. Mergo, M. Makara, T. Nasilowski, F. Berghmans, and H. Thienpont, “Measurements of sensitivity to hydrostatic pressure and temperature in highly birefringent photonic crystal fibers,” Opt. Quantum Electron. 39(4-6), 481–489 (2007). 14. P. Hlubina, “Measuring dispersion between modes of an optical fibre using time-domain and spectral-domain low-coherence interferometry,” J. Mod. Opt. 45(9), 1767–1774 (1998). 15. P. Hlubina, D. Ciprian, and M. Kadulova, “Wide spectral range measurement of modal birefringence in polarization-maintaining fibres,” Meas. Sci. Technol. 20(2), 025301 (2009). 16. S. C. Rashleigh, “Measurement of fiber birefringence by wavelength scanning: effect of dispersion,” Opt. Lett. 8(6), 336–338 (1983). 17. P. Hlubina, T. Martynkien, and W. Urbańczyk, “Dispersion of group and phase modal birefringence in ellipticalcore fiber measured by white-light spectral interferometry,” Opt. Express 11(22), 2793–2798 (2003). #165074 $15.00 USD Received 26 Mar 2012; revised 18 Apr 2012; accepted 24 Apr 2012; published 7 Jun 2012 (C) 2012 OSA 18 June 2012 / Vol. 20, No. 13 / OPTICS EXPRESS 13878 18. A. Ortigosa-Blanch, J. C. Knight, W. J. Wadsworth, J. Arriaga, B. J. Mangan, T. A. Birks, and P. St. J. Russell, “Highly birefringent photonic crystal fibers,” Opt. Lett. 25(18), 1325–1327 (2000). 19. J. Wojcik, P. Mergo, M. Makara, K. Poturaj, L. Czyzewska, J. Klimek, and A. Walewski, “Technology of suspended core microstructured optical fibers for evanescent wave and plasmon resonance optical fiber sensor,” Proc. SPIE 6990, 6990T (2008). 20. C. S. Burrus, R. A. Gopinath, and H. Guo, Introduction to Wavelets and Wavelet Transforms: A Primer, (Prentice Hall, 1997) 21. O. Koysal, S. E. San, S. Ozder, and F. N. Ecevit “A novel approach for the determination of birefringence dispersion in nematic liquid crystals by using the continuous wavelet transform,” Meas. Sci. Technol. 14(6), 790–


Introduction
Over the last decade, microstructured optical fibers (MOFs) have been a subject of many theoretical and experimental investigations.The interest in the fibers is mainly motivated by flexibility of engineering their structures enabling control over their optical properties.Until now, tailoring of such properties as dispersion [1], nonlinearity [2], beam shape [3], and birefringence [4] was demonstrated and specific fiber structures optimized for particular applications, e.g., endlessly single-mode transmission [5], generation of ultra-high power light [6], and chemical sensing [7] were elaborated.
In this paper, we focus on a particular class of MOFs -the suspended-core fibers (SCF).Such fiber consists of three air holes/channels surrounding the core suspended on three thin struts.Typically, a diameter of the light-guiding core is of the order of a micrometer, which makes light to penetrate the channels in a form of the evanescent field [8].The channels may be filled by a foreign medium: gas or liquid.This enables studies of the interaction between light and the medium under specific experimental conditions and opens up an avenue for applications of this type of fibers in chemical and biological sensing [8].However, in order to obtain quantitative information about the filling medium, e.g., concentration of a specific compound, a detailed knowledge of the fiber properties is required.Recently, we have investigated a number of optical properties of SCFs such as mode structure, numerical amplitude, etc [9].Herein, we concentrate on another important property of the fiber, the birefringence.
Breaking rotation symmetry of a fiber structure is one of the techniques for introducing birefringence in optical fibers.However, it was demonstrated theoretically that only fibers with a rotational symmetry of the order lower than three may reveal such structural birefringence [10].Thus, a fiber like SCF should not possess any birefringence.Nonetheless, in reality, deviations from the perfect structures are always observed and result in residual birefringence of optical fibers.
In optical fibers two parameters may be used to describe birefringence: the phase and the group birefringence.The two properties are closely related to each other and it is generally possible to calculate one of them based on spectral dependence of the other [11].The group birefringence is also related to other important characteristics of a fiber -the differential group delay and polarization mode dispersion [11].
There are many techniques used to determine birefringence and/or polarization-mode dispersion of optical fibers.One may, for instance, investigate the change of a polarization state of light departing a fiber as a function of its length (the cutting-end method) [12].This approach, however, requires destruction of the fiber.A different technique is based on mechanical coupling of different polarization modes.In such an approach the state of light at the fiber end is monitored as a function of the position of a mechanical stress-point.i.e., the length of co-propagation of orthogonal polarization modes [13].Both these techniques enable to determine phase birefringence of a fiber.Also the white-light (low-coherence) techniques can be applied for measurements of the polarization-mode dispersion [14].The most common approach utilizes an interferometer which compensates the optical-path difference introduced by the fiber at certain wavelengths and enables to determine the dispersion [15].
In this paper we use another approach to analyze birefringence, -the wavelength-scanning method (WSM) [16].It combines different aspects of all the above mentioned techniques.In the method a fiber placed between two crossed polarizers is illuminated by a broadband light source.The first polarizer enables to excite orthogonal polarization modes in the fiber, while the second polarizer allows one to observe interference between the modes.Application of a broadband light makes it possible to observe the interference in a spectral domain based on relation [17] for the group birefringence: where λ is the wavelength of light at which oscillation is measured, λ ∆ is the oscillation period, and L is the fiber length.WSM offers a number of advantages to determine optical fiber birefringence.It is simple, non-destructive, and allows to simultaneously measure birefringence in wide spectral range.It seems to be particularly well suited for investigation of MOFs as it does not require application of mechanical disturbance to the fiber; it was shown in Ref [18].that propagation of light in SCFs is hardly affected by mechanical stress as the channels surrounding the core shield it from external disturbances.
Herein, in order to retrieve information about birefringence of SCF, data recorded using WSM are analyzed using the continuous wavelet transform.It allows us to determine the birefringence of the fundamental mode but also to determine the birefringence of higher-order modes, which, to our best knowledge, was done for the first time.We show that the higherorder mode-birefringence strongly depends on the wavelength, which is an obstacle in its determination using traditional approaches, e.g. the Fourier transform.A possibility of measuring birefringence of many modes is important in analysis of few-mode optical fibers, in particular SCFs.

Experimental apparatus and results
Our experimental setup is shown in Fig. 1.A halogen quartz lamp (StallarNet SL-1) serves as a broadband light source and a fiber optic spectrometer (AvaSpec-3648) is used for transmission measurements in the 500-1000-nm range with a resolution of 0.5 nm.To couple broadband light into the fiber, we use a system of several spherical mirrors and optical fibers denoted in Fig. 1 as a beam-shaping system.The investigated fiber is situated between two linear polarizers.The first linearly polarizes the beam and excites different polarization modes in the fiber (incident polarization at 45° with respect to the optical axis of the fiber), while the second is used to observe the mode interference.The axes of both polarizers are crossed, which provides the best signal-to-noise ratio, note, however, that other angles are also possible and, in some situations, are worth considering (see below).
The SCF we examine in this work (Fig. 2) was fabricated in the Laboratory of Optical Fibre Technology at the Marie Curie-Sklodowska University in Lublin, Poland [19].This silica fiber possess a small, triangular-like silica core (800 nm diameter of the inscribed circle) attached to three thin struts (length of 9 µm and a thickness of 70 nm) to the solid fiber cladding.The fiber supports propagation of a number of modes within a range from 500 nm to 1000 nm.The SCF structures which have three-fold rotation symmetry are not birefringent and all their modes either have the same symmetry or exist in degenerated doublets.Thus any imperfection of the structure which breaks the symmetry will lead to a very high birefringence between higher-order modes [10].A typical signal measured with the WSM is shown in Fig. 3.It was obtained with the spectrometer detecting the light exiting the 290-mm long SCF and interfering at the output polarizer.Normalization of the spectra eliminates influence of spectral properties of the fiber, halogen lamp, spectrometer and other optical devices, hence the recorded spectra demonstrate only the mode interference.Besides the clearly visible oscillations of about 25-nm periodicity, the spectrum shown in Fig. 3 reveals higher-frequency components of the signal, most distinctly at longer wavelengths.It is important to note that this high frequency modulation is stable in time, hence is not noise but, as discussed in Sec.3.3, it originates from group birefringence of the fiber.

Continuous wavelet transform
A successful technique that typically enables to distinguish different components of multifrequency signal is the Fourier transform.However, application of the technique in the considered case is not possible due to the strong dependence of the oscillations on the wavelength of light (see Fig. 3).The case, however, may be conveniently analyzed with the help of the Continuous Wavelet Transform (CWT) [20].Note that CWT has already proved its usefulness for analysis of liquid crystal birefringence [21], but to our best knowledge, has never been applied for investigation of optical fibers.

Continuous wavelet transform of a given function ( )
Function φ(x) is the, so-called, mother wavelet which is the basis for constructing the dilated and translated wavelets by using parameters a and b corresponding to the frequency and the central position of the wavelet, respectively.Choosing the appropriate wavelet for a given application is not trivial.In our case of fiber birefringence we apply the Morlet wavelet, which is a Gaussian function modulated by a plane wave [20]:  Figure 4(a) presents the two-dimensional wavelet spectrum of the data shown in Fig. 3 (note that the plots were adapted for the considered case, i.e., wavelet frequency was replaced by its period and position of the wavelet is labeled as a wavelength).The spectrum reveals three well visible traces marked as A, B, and C. The strongest trace corresponds to the oscillation period of about 25 nm that does not significantly vary with the wavelength.This result may be almost completely reproduced by application of the Fourier transform [Fig.4(b)].However, Fig. 4(a) reveals other important features.For example, the other trace marked as B is visible in the range of 750-900 nm but its period decreases from ~10 nm for shorter wavelengths to ~2 nm for longer wavelengths.The last trace, marked as C, is observed for the 600-750 nm range.Its period begins at ~3 nm at 550 nm and rises until ~10 nm at 650 nm where it start decreasing back to ~2 nm at 750 nm.Since none of these traces is observed using the Fourier transform, Fig. 4(a) demonstrates well the potential of the CWT method.

Data filtering
In order to identify the signal associated with the fiber from the overall background, the spectrum shown in Fig. 4(a) is analyzed by searching for local extremes.It allows one to obtain two-dimensional point plot [Fig.5(a)], which, however, due to simplicity of the applied identification technique, contains many points not associated with the fiber properties.The signal is hence processed further by application of three main criteria.The first two are set by the resolution of the spectrometer which limits the lowest and the highest periods that can be determined.Somehow arbitrary we assume that no less than four points per period are required to detect the oscillations (in fact, two points per period are enough, but applying a more restrictive criterion increases credibility of the results).Thus the applied limits are: 2 nm (four times the interferometer resolution of 0.5 nm) and ~100 nm (about one quarter of the 450 nm interferometer bandwidth).The third criterion is related to the threshold value of set to 10% of the maximum value available with the normalized spectrum.It needs to be stressed that the threshold has to be set very carefully as too high threshold causes loss of data points, while too low threshold level leaves too much noise.The raw data and the data processed using the criteria discussed above are presented in Fig. 5. Since the applied filtering technique is rather simple, some points originating from the fiber properties may be lost.Nonetheless even such rough filtering yields valuable information about fiber (see below).In particular, the points depicted in Fig. 5 group in a way which reveals dispersion of the resolved modes.Such analysis is was not possible with the regular standard FFT method (compare with Fig. 4(b)).

Group birefringence
The data shown in Fig. 5(b) may be recalculated into the group birefringence based on Eq. (1).In order to verify that they are associated with the fiber properties and are not caused by any apparatus artifacts, additional cross-checks are performed.The most basic test is to measure pieces of the fiber of different lengths.If certain points are related to the fiber properties, they should lay along a given curve after length normalization.At the same time the points which are not associated with the birefringence are just randomly scattered over the whole plot area and represent a uniform background.The purpose of this procedure is to show that obtained points are in fact originated by the fiber birefringence.In case of typical birefringence measurement, such step is not required and analysis of single interferogram is sufficient.
Other test consisted in modification of the measurement geometry, i.e. rotation of the optical axes of the polarizers.Such change deteriorated the contrast of the observed signal, but also allowed us to observe signals for different light coupling conditions.This test was very important, since it revealed birefringence of the higher-order modes where proper polarizer arrangement was not easy to accomplish.Fig. 6.Group birefringence values collected from measurements done for different arrangement of polarizers and different fiber lengths.Dashed lines limit the birefringence range accessible in measurements using our spectrometer.
The data taken with different fiber lengths and polarization arrangements are brought together in Fig. 6.It is clearly visible that most of the collected points form four strongly wavelength-dependent traces.Points originating from the birefringence of the fundamental modes lay on a curve which rises monotonously with λ from about respectively.For longer wavelength their birefringence decreases.The behavior of the curves around 600-650 nm and 900-950 nm may suggest that the group birefringence changes its sign in these regions and for certain wavelengths is equal zero.However, since in such measurements only the absolute value of the group birefringence is measured, the sign of birefringence is not known.The fourth trace starts at about 725 nm and ranges as far as near 950 nm.
It is worth noting that the resolution of a spectrometer limits the range of the birefringence that can be measured.High birefringence yields oscillations with a short period [see Eq. ( 1)] and requires a spectrometer with high resolution.The range of birefringence that can be measured using a spectrometer of a given resolution changes with the wavelength of light and the length of a fiber.The limits imposed by the spectrometer in presented measurements are marked in Fig. 6 by two dashed lines.

Conclusions
We described the approach for determination of the birefringence dispersion of few-mode optical fibers.It is based on the wavelength scanning technique and continuous wavelet transform.Application of these techniques allowed us to measure group-birefringence dispersion of the suspended-core fiber in a wide spectral range.This was done despite of a very strong wavelength dependence of the birefringence signals.This approach can be improved by more careful data filtering for better compromise between the number of eliminated data points and noise.The range of birefringence values that can be measured is limited by the range and accuracy of the spectrometer.
To the best to our knowledge, group birefringence of the higher-order modes has been measured for the first time.We expect the described method will prove useful for determination of polarization modes associated with the measured signals and for verification of various models of the suspended-core fibers that we planned to used for chemical sensing.

Fig. 1 .
Fig. 1.Experimental setup for birefringence measurements using the wavelength scanning method.

Fig. 3 .
Fig. 3. Typical transmission spectrum obtained for a 290 mm long piece of a suspended-core fiber.

.
Calculating CWT of a signal allows one to generate a twodimensional spectrum, where each point refers to a certain pair (a, b).

Fig. 4 .
Fig. 4. (a) Real part of the wavelet spectrum of the signal shown in Fig. 3. Three distinct visible traces are marked as A, B, C, (b) FFT spectrum of the same signal.Besides of a large zerofrequency component, only the primary birefringence mode is visible (at frequency close to 0.04 nm −1 ).Other features relevant for CWT trances B and C should be located in the 0.1-0.5 nm −1 frequency range, but are not visible.

Fig. 5 .
Fig. 5. (a) Local extremes of the real part of the wavelet spectrum without filtering, (b) Signal filtered out using the criteria described in the text.Each point of the plots is related to an extremum in the CWT spectrum in Fig. 4(a) and defines a pair (λ, ∆λ).The pairs group into bent lines which reflect the systematic wavelength dependence of the group birefringence of the resolved modes.

−
at 900 nm.This is a very high group birefringence for the fiber which was fabricated without intentionally introduced birefringence.Other three traces are connected with higher-order modes.Two of them are rising with λ until ~800 nm where they reach their maximum values, x is expressed as: #165074 -$15.00USD Received 26 Mar 2012; revised 18 Apr 2012; accepted 24 Apr 2012; published 7 Jun 2012