The spectral characteristics of a fiber Bragg grating (FBG) with a transversely inhomogeneous refractive index profile, differs considerably from that of a transversely uniform one. Transmission spectra of inhomogeneous and asymmetric FBGs that have been inscribed with focused ultrashort pulses with the so-called point-by-point technique are investigated. The cladding mode resonances of such FBGs can span a full octave in the spectrum and are very pronounced (deeper than 20dB). Using a coupled-mode approach, we compute the strength of resonant coupling and find that coupling into cladding modes of higher azimuthal order is very sensitive to the position of the modification in the core. Exploiting these properties allows precise control of such reflections and may lead to many new sensing applications.
© 2011 Optical Society of America
In a typical transmission spectrum of fiber Bragg gratings (FBGs), it is common to observe a comb of resonances on the shorter wavelength side of the Bragg peak. These dips result from the resonant coupling between the forward-propagating core mode and modes that propagate backwards within the cladding of the fiber [1, 2]. Since modes guided in the cladding of a fiber are very sensitive to bending of the fiber as well as the surrounding environment, they are widely exploited for sensing applications. In such sensor arrangements, FBGs and long period gratings (LPGs) are used as core to cladding mode converters . However, conventional FBGs, in which the cross-section of the core is uniformly modified, couple light only into cladding modes of the lowest azimuthal order . In general, coupling into higher azimuthal fiber modes can be realized by FBGs that are tilted [4–6] or have an asymmetric transverse index profile . More elaborate in-fiber devices to achieve strong cladding mode resonances involve tapered fiber sections  or deliberately misaligned fiber splices .
In this paper, we investigate how strong cladding mode coupling can be achieved and controlled with highly localized FBGs. Such FBGs can be fabricated using focused ultrashort lasers operating in the near infrared (NIR) [10, 11]. Due to the nonlinear nature of the absorption process of light pulses in the femtosecond regime, highly localized refractive index profiles can be obtained in non-photosensitive materials without damaging the surrounding fiber material. The induced refractive index change is no longer confined to the photosensitive region of the fiber (as is the case for FBGs written with ultra-violet lasers), but is determined solely by the focusing geometry of the ultrashort laser. Both the size and position of the modifications depend heavily on the alignment and focusing of the writing beam [12, 13]. In this paper we demonstrate fs-pulse written, untilted FBG with very pronounced (up to 25dB) cladding mode coupling. Our aim is to investigate and understand these cladding mode spectra in detail in terms of a full vectorial coupled mode theory. Within the three-layer model, a consistent mode classification is proposed that helps to distinguish cladding-modes with high core-mode overlap. Our work shows that excitation of cladding-modes with strong selectivity and precision is possible, detailing how the coupling behavior is determined by the fiber geometry and to which extend it can be manipulated with the cross-section of the FBG.
2. Realization and spectra of highly localized FBG
2.1. Fabrication and characterization
The details of the grating inscription procedure closely follow that, which has been reported in [12, 14]. A femtosecond laser (Spectra-Physics, Hurricane) that delivered pulses shorter than 110 fs at a wavelength of 800 nm was utilized for grating inscription. A 20× oil immersion objective (NA = 0.8) was used for focusing the ultrashort pulses into the core of the fiber, as well as for imaging the fiber core before and after inscription of the gratings. All experiments were performed with a standard single mode fiber (Corning SMF-28e). The fiber was completely stripped of its polymer cladding prior to grating inscription and subsequently threaded through a glass ferrule which was placed in close proximity to the objective lens. The lens, ferrule and fiber were immersed in an index matching oil in order to minimize refraction as the pulse propagates towards the focal spot. The ferrule was positioned such that the core of the fiber was located at the focus of the laser beam and the fiber was then pulled through the ferrule at a constant velocity while the grating was inscribed. By adjusting the ferrule’s location with respect to the focus of the lens, the transverse position of the modifications of the FBG could be controlled to micrometer accuracy within the fiber core. The error in the longitudinal position of the modifications was much smaller, because the fiber was pulled by an air-bearing translation stage designed for constant velocity traversals, and the repetition rate of the laser was constant. All gratings were inscribed in SMF-28e fibers and are 20 mm long. Except where stated otherwise, the gratings have a period Λ = 1.062 μm that reflects 1.5 μm light in the second order.
Pulse energies between 200 and 275 nJ were used, which provided enough irradiance within the focal volume to cause a micro-explosion in the fiber, leaving behind a micro-void and hence forming a type II-IR Point-by-Point (PbP) FBG . The micro-voids are surrounded by a compressed region of increased refractive index . In previous work, we have demonstrated that the net effective refractive index change is negative, although this is not significant in the present work . Microscope images of the grating modifications taken from the top (a) and side (b) views are shown in Fig. 1. The micro-void regions, which are highlighted in Fig. 1, measure 0.4 μm in width by 1.9 μm in height. They are enclosed by a densified Type I shell of width approximately 1 μm and height 8 μm, where the index is slightly increased. The modifications are not necessarily centered, thus having an offsets of dx and dy.
2.2. Spectral response
The measurement of the transmission and reflection properties was carried out with a swept wavelength system (SWS) (JDS Uniphase), that consisted of a narrow line-width external cavity laser diode. The probing wavelength of the laser diode can be continuously swept from λ = 1520 to 1570 nm with a resolution of 3 pm. A typical transmission spectrum is shown in Fig. 2. In contrast to FBGs written with conventional techniques, the cladding mode structure is very pronounced: the reflectivity of core light into cladding modes is higher than 70 percent, even at wavelengths 20 nm shorter than the Bragg wavelength λB = 1540 nm. For gratings at 1555 nm, we see strong resonances 35 nm below the Bragg peak, limited only by the wavelength range of our measurement system. Comparable cladding mode spectra can be obtained with conventional techniques only with blazed FBGs with large tilt angles, e.g. 80° , whereas our gratings are untilted.
As detailed in section 3, each resonance can be attributed to the coupling of the core mode to one or more definite cladding modes characterized by their radial and azimuthal mode numbers. It can be seen that the spectrum contains three sets of resonances with distinct well-defined envelopes. In section 4.1 we show that the envelopes correspond to modes with even (green dots) or odd (red dots) azimuthal mode numbers.
Close inspection reveals that many of the resonances, especially at wavelengths lower than 1530 nm (see lower left panel of Fig. 2), actually consist of doublets (marked by red and blue dots respectively). No such splitting is observed at the low order (long wavelength) resonances. These doublets are a direct consequence of the vectorial nature of the cladding guided modes, and are the three-layer analog to the standard EH/HE modes of the fiber which will be discussed in detail in section 3.
Finally, the lower right panel of Fig. 2 provides a striking illustration of the strength of the cladding mode coupling in these gratings. It shows the response at wavelengths longer than the second-order Bragg peak at 1540 nm. A sequence of regularly-spaced pairs of resonances 1–2 dB deep, each themselves with a recurring splitting pattern, is clearly visible. We interpret these peaks as the tail of the first-order cladding mode spectrum, and they therefore correspond to resonances many hundreds of modes below the fundamental. The actual first order Bragg peak would occur in the mid-IR around 3 μm and is beyond our observable range due to losses and the spectral range of our near-IR characterization setup.
In the following sections, we perform calculations to understand these types of cladding mode spectra and explain their features in precise detail.
3. Classification and properties of modes in three layer fibers
For the purpose of calculating and classifying the cladding modes we consider a standard, cylindrical step-index fiber supporting a single core-guided mode over the wavelength range of interest. The elliptical modifications introduce a small anisotropy to the core and the modes are technically birefringent, but this effect is insignificant for the cladding modes. We therefore treat the fiber as having perfect circular symmetry, with core radius a1 and cladding radius a2. The refractive indices are n1, n2, respectively.
Exploiting the small difference in n1 and n2, cladding modes are most often calculated with a two-layer model consisting of a solid rod of index n2 and radius a2 surrounded by an air cladding with index 1. They are further often treated in the approximation of weak-guidance or linearly polarized (LP) modes . This approximation is usually sufficient for computing the position of the resonances in the spectrum [17, 18], since the propagation constants of the true HE and EH hybrid modes of the step index fiber differ only by a very small amount Δβ .
It has previously been found, however, that in order to accurately calculate coupling coefficients between core and cladding modes, both these approximations must be avoided . The core-cladding boundary influences the electric field distribution of the lower order modes especially. Hence it is necessary to solve the cladding modes for the full three-layer structure . Moreover, because fiber modes of higher azimuthal order are not linearly polarized, we need a full vectorial approach  to accurately explain splittings in the cladding-mode spectra that do not appear in the LP picture.
In the remainder of the paper, the core radius a1 = 4.15 μm and the cladding radius a2 = 62.5 μm. The refractive indices of the core, cladding and surrounding air are n1 = 1.4670, n2 = 1.4618 and n3 = 1.0 respectively. Modes travel in the z direction with the propagation constant β = 2πn̄/λ, where n̄ is the effective index. A mode is guided within the core if n̄ > n2, and in the cladding if n2 > n̄ > n3. To excite the system, we launch a single linearly polarized HE1,1 core mode.
3.1. Field expressions for the modes
In our analysis, we make extensive use of the three-layer fiber model in the form laid out by Erdogan . He adapted the full analytical solution of the three-layer treatment of Tsao , which yields all hybrid modes with the azimuthal and radial integer indices l and m. However, in , Erdogan focused solely on cladding mode reflections of conventional fiber gratings, for which only l = 1 resonances occur. In our work, higher order modes with l > 1 are needed. We use the notation employed by Erdogan, but express the azimuthal dependence in trigonometric rather than exponential form. In cylindrical coordinates (r, ϕ, z), the electric E and magnetic H fields of the cladding modes inside the core (r < a1) can be expressed in terms of Bessel functions Jn of the first kind
Note that Eqs. (1) represent two orthogonal sets of solutions, distinguished by the rotation angle φ, which we take as φ = 0 or φ = −π/2. Thus, all hybrid mode solutions of the form in Eqs. (1) appear as degenerate pairs of fields with orthogonal polarization states. Since in our experiments the input mode HE1,1 is y-polarized, it is represented purely by the φ = 0 state. In general, a grating may produce coupling to cladding modes of both φ = 0 and φ = −π/2. As discussed below however, careful alignment of the polarization as well as the defect orientation and location eliminates coupling to the perpendicular φ = −π/2 modes.
3.2. Classification of the modes
It is useful to recall some standard properties of the core-guided hybrid modes of the two-layer model. In that system, the modes are classified as HE or EH according to the value of Plm. The HE modes have −1 < Plm < 0 and the EH modes have Plm > 1, with Plm ≈ ±1 for the strongly bound modes, (and Plm = ±1 in the LP limit). For azimuthal number l, the transverse field components of the HE modes have predominantly a Jl–1(u1r) radial dependence, whereas the EH modes are dominated by the Jl+1(u1r) Bessel function. Consequently, a mode HEl,m has a time-averaged Poynting flux 〈Sz〉 with peaks closer to the origin than the corresponding EHl,m mode. For fixed l, the modes occur in a strict alternating sequence of HE and EH, (excepting the fundamental HE1,1).
The mode parameter Plm also turns out to be very convenient for the classification of cladding modes in the three-layer fiber. However, the HE/EH criterion has to be modified. For further discussion, we define the projected normalized wavevector
In Fig. 3, we plot the hybrid mode parameter Plm against the normalized wavevector for the lowest six azimuthal orders. The plots show that Plm has a complex dependence and is not restricted to the ranges −1 < P < 0 and P > 1. Adjusting the two-layer fiber definitions, we designate the modes with |Plm| < 1 as HE (colored blue) and those with |Plm| > 1 as EH (colored green). Note that with this definition, we retain consistency with the classification of core modes in the two-layer picture, and the modes again appear in a strict sequence alternating between HE and EH. Now in Fig. 3 we find that for each value of l, the diagrams can be divided into regions according to the sign of Plm for either set of modes. In fact, the transitions, indicated by vertical dashed lines, correspond to mode indices m at which an additional intensity ring enters the core. We refer to these discontinuities as “virtual cut-offs” Ũlm′, where the virtual radial number m′ is to be associated with the number of maxima within the core. The cutoff value Ũlm′ is the value of Ulm at which the m′ th maximum appears in the core. The intensity distributions near and in the fiber core for a number of modes are plotted as insets. It can indeed be seen that a further intensity ring appears within the core each time there is a discontinuity in the Plm plot.
The character of the modes changes qualitatively at the first few virtual cutoffs. Cladding modes with Ulm < Ũl,1′ have essentially no intensity within the core. (As with the classic two layer model, Ũ1,1′ ≡ 0. Therefore, l = 1 modes always guide some light within the core region.) When the normalized wave-vector satisfies Ũl,1′ < U < Ũl,2′, HE modes (Plm ≈ −1) have one maximum within the core. In contrast, EH modes (Plm ≈ 1) have a minimum at the center and have significant field contributions only near the core-cladding boundary. Thus the spatial distribution of energy is distinct for HE versus EH. The polarization of each mode in this regime is quite complex and far from uniform across the fiber.
The propagation of modes with U > Ũl,2 is more determined by total internal reflection at the cladding/surrounding interface than the core/cladding interface. In this region we have either |Plm| ≈ 0 (HE), or |Plm| ≫ 1 (EH), which is indicative of modes for which the energy in the longitudinal components resides almost purely in Hz (HE) or purely in Ez (EH). As one or other longitudinal component is close to zero, consequently the transverse field of the HE modes eventually become close to purely azimuthally polarized (quasi-TE) and that of the EH modes becomes strongly radially polarized (quasi-TM) .
4. Calculation of the spectra
We now turn to the detailed explanation of the form of the cladding mode spectra, calculating the resonant wavelengths, the coupling coefficients and finally the complete transmission spectra.
4.1. Resonance wavelengths
One of the challenges of computing FBG spectra of strong cladding mode coupling is that it usually occurs over large bandwidth ranging from 20 to several hundred nanometers. All resonances that result from the reflection of the incoming mode HE1,1 into cladding modes (l,m) by the grating of period Λ are defined by the phase matching condition
Since Eq. (4) involves the wavelength implicitly through the effective indices of the modes, computing the resonance wavelengths requires solving the dispersion relation for hundreds of modes at finely discretized wavelengths. Within a bandwidth of 30 nm, the material dispersion of the fiber can usually be neglected as it is of the order of 10−4 and affects both core and cladding refractive indices in a similar fashion. Therefore the index contrast n1 – n2 varies only weakly. However, the waveguide dispersion is significant and must be accounted for. Solving Eq. (4) exactly is relatively straightforward, but considerable care is required in solving the dispersion relation (11) to find all the required effective indices n̄lm at all wavelengths, since some of the EH modes appear as resonances as narrow as one part in a million.
If the intention is to identify particular modes in a measured spectrum rather than precisely predict their location in advance, a more convenient approximate approach that provides additional insight is to solve the dispersion relation at the Bragg wavelength λ0 and extrapolate the waveguide dispersion of the propagation constants using the relation 
In order to precisely predict the resonance wavelengths using Eq. (5), for each mode the dispersion free constantEq. (5) and combining that equation with Eq. (4), we obtain a quadratic equation for the wavelength of the (l,m) resonance Figure 4 depicts a similar transmission spectrum to that in Fig. 2, together with the three-layer multimode resonances, computed according to the steps just described. The modes were calculated at λ0 = 1555 nm. As expected from the LP approximation in Eq. (5), the numerical solution also yields two sets of nearly-degenerate resonance doublets, one for odd l (red lines) and one for even l (blue lines). For given l, each doublet corresponds to an HE/EH pair. The alignment of the measured and calculated resonances is good. For low wavelengths (Fig. 4(b)), there is a small wavelength shift compared to the measured spectra which arises because of the neglected material dispersion, however the splitting of the HE/EH doublets is in very close agreement. In this range, the small difference in the propagation constant βHE – βEH of the HE/EH doublets is constant, as expected from the standard 2-layer hybrid mode solution applied to the cladding-air boundary . This is not the case for cladding resonances closer to the Bragg peak (Fig. 4(c)): βHE – βEH becomes very small in the vicinity of a virtual cut-off. It can also be seen that in this regime, the doublets of a different azimuthal order are shifted such that the EH resonance of a lower l is in most but not all cases degenerate to the HE resonance for l + 2. The use of LP modes LPℓ,m = HEℓ+1,m + EHℓ–1,m still leads to surprisingly good results for coupling strength of lower m cladding modes , since the contribution of higher l modes cancels out in transversally homogeneous FBGs as shown in the following sections.
4.2. Computation of coupling constants and spectral response
The interaction of the modes with the FBG can be computed by solving the coupled mode equations as established by Kogelnik . The unperturbed mode fields of the fiber and their propagation constants have been calculated beforehand as described in the preceding sections. Here, we use the coupled mode equations adapted for FBGs by Erdogan . This set of differential equations is determined by the coupling constants, which are obtained by evaluating an overlap integral of the interacting modes with the perturbation. In the case of FBGs, the transverse coupling constant that governs the coupling of the incoming fundamental mode HE1,1 to any other reflected mode (l,m) is computed with1, 4].
4.3. Fundamental coupling properties of the modes
Before characterizing the modification term Δɛ in Eq. (8), it is useful to elaborate on the kernel of the coupling constant integral, since this can be evaluated before integrating with the help of the results of section 3. Here, the incoming mode is always the HE1,1 core mode with φ = 0. This mode couples with any other other hybrid mode, characterized by Plm and with
If the FBG extends homogeneously over the whole cross section of the fiber core and its surrounding, Eq. (8) would be equivalent to the orthogonality relation of the modes, evaluating to zero for any mode (l,m) other than the core mode HE1,1. In general, coupling between modes of different indices occurs if these conditions are violated. In conventional FBGs for example, the grating is only formed in the photosensitive core and not the cladding. This case corresponds to a trivial integration over ϕ and any terms with a sine or cosine cancel out. Thus, only the first term of Eq. (9) remains for l = 1 and arbitrary m, while the perpendicular kernel K⊥ always sums to zero. This case is described in detail in . In particular, coupling to higher order modes with l > 1 is only possible if the refractive index of the FBG is inhomogeneous or asymmetric. The strength of such coupling is influenced by the cross section of the FBG through its non-uniformity and localization.
The classification of the fiber modes as discussed in section 3.2 is reflected in the kernel, too. For the cladding modes below the Ũl,2′ cutoff, the strength of the coupling coefficient has a pronounced alternating character as a function of m due to the hybrid mode parameter Plm ≈ ±1: The (1 – P1m) term in Eq. (9) is responsible for the strong selective character. While for HE modes the Jl–1(ulmr) term dominates, it is the Jl+1(ulmr) term for EH. In contrast, cladding modes beyond the Ũl,2′ cut-off have quasi-TM/TE core fields. Because of this, Plm ≠ ±1 and HE and EH modes differ less in their coupling behavior.
4.4. Coupling coefficients and spectra
We now proceed to calculate the coupling coefficients and complete transmission spectra for highly-localized FBGs. For such gratings the position of the point-modification usually has a bigger impact on the coupling behavior than its shape and refractive index profile. Important exceptions are discussed at the end of section 5. In our experiments, the FBG consisted of femtosecond pulse induced micro voids which are approximately prolate spheroids (see Fig. 1 and Ref. [11, 14]). They are surrounded by a densified shell of positive refractive index . Quantitative measurements of its refractive index suffer from a high uncertainty: values in the literature range from 1 · 10−3  up to 5 · 10−3 in one case . The impact of the densified shell on the coupling constants has been calculated for the fundamental mode, using finite element analysis: For an increasing positive refractive index change, the coupling strength monotonically rises , but the overall impact is small compared to the effect of the general uncertainty in the dimensions of the micro-void, and so we do not explicitly include the densified region in our model.
In the following, the transversal refractive index profile of the micro-void is approximated by a homogeneous ellipse (with semi- major axis 0.95 μm and semi-minor axis 0.2 μm). Any further asymmetries of the void and its densified shell are therefore neglected, accepting that the resulting coupling constants could be systematically too small. To reduce the geometry to the form of Eq. (8), we treat the grating as a periodic square modulation with an elliptical cross-section and an index contrast obtained by averaging the local index of the true structure over a whole period. Under these assumptions, the perturbation of the dielectric constant can be approximated with14].
In our fabricated gratings we can position the voids arbitrarily in the core (although there is some variation in the void displacement along the core due to vibration in the ferrule during the writing process). The strength of the cladding mode coupling strongly depends on small displacements of the micro voids from the central axis of the fiber (results presented below in Fig. 6). For inhomogeneous FBGs, the launch polarization has a significant influence on the overall coupling strength because of coupling to both sets of modes (Eq. (9)). These polarization dependent effects will be discussed in more detail in a forthcoming publication. In the following calculations, we emphasize more the dependence of the coupling to the radial position of the highly localized FBG. We therefore choose the special case for which the launch polarization is aligned with the long axis of the modifications and the modifications are only displaced along the x-axis. The perpendicular kernels K⊥ therefore cancel out and can be neglected.
The coupling integral was evaluated for displacements dx in the x-direction, ranging from 0 to the core cladding boundary at 4.15 μm. The resulting coupling constants for the parallel kernel K|| are plotted in Fig. 5(a) as a function of radial mode order m and displacement dx for azimuthal orders l =1–6. The horizontal white dotted lines show the locations of the virtual cutoffs for each value of l. Since resonances above Ũ7,1′ lie beyond the wavelength range of our measured spectra, we expect negligible coupling to modes of l > 6, so that the range l =1–6 is appropriate. All of the concepts of the preceding section are present in these plots: there is almost no coupling below the Ũl,1′ cutoff, pure HE coupling for modes below Ũl,2′ (as indicated by the bright red bands for odd m (HE) and dark bands for the even m (EH) modes); and reduced contrast between the HE and EH selective coupling above Ũl,2′, when the nature of the modes changes to quasi-TE and quasi-TM. In agreement with  and the earlier symmetry argument, there is no coupling into modes of l ≠ 1 if the modification is centered (dx = 0). Figure 5(b) displays the corresponding spectra calculated from the coupling constants of 5(a) as outlined in section 4.2. The even l resonances appear in the spectra in between the odd l resonances (see Eq. (5)). Most striking is the sensitivity of coupling to l = 2 modes: even a small displacement dx ≈ 0.4 μm of the modification suffices for measurable coupling strength. The transmission spectra do not allow us to distinguish the strength of coupling to modes with l of the same parity, but the maps of coupling coefficient in (Fig 5(a)) indicate that measurable coupling to l = 3 and l = 4 modes can only be observed if the modification is at least 1 μm from the center.
In Figure 6 we compare measured and calculated transmission spectra for one FBG with the light launched perpendicular polarized to the fiber axis. The positions of the FBGs within the fiber core were estimated from microscope images. Due to fiber vibration in the ferrule, the micro voids are not uniformly distributed along the FBG, and the displacement values used in calculations are the maximum measured values. In the first case (Fig. 6(a) and 6(b), the spectrum was simulated for dx = 0.8 μm and dy = 0.7 μm. Both in experiment and simulation, the polarization direction was along the long axis of the micro void, parallel to the y- axis. Because of the displacement in y-direction, the perpendicular coupling constants (φ = −π/2 in Eq. (9)) had to be considered as well. The computed spectra closely agree with the experimentally obtained spectra. The slight difference in the envelope can be attributed to the restriction of the model to account only for uniform displacement. (Note that we have also performed coupling calculations using a core mode that was calculated using the finite element method to include the mean index perturbation associated with the grating modifications. This makes a negligible change to the transmission spectra.)
Due to the displacement from the center and the therefore asymmetric cross-section of the FBG, coupling to higher order modes occurs. This is clearly visible in the second comb of resonances in between the l = 1 resonances, which corresponds to the even l modes [see Eq. (5)]. Additionally, with the aid of imaging the reflected cladding modes, we were able to determine that interference between degenerate modes of different azimuthal order took place and therefore that higher order azimuthal modes were significantly excited. This will be the subject of a future publication.
Consistent with the results in Fig. 4, the low order resonances in Fig. 6(a) appear as single peaks but doublets appear at higher m (shorter wavelengths). The single peaks correspond to pure HE modes. The EH modes are formally present in these peaks, but as seen in Fig. 5(a), the coupling is too weak for them to be observed. While the single peaks show almost no polarization dependence, the strength of the doublet peaks heavily depends on how the axis of polarization is orientated with respect to FBG. The experimental and theoretical results concerning this aspect will also be presented in a future publication.
Aside from accurately rendering the spectral properties, our calculations (Fig. 5(a)) show the possibility of controlling coupling to certain classes of cladding modes by the position of the localized FBG. Suppressing coupling to certain modes is feasible by placing the micro-void on the symmetry axes of the kernel (Eq. (9)). However, the mode-canceling position might require a small off-set, depending on the asymmetries of the micro-void and its surrounding shell. Its exact value can be calculated by replacing the simplified void-model (Eq. (10)) with a more detailed refractive index profile.
In this paper we investigated the core-cladding coupling behavior of FBGs that have transversely inhomogeneous refractive index profiles. For that purpose, we wrote intra core FBGs in step index fibers based on the point-by-point approach. The micro voids formed within the index modifications that the grating comprises, are much smaller than the core and therefore allow for an accurate excitation of particular modes. Very strong coupling to cladding modes, exceeding 99 percent, could be achieved with such micro-void FBG. Up to now, efficiencies of this order of magnitude could only be obtained with heavily tilted FBG .
Furthermore, the high localization of the modification allows us to control the asymmetry of the FBG cross-section. These possibilities are evaluated within the three-layer-model of a circular fiber. Based on this model, we introduce a classification scheme of the cladding modes based on their overlap with the core mode. The concept of virtual cut-offs enables us to set a limit to the azimuthal mode number, depending on how broad a spectrum one wishes to investigate. Thus, multimode-spectra can be computed with a finite number of degenerate cladding modes. Even l resonances served as a strong test for the centeredness of a FBG. Our findings were validated by a rigorous calculation of the coupling constants. Accurate spectra of transversely inhomogeneous FBGs could be reproduced, which predicted many of the transmission characteristics of such FBGs, that have only been described qualitatively before. The necessary theoretical foundation for a detailed investigation of polarization dependent features could be established.
We have demonstrated efficient mode converters to odd and even l modes, realized with off-centered, highly localized FBGs. Such intra-core gratings allow for specifically targeting coupling to cladding modes of higher azimuthal index l. We believe these findings will enable many new fiber sensor designs, that rely on asymmetrically structured fiber cores. The efficient core-cladding mode conversion presented here can also be exploited for in-fiber elements that allow for the tailored excitation of higher-order mode fibers . Our extensive treatment can be easily transferred to fibers that carry multiple modes in the core for design and calculation of core-mode converting gratings.
A. Dispersion relation
Here we provide the expressions for the fields in the cladding and air regions. Inside the cladding (a1 ≤ r ≤ a2):
The financial support by the German Federal Ministry of Education and Research (BMBF) under Contract No. 13N9687 is gratefully acknowledged. Jens Thomas also acknowledges funding by the DAAD, grant D/0846673. CUDOS is an Australian Research Council Centre of Excellence.
References and links
1. T. Erdogan, “Cladding-mode resonances in short- and long-period fiber grating filters,” J. Opt. Soc. Am. A 14, 1760–1773 (1997). [CrossRef]
2. V. Mizrahi and J. E. Sipe, “Optical properties of photosensitive fiber phase gratings,” J. Lightwave Technol. 11, 1513–1517 (1993). [CrossRef]
4. T. Erdogan and J. E. Sipe, “Tilted fiber phase gratings,” J. Opt. Soc. Am. A 13, 296–313 (1996). [CrossRef]
5. C.-F. Chan, C. Chen, A. Jafari, A. Laronche, D. J. Thomson, and J. Albert, “Optical fiber refractometer using narrowband cladding-mode resonance shifts,” Appl. Opt. 46, 1142–1149 (2007). [CrossRef] [PubMed]
6. K. Zhou, L. Zhang, X. Chen, and I. Bennion, “Optic sensors of high refractive-index responsivity and low thermal cross sensitivity that use fiber Bragg gratings of 80° tilted structures,” Opt. Lett. 31, 1193–1195 (2006). [CrossRef] [PubMed]
7. T. Mizunami, T. Djambova, T. Niiho, and S. Gupta, “Bragg gratings in multimode and few-mode optical fibers,” J. Lightwave Technol. 18, 230–235 (2000). [CrossRef]
8. T. Guo, L. Shao, H.-Y. Tam, P. A. Krug, and J. Albert, “Tilted fiber grating accelerometer incorporating an abrupt biconical taper for cladding to core recoupling,” Opt. Express 17, 20651–20660 (2009). [CrossRef] [PubMed]
10. A. Martinez, M. Dubov, I. Khrushchev, and I. Bennion, “Direct writing of fibre Bragg gratings by femtosecond laser,” Electron. Lett. 40, 1170–1172 (2004). [CrossRef]
11. A. Martinez, M. Dubov, I. Khrushchev, and I. Bennion, “Photoinduced modifications in fiber gratings inscribed directly by infrared femtosecond irradiation,” IEEE Photon. Technol. Lett. 18, 2266–2268 (2006). [CrossRef]
12. G. Marshall, R. Williams, N. Jovanovic, M. J. Steel, and M. J. Withford, “Point-by-point written fiber-Bragg gratings and their application in complex grating designs,” Opt. Express 18, 19844–19859 (2010). [CrossRef] [PubMed]
13. J. Thomas, E. Wikszak, T. Clausnitzer, and U. Fuchs, “Inscription of fiber Bragg gratings with femtosecond pulses using a phase mask scanning technique,” Appl. Phys., A Mater. Sci. Process. 86, 153–157 (2007).
14. N. Jovanovic, J. Thomas, R. J. Williams, M. J. Steel, G. D. Marshall, A. Fuerbach, S. Nolte, A. Tünnermann, and M. J. Withford, “Polarization-dependent effects in point-by-point fiber Bragg gratings enable simple, linearly polarized fiber lasers,” Opt. Express 17, 6082–6095 (2009). [CrossRef] [PubMed]
16. A. Snyder and W. Young, “Modes of optical waveguides,” J. Opt. Soc. Am. 68, 297–309 (1978). [CrossRef]
17. S. J. Hewlett, J. D. Love, G. Meltz, T. J. Bailey, and W. W. Morey, “Cladding-mode coupling characteristics of Bragg gratings in depressed-cladding fibre,” Electron. Lett. 31, 820–822 (1995). [CrossRef]
18. S. Hewlett, J. Love, G. Meltz, T. Bailey, and W. Morey, “Coupling characteristics of photo-induced Bragg gratings in depressed-and matched-cladding fibre,” Opt. Quantum Electron. 28, 1641–1654 (1996). [CrossRef]
19. C. Tsao, D. Payne, and W. Gambling, “Modal characteristics of three-layered optical fiber waveguides: a modified approach,” J. Opt. Soc. Am. A 6, 555–563 (1989). [CrossRef]
21. B. E. A. Saleh, M. C. Teich, and J. W. Goodman, Fundamentals of photonics, (Wiley, 1991) pp. 272–309. [CrossRef]
22. H. Kogelnik, “Theory of dielectric waveguides,” Top. Appl. Phys. 7, 15–83 (1979).
23. D. Homoelle, S. Wielandy, A. L. Gaeta, N. F. Borrelli, and C. Smith, “Infraredphotosensitivity in silica glasses exposed to femtosecondlaser pulses,” Opt. Lett. 24, 1311–1313 (1999). [CrossRef]
24. S. Ramachandran, J. Fini, M. Mermelstein, J. W. Nicholson, S. Ghalmi, and M. F. Yan, “Ultra-large effective-area, higher-order mode fibers: A new strategy for high-power lasers,” Laser Photon. Rev. 2, 429–447 (2008). [CrossRef]