1. Introduction
Long period fiber gratings (LPGs) have attracted much attention in recent years and become the subject of intense research for sensing and communication due to peculiar advantages such as immunity to electromagnetic interference, high sensing resolution, low insertion loss, low back-reflections, compactness, and large-scale and low-cost fabrication, among many others [
1,
2,
3,
4,
5,
6,
7]. Currently, LPGs are important optical fiber passive components, and the accuracy of the numerical tools used for the calculation of wavelength-dependent transmittance has become a key factor, as the LPG operating sensing principle is closely related to a wavelength readout scheme [
8,
9].
In this paper we propose, for the first time to the best of our knowledge, the mathematical model of a device named as long period grating ring resonator (LPGRR), where the LPG is assumed to be inscribed in a fiber ring resonator (RR). In fact, LPGs have been always studied and discussed as stand-alone devices, in open-loop configurations, to operate as strain, temperature, and refractive index sensors with maximum sensitivities equal to 30.31 nm/1000 µε, 0.861 nm/°C, [
10], and 10
3–10
4 nm/RIU [
11,
12], respectively or as band-rejection filters [
13]. The only exception consists in a microwave photonic filter based on a fiber cladding-mode coupler, where the cladding-mode generated by an external LPG is injected into a RR and coupled out from the device where another identical LPG is placed outside the RR and used to convert the cladding mode into a core mode [
14].
The transfer matrix method (TMM) has been implemented to model and simulate the LPGRR assumed to be based on a three-layer step-index fiber geometry, in agreement with Erdogan’s theory used in this paper to calculate the cladding and core optical field components and core-to-cladding coupling coefficients in the LPG fiber [
15]. The mathematical model is demonstrated to be flexible, fast and accurate, and allows the coupled-mode theory (CMT) to be implemented for the design of directional couplers characterized by fiberoptic input/output bus and the LPGRR.
The effective refractive indexes of cladding and core fiberoptic modes have been calculated by using the finite element method (FEM) [
16], instead of the graphical method based on Bessel functions widely used in the literature [
15].
The paper is organized as follows. In
Section 2, the TMM-based modeling of the LPGRR is presented. In
Section 3, the simulation results are shown and discussed describing the LPGRR operation and unique spectral peculiarities. Further, LPG geometrical and optical features such as the length and the core-to-cladding coupling coefficients have been varied to perform parametric simulations. Theoretical results have been compared to an open-loop LPG experimental spectrum, demonstrating a very good agreement.
Section 4 investigates the LPGRR-based strain sensing performances as a function of different environmental conditions simulated as cover refractive index variations.
Section 5 summarizes the conclusions.
2. Modeling
The fiberoptic resonant architecture consists of a fiber racetrack cavity in which a conventional LPGRR is included. In particular, LPGRR can be simply modeled by the TMM, scattering matrix and CMT. Further, a linearly polarized field that corresponds to the core mode is coupled to the RR by means of a fiber directional coupler and to the cladding mode by the LPG, provided that the phase matching condition (PMC) is satisfied [
13,
15].
The fiber RR is modeled in
Figure 1, where the arrows indicate the optical signal propagation path,
R is the ring radius, and
τ and
κ are the directional coupler transmission and coupling factors, respectively. In particular, several assumptions can be made to investigate the device operation and performance. For example, the excitation of a propagating unidirectional resonant single mode within only one polarization state, as well as lossless coupling, is assumed. Further, it is supposed that none of the waveguide segments and coupler elements can couple waves of different polarization states, and that all the loss contributions occurring along the light propagation in the ring resonator filter are included in the attenuation constant. In particular, the field profiles and the effective refractive indexes of the supermodes associated with the fundamental and selected cladding modes allows the analytical calculation of the complex transmission and coupling factors
τ and
κ to be performed by evaluating the overlap integrals at their respective ports [
17,
18,
19,
20].
Figure 2 shows the electric components of the TE symmetric and antisymmetric supermodes of the directional coupler calculated by a full-vectorial finite element method and compares the supermode profiles obtained directly from the mode solver for
LP01 core and
LP08 cladding modes. Further details on fiberoptics geometrical and optical parameters are reported in
Table 1 (
Section 3).
The total transfer matrix is given by the product of the matrices corresponding to each section in which the device has been divided, in the reverse numbered order as indicated by the subscripts in Equation (1) and shown in
Figure 3:
In Equation (1),
T2,lpg is the LPG scattering matrix and
T1 =
T3 are the scattering matrices associated with the highlighted sections of the cavity (see
Figure 3).
Figure 3 also shows the four-port representation of the LPGRR with inputs and outputs field components used to develop the analysis, i.e.,
ei and
fi the core components,
gi and
hi the cladding components at the
i-th port (
i = 1, 2, 3, 4) with the directions of propagation indicated by the arrows.
The individual scattering matrices are defined as in Equations (2) and (3):
where
O is a 2 × 2 zero matrix,
TLPG is a 2 × 2 matrix defined as in Equation (4):
and
In Equations (4) and (5),
is the grating phase,
, with
δ the detuning parameter,
Llpg the LPG length,
kco-cl is the core-to-cladding mode coupling coefficient,
βco and
βcl are the propagation constants of core and cladding modes, respectively [
15]. Moreover, Λ is the grating period,
α represents the cavity losses and
L is the length of the orange and yellow colored bend and straight symmetric sections calculated as in Equation (6):
In Equation (6),
Ls is the physical length of the coupler,
Lc is the coupling length, and
R is the racetrack radius as already defined in
Figure 3.
In the coupling region, output fields (i.e.,
fi and
hi) are related to the input fields (i.e.,
ei and
gi) by the transmission and coupling coefficients as represented by the Equations (7)–(10), where τ, τ
c and κ, κ
c are the transmission and coupling coefficients for core and cladding modes, respectively.
In particular, by defining
, and
, and considering the definition of the LPGRR transfer matrix (i.e., Equation (1)), the incoming and outgoing fields at the 3rd port can be calculated as in Equations (11) and (12), respectively:
Then, by substituting Equation (10) in Equations (11) and (12), Equations (13) and (14) can be derived:
Finally, Equation (15) is also reported as calculated by Equation (14):
where
,
,
, resulting in
and
, where
, and
,
.
By imposing the initial conditions
e1 = 1,
e2 =
g1 =
g2 = 0 (note that because of these assumptions
Tlpg,
Ti,
Tfiber and the matrices in Equations (7)–(10) are all diagonal), which means that there is only one input port and the input field is only that of the core mode (the cladding mode is generated by the LPG), the reflected components are given by Equation (16):
The transmitted components can be calculated as in Equation (17):
3. Results
A 3 cm-long LPG with a period of 311 µm and an average index change of the grating area equal to σ = 4 × 10
−4 is assumed to be written in a quarter of a single-mode optical fiber ring resonator. This period satisfies the PMC [
15] for LP
08 cladding mode at the resonance wavelength of 1539.8 nm. The SMF-28 fiber and LPGRR design parameters used in the simulations are summarized in
Table 1, while the transmission and coupling coefficients for both modes are obtained by using the Supermodes Theory in a 6.55 cm-long coupling region and listed in
Table 2. The LPG has been written in the SMF-28 fiber by classical UV laser exposure through an amplitude mask. Experimental results have been achieved by measuring the transmission spectrum by an OSA, are plotted in
Figure 4 and compared with theory.
A good agreement between theoretical and experimental spectral distributions of the open-loop LPG has been achieved as demonstrated in
Figure 4, where a fitting core-to-cladding coefficient equal to 50 m
−1 has been used in accordance with the LPG fabrication. Actually, a disalignment between experimental and theoretical LPG spectral resonances can be observed. However, some reasons for this behaviour can be identified. First, the theoretical calculations are based on the weakly guiding approximation that predicts the dispersion of the core and cladding modes with a limited accuracy. Second, the average value of the UV-induced grating was not rigorously constant, as assumed in the calculations, due to slightly different UV exposure times. Nevertheless, the experimental agreement with the predicted trends is quite good, with only very small discrepancies, quantifiable in a Δλ respectively equal to ~6 nm for the short wavelengths (e.g., 1250 and 1350 nm) and ~10.2 nm at 1550 nm.
Under these assumptions, the core and cladding spectral responses at the output of the LPGRR can be observed in
Figure 5a,b.
The peculiarity of this device is the presence of two resonance peaks, one deeper than the other in the core output spectrum and, dually, one higher than the other in the cladding one. This spectral feature is indicated as splitting resonances (SR), namely due to the combination of resonance effects of LPG and RR on the input signal. This effect can be clearly seen in
Figure 6, where the LPGRR core and cladding transmission spectra are zoomed in around 1539 nm. The specific transmission spectra originate from the combination of the resonance characteristics of both RR and LPG. Note that the core and cladding modes outputs are dual because of their mutual power exchange.
RR filters can be described by specific figures of merit, such as the resonance width, defined as the full-width-at-half-maximum (FWHM) and the distance between resonance peaks, i.e., the free spectral range (FSR). A graphical estimation of these parameters has been performed for the LPGRR (see
Figure 5a) and for a standard RR characterized by the same length. Numerical parameters are listed in
Table 3. In particular, referring to the FWHM values, the width of the shallowest LPGRR core spectrum splitting resonance peak (
Figure 5a) is 39.5% thinner than the standard RR resonances without any LPG. Furthermore, the width of the deepest LPGRR core spectrum splitting resonance peak is 44.2% thinner than standard RR resonances without LPG. Finally, the distance between the splitting resonances, i.e., labeled as splitting distance in
Figure 5a,b, is equal to 4.03 pm.
Parametric Analysis
It is worth noting that the spectral responses, in particular the distance between the splitting peaks and their amplitudes, change by varying the LPG length and the core-to-cladding coupling coefficient, k
co-cl. The simulations reported in
Figure 7 and
Figure 8 refer to the core spectrum centered at ~1539 nm, but still remain valid in the case of the cladding mode.
Table 4 lists the calculated splitting resonance distance when the grating parameters are changed.
In
Figure 7 and
Table 4, the grating length of the LPG influences the distance between the splitting resonances,
dspl; it increases linearly by increasing the grating length. Instead, a linear trend cannot be observed when the core and cladding mode splitting distance and coupling coefficient are considered, as also evidenced in
Figure 8 and
Table 4.
5. Conclusions
This paper explores, for the first time to the best of our knowledge, an optical device based on a long period fiber grating (LPG) inscribed in a section of a conventional fiberoptic racetrack resonator, namely an LPGRR. In particular, the LPGRR is fully and generically modelled by means of a TMM-based mathematical model which is a flexible and robust tool that allows optical and structural LPGRR changes to be easily simulated.
Spectral distributions of an open-loop LPG with a length of 3 cm and a period of 311 µm have been demonstrated, achieving a very good agreement between theory and experiments. Furthermore, the combination of the RR and the LPG generates specific transmittance spectra characterized by the presence of two resonance peaks, indicated as splitting resonances, separated by a spectral distance of 4.03 pm at ~1539 nm. This unique spectral operation has been simulated by varying LPG optical and geometrical features in order to perform a parametric analysis. In particular, the LPG length varied from 2 cm to 5.5 cm, resulting in an increase of the distance between the splitting resonances from 3.75 pm to 4.43 pm, respectively. Furthermore, by varying the core-to-cladding coupling coefficient from 50 m−1 to 300 m−1 the splitting resonances distance varies in a nonlinear way but in the proximity of 4 pm.
The LPGRR spectral response can be the subject of future studies and be employed for sensing applications. For example, the distance between the splitting resonances or the attenuation band shift can be used to perform sensing functionalities. In particular, the LPGRR-based strain sensing performances have been investigated in different environmental conditions with a core-to-cladding coupling coefficient varied as a function of the cover RI variations. A strain sensitivity of ~1.8 × 10−3 nm/µɛ can be theoretically achieved with an applied strain ranging from 0 μɛ to 1 μɛ. Finally, as negligible slope variations of the sensitivity curves have resulted, the sensing performances can be still preserved when the LPGRR operates in different locations. In conclusion, the LPGRR strain sensor can exhibit a linear operation and immunity from external environment conditions which can be extremely advantageous for several sensing applications.