Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Enhanced refractometer for aqueous solutions based on perfluorinated polymer optical fibres

Open Access Open Access

Abstract

The use of the new CYTOP (Cyclized Transparent Optical Polymer) fibres for the inscription of optical structures and the detection of different parameters has started to gain importance in the past decade. This work presents the design, simulation and manufacture of a CYTOP-based surrounding refractive index sensor for aqueous solutions, given its high sensitivity in the range 1.315 − 1.333 (at 1550 nm wavelength). The structure is based on a bent and polished fibre (in order to increase its sensitivity), the polished area being the surface on which a diffraction grating is inscribed with a femtosecond laser. The interaction of the field propagated by the fibre with the grating causes diffraction of certain orders towards the outside, depending, among other things, on the refractive index of the fluid. In addition to a maximum sensitivity of −208.8 nm/RIU and a remarkable insensitivity to temperature, it offers a spectral fingerprint of each sensed fluid.

© 2022 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement

1. Introduction

A great attention and research effort have been drawn to the transducing properties of optical fibres, which has made possible a breakthrough in the last $50$ years with respect to optical fibre sensors (OFSs), whether in glass fibre [1] or plastic optical fibre (POF) [2,3]. Its electrical immunity, small transmission loss, remote sensing, and high sensitivity, outperform other more traditional schemes. Moreover, the discovery of new materials in fibres, which improve their optical properties, allow their application in increasingly varied fields, whether in the scope of engineering, industry, medicine, or biochemistry, among others.

Regarding plastic optical fibres, poly(methyl methacrylate) (PMMA) has traditionally been used. The first fibres based on PMMA cores emerged in the $1960$s, and their growth has been extraordinary to this day. In contrast to glass fibres, POF stands out for its ease of handling, flexibility, and reduced price [2,4]. Likewise, the improvement of the manufacturing process has considerably limited their attenuation, and they can be used for short-haul communications or optical fibre sensors.

However, in the late $1990$s Asahi Glass Co (currently ACG Inc.) developed a special type of POF, based on an amorphous fluorinated polymer called CYTOP [5,6]. It is poly(perfluoro-butenylvinyl ether). Characterized by a refractive index of $\sim 1.34$, it stands out for its compatibility with the aqueous index regime, an important feature for biosensing [7]. Moreover, it has a very low near-infrared attenuation ($\sim 10$ dB/km) and a transmission bandwidth of up to $1300$ nm. It is a graded-index polymer optical fibre (GI-POF), with less modal noise and material dispersion than PMMA fibers [2]. The use of CYTOP fibres began to gain important relevance in the $2010$s (Table 1), so it is a field with a great future projection. In recent years, the inscription of optical structures using femtosecond lasers in CYTOP fibres for the development of OFSs stands out [8]. Table 1 lists the most important works in the literature. As it can be seen, the works are recent, and most of them use fibre Bragg gratings (FBGs) as the sensing structure.

Tables Icon

Table 1. Works related to the manufacture of optical structures in CYTOP fibres by means of femtosecond laser.

In this work, an optical structure based on a perfluorinated POF is designed, simulated and manufactured using a femtosecond laser. Taking advantage of its infrared transmission, the device is used as a highly-sensitive refractive index sensor in aqueous solutions, thus having potential applicability in the field of biomedicine [16]. To this end, in addition to increasing sensitivity to the external medium by bending the fibre and polishing its cladding, a diffraction grating is inscribed into the previously polished area. The variability of the transmitted diffraction orders provides a spectral fingerprint of the different aqueous solutions, which can be used for the sensing of their refractive index, among others.

2. Optical structure design

The designed optical structure is the one depicted in Fig. 1. The dimensions of the CYTOP fibre used are listed below. The graded-index core has a diameter $d_{core}=120\;\mathrm{\mu }$m, while the cladding has a thickness of $20\;\mathrm{\mu }$m ($d_{clad}=160\;\mathrm{\mu }$m). Likewise, there is a polycarbonate protection overcladding with a total diameter $d_{oc}=490\;\mathrm{\mu }$m. The fibre, measured in transmission, was placed with a curvature defined by a radius of curvature (ROC) $R$. To fix the fibre, two PMMA bulks were manufactured with the groove where the fibre was placed, closing the sandwich-structure with epoxy resin.

 figure: Fig. 1.

Fig. 1. Diagram of the CYTOP-based optical structure.

Download Full Size | PDF

Subsequently, the PMMA bulk corresponding to $\Psi =90^\circ$ started to be polished (first with ultra-fine grit sandpaper, P4000; and then with $12\; \mu$m lapping film). $\epsilon$ refers to the depth of polish. At the moment in which the outer interface of the fibre is reached, $\epsilon =0$. When polishing continued ($\epsilon >0$), the polished part of the fibre began to generate an elliptical shape in the $XY$ plane. To determine the polishing depth ($\epsilon$), the dimensions of the ellipse generated were considered, specifically their minor and major axes. The geometric and trigonometric relationships that relate the different design parameters are presented in Fig. 2(a) (for the major axis) and Fig. 2(b) (minor axis) [17]. The expressions that relate $\epsilon$ and the axes of the ellipse ($x$, $y$) are shown below:

$$y^2+\left(R+\frac{d_{oc}}{2}-\epsilon\right)^2=\left(R+\frac{d_{oc}}{2}\right)^2,$$
$$x^2+\left(\frac{d_{oc}}{2}-\epsilon\right)^2=\left(\frac{d_{oc}}{2}\right)^2.$$

 figure: Fig. 2.

Fig. 2. Geometric parameters that define the (a) major axis and (b) minor axis of the ellipse generated in the XY plane (Fig. 1).

Download Full Size | PDF

3. BPM simulations

Diffraction grating aside, there are two design parameters: $R$ and $\epsilon$. To obtain the optimal values, finite element simulations are carried out using the beam propagation method (BPM) with the commercial software RSoft. First, a 3-D simulation (in order to obtain the spatial-pattern of the cross-section) is performed by sweeping the radius of curvature ($R$). Subsequently, once the $R$ value has been chosen, the sweep is performed on the polishing depth ($\epsilon$). In the latter case, by focusing only on the final fibre losses, a 2-D simulation is sufficient. The simulation parameters used in both cases are listed in Table 2.

Tables Icon

Table 2. Simulation parameters used in RSoft for both cases ($R-$sweep and $\epsilon -$sweep).a

It should be noted that in both cases, in order to make the bending simulation feasible in BPM, a mathematical coordinate transformation called conformal mapping is used [18,19]. The transformation used, that preserves the angle between oriented curves, is as follows:

$$W=u+iv=f(Z)=f(x+iy)=R\cdot\ln\left(\frac{Z}{R}\right),$$
for which
$$\left|\frac{\textrm{d}{Z}}{\textrm{d}{W}}\right|=\exp\left(\frac{u}{R}\right).$$

In this way, the fibre curvature can be modelled from a tilted refractive index profile, achieving a translation invariance in light propagation direction. The tilted RI profile ($n_{tilt}$) as a consequence of the transformation is obtained by multiplying the actual RI of the material ($n_{real}$) by Eq. (4):

$$\begin{aligned}n_{tilt}(u,v)=n_{real}(u,v)\cdot\exp\left(\frac{u}{R}\right)&{\approx}n_{real}(u,v)\left(1+\frac{u}{R}\right). \\ &\uparrow \\ d_{clad}&\ll R\end{aligned}$$

Figure 3 depicts, for both unpolished ($R-$sweep, Fig. 3(a)) and polished ($\epsilon -$sweep, Fig. 3(b)) fibre, a schematic of the fibre in the $Z$-plane (no transformation) and in the $W$-plane (conformal mapping applied), as well as the corresponding refractive index profiles. The RI profile in Fig. 3(b) corresponds to a polished fibre in the $\Psi =90^\circ$ position. Fibre polishing centered at $\Psi =90^\circ$ is mathematically associated with the function shown below (defined in the Z-plane), in polar coordinates. The limits of $\theta$ are defined from trigonometric relations obtained according to Fig. 2(a).

$$r(\theta)=\frac{R+\frac{d_{clad}}{2}-\epsilon}{\sin\theta},\quad\quad \theta\in\left[\frac{\pi}{2}-\alpha,\frac{\pi}{2}+\alpha\right],$$
with $\alpha =\cos ^{-1}\left (\frac {R+\frac {d_{clad}}{2}-\epsilon }{R+\frac {d_{clad}}{2}}\right )$. Accordingly, $\delta =2R\alpha$.

 figure: Fig. 3.

Fig. 3. Schematic of the fibre in the $Z$-plane (no transformation) and in the $W$-plane (conformal mapping applied) for (a) unpolished ($R-$sweep simulation), (b) and polished (at $\Psi =90^\circ$) fibre ($\epsilon -$sweep simulation). Refractive index profiles are depicted for each case (with and without transformation).

Download Full Size | PDF

Once conformal mapping has been defined, the simulations are carried out. First of all, regarding the ROC, a value of $R$ must be chosen such that it causes relevant bending losses. In this case, an increase in attenuation is the basis of many optical sensors, due to the high sensitivity achieved to detect variations in the surrounding medium. Figure 4 depicts the simulated cross-section of the fibre (core and cladding) at $\Psi =90^\circ$ (polishing position) for different values of $R$.

 figure: Fig. 4.

Fig. 4. RSoft simulation of the resulting field in the fibre cross-section at $\Psi =90^\circ$ for bending radii of (a) $2$, (b) $4$, (c) $5$, (d) $6$, (e) $10$ and (f) $20$ mm. The attenuation obtained at the end of the curvature ($\Psi =180^\circ$) for each radius is $31$, $16.3$ $12.7$, $9.7$, $5.1$ and $2.7$ dB, respectively. The boundaries corresponding to core (dotted line) and cladding (solid line) are also depicted.

Download Full Size | PDF

It can be seen that extremely small $R$ values (e.g. $2$ mm) result in excessive sensing attenuation ($>15$ dB), as well as greater ease of fibre breakage. On the other hand, bending radii larger than $6$ mm make the power propagate more easily along the axial axis of the fibre (GI-POF), with less dependence of the bending losses on the outer refractive index. Consequently, the optimum bending radii are around $5$ mm, which is the value chosen. For this value, a noticeable presence of propagated cladding modes can be observed, which, after polishing, will be dependent on the surrounding refractive index (SRI).

For the chosen bending radius value ($R=5$ mm), the fibre losses $L$ (at $\Psi =180^\circ$) are simulated for the values given in the corresponding part of Table 2. Specifically, a sweep is performed on the surrounding refractive index (SRI between $1$ and $1.36$) for different values of polishing depth ($\epsilon$ between $0\;\mathrm{\mu}$m and $60\;\mathrm{\mu }$m with steps of $10\;\mathrm{\mu }$m). The results are plotted in Fig. 5.

 figure: Fig. 5.

Fig. 5. (a) RSoft simulation of fibre losses ($L$) due to curvature ($R=5$ mm) as a function of polishing depth ($\epsilon$) and surrounding refractive index. (b) $\partial L/\partial \textrm {SRI}$ of the values represented in (a). The SRI range between $1.318$ and $1.336$ is highlighted, due to its remarkable sensitivity.

Download Full Size | PDF

It can be seen that, together with the fibre curvature itself, the polishing of the cladding contributes to increase the sensitivity to the external medium. However, due to its graded-index refractive index profile, the transverse field travels more confined to the axial axis of the fibre. Consequently, there are no evident differences with increasing polishing depth (only a slight increase in losses), even going into the core ($\epsilon \geq 20\;\mathrm{\mu }$m). This differs from using a step-index fibre [17]. A polish of about $50\;\mathrm{\mu }$m, slightly inside the core, will be a good value to be able to inscribe a sufficient number of periods on the diffraction grating, as detailed in the next section.

4. Diffraction grating operating principle

A diffraction grating will be inscribed on the polished surface. The angles of the diffracted modes ($\theta _d$ for transmission, and $\theta _r$ for reflection) are related to the wavelength ($\lambda$) and grating period ($\Lambda$) through the grating equation, both in the transmission (Eq. (7)) and the reflection region (Eq. (8)).

$$n_L\sin\theta_d=n_{core}\sin\theta_i+\frac{\lambda}{\Lambda}p,\quad\quad p\in\mathbb{Z}.$$
$$n_{core}\sin\theta_r=n_{core}\sin\theta_i+\frac{\lambda}{\Lambda}p,\quad\quad p\in\mathbb{Z}.$$

$\theta _i$ refers to the angle of the incident wave, $p$ to the diffraction order, $n_{core}$ to the core RI, and $n_L$ to the SRI. These equations only predicts the direction of the modes, not how much power is in them. Figure 6 depicts a schematic of the above.

 figure: Fig. 6.

Fig. 6. Diagram of 1-D diffraction grating. The diffraction orders in transmission and reflection are depicted. It is assumed that there is no cladding (removed in polishing).

Download Full Size | PDF

The grating period chosen must be such as to ensure a low-order grating ($\Lambda >\lambda /n_{core}$) with relevant efficiency in the diffracted orders. For the spectral range to be used ($\lambda \in [1500-1600]\textrm { nm}$), a grating period $\Lambda =2\;\mathrm{\mu }$m will be used. According to the selected value $\epsilon =50\;\mathrm{\mu }$m, it is obtained an ellipse with major and minor axes of $1.42$ mm and $148\;\mathrm{\mu }$m, respectively. This results in the inscription of a grating with more than $700$ periods, thus improving the efficiency of the different diffraction orders.

In this case, the appearance of new diffraction orders (as a function of $\lambda$ and $n_L$) causes a transfer of energy from the old to the new diffraction orders (Rayleigh anomalies [20,21]). Consequently, there will be transitions or power variations (absorption peaks in short) in the spectrum transmitted and collected in an optical spectrum analyser (OSA). This is the main feature to be exploited in the measurement process.

5. Experimental results and discussion

5.1 Fs laser setup and inscription parameters

The inscription of the diffraction grating has been carried out using the optical setup depicted in Fig. 7. A commercial fibre laser chirped pulse amplifier (FLCPA) from CALMAR laser has been applied. It has a $1030$ nm wavelength, and $370$ fs pulse duration. However, the laser beam is passed through a barium borate crystal (BBO) in order to generate a $515$ nm laser beam (Second Harmonic Generation, SHG). After harmonic separation is performed, the previously reduced beam (higher irradiance in the BBO results in higher conversion efficiency) is expanded. The pulse energy ($E_p$) is controlled by the combination of $\lambda /2$ plate and linear polarizer. Subsequently, a half-wave plate and a quarter wave-plate are arranged in order to have total control of the polarization state of the laser beam. Subsequently, an iris diaphragm is placed to provide a clean wavefront. The laser pulses are tightly focused through a $\textrm {NA}=0.5$, $\times 100$ Mitutoyo objective lens. The sample is placed on a nano-resolution XYZ motor stage from Aerotech. Finally, the monitoring and visualization of the inscription process is carried out with a CMOS camera.

 figure: Fig. 7.

Fig. 7. Femtosecond laser inscription setup.

Download Full Size | PDF

Once the structure has been polished with ultra-fine vegetable fibre sandpaper, the surface depicted in Fig. 8(a) is obtained. The area corresponding to the cladding corresponds to an ellipse with major and minor axes of $1.41$ mm and $148\;\mathrm{\mu }$m, respectively. According to Eq. (1) and Eq. (2), the $\epsilon$ value is $49.2\;\mathrm{\mu }$m. Fs laser inscription parameters of the diffraction grating are $\Lambda =2\;\mathrm{\mu }$m, grating length of $1.45$ mm, grating line length of $160\;\mathrm{\mu }$m, $E_p=220$ nJ, pulse repetition rate (PRR) of $30$ kHz and $1000$ pulses/µm. Figure 8(c) depicts an image of the diffraction grating, with the final complete structure in Fig. 8(d). With the indicated protocol, the repeatability of the sensor is assured.

 figure: Fig. 8.

Fig. 8. (a) $5\times$ microscope image of the ellipse generated as a result of polishing, showing both cladding and core, as well as overcladding. $20\times$ microscope images (b) before (c) and after diffraction grating inscription of the region indicated in (a) are also shown. (d) Image of the complete CYTOP-based optical structure.

Download Full Size | PDF

5.2 Sensor application

The applicability of the fabricated structure as a refractometer for aqueous solutions is then verified. For this purpose, transmission measurements are performed using the super-luminescent diode (SLD) light source S5FC1005S (Thorlabs), and the OSA MS9740A (Anritsu), with $30$ pm resolution. There is approximately $\sim 1$ m of CYTOP fibre, with single-mode fibre (SMF) connections at both ends.

The liquids used in the measurement process are obtained using different ratios of sucrose-aqueous solution. The RI determination of those liquids has been carried out with the Krüss DR6100 refractometer, at a temperature of $23.9\pm 0.4^\circ$C. The RI values, both measured in the yellow doublet D-line of sodium ($589.3$ nm) and estimated in C-Band ($1550$ nm), are presented in Table 3. The estimation is made by combining existing literature data [22,23], and Sellmeier coefficients in water [24]. According to the simulation results in Fig. 5, the liquids used are in the RI range of highest sensitivity.

Tables Icon

Table 3. Solute concentration, $n_D$ and $n(1550\textrm { nm})$ of the liquids used in the experiments.

Placing the manufactured structure in a fixed position, it is immersed in the different liquids, obtaining the transmission spectra depicted in Fig. 9(a). The measurements are repeated five times, with the uncertainty ($95\%$ confidence level, CL) depicted in the shaded region of each curve. Processing the spectra, it is observed that there are six points whose tracking consistently gives a discrimination of $n_L$ with a high sensitivity in the range $1.315-1.33$. These are valleys (absorption peaks) $2$, $4$, $5$ and $6$ in the figure, and peaks $1$ and $3$. Figure 9(b) contains the wavelength shift ($\Delta \lambda$) of each of these peaks/valleys for the six liquids used, with the corresponding quadratic fit. In addition to the high coefficient of determination ($R^2$) of the fits, it is worth noting that there are three obvious trends, i.e. peaks with similar behaviour: $1$, $2$ and $6$ have a higher sensitivity in the range $n_L\in [1.325-1.33]$, $4$ has its highest sensitivity in the range $1.315-1.32$, while $3$ and $5$, although they have a lower sensitivity, have a higher linearity over the whole range $1.315-1.33$ (Table 4).

 figure: Fig. 9.

Fig. 9. (a) Transmission spectra when the structure is immersed in liquids with different RI. The $6$ peaks and valleys on which the tracking in sensing is performed are highlighted. (b) Wavelength shift ($\Delta \lambda$) of each of the $6$ peaks/valleys for the liquids used, with their quadratic fit and error bars ($\textrm {CL}=95\%$). (c) Simulation and experimental transmission losses for $\lambda =1550$ nm, as well as experimental results for $\lambda =1500$ nm and $\lambda =1600$ nm, with the associated error bars ($\textrm {CL}=95\%$)

Download Full Size | PDF

Tables Icon

Table 4. Parameters of the quadratic and linear fit of the six curves depicted in Fig. 9(b).

On the other hand, the fibre losses are also depicted in Fig. 9(c). For $\lambda =1550$ nm, both the experimental results and their quadratic fit, as well as the simulation results (Fig. 5(a)), are shown. Taking into account that the intrinsic losses of the CYTOP fibre have not been taken into account in the simulation, the similarity is remarkable. In fact, the Pearson correlation coefficient between the simulation curve for $n_L\in [1.315-1.333]$ and the quadratic fit of the experimental results is $\rho =0.9799$. In this range, the sensitivities are $695.9$ dB/RIU (experimental measurements) and $422.5$ dB/RIU (simulation). Finally, the experimentally obtained losses for the $1500$ and $1600$ nm wavelengths are also depicted in order to observe the trend. Considering that sensitivities above $120$ nm/RIU are obtained in any of the three defined ranges (Table 4), it is possible to sense changes in SRI when they are in the order of $10^{-3}$ ($\Delta \lambda \geq 0.12$ nm) or even $2.5\cdot 10^{-4}$ ($\Delta \lambda \geq 30$ pm). Only peak $\#4$ presents an uncertainty higher than the OSA resolution ($30$ pm), especially in the first of the linear sections defined in Table 4. In that range ($n_L\in [1.315-1.32]$), the detection occurs for $\Delta \lambda \geq 57$ pm, but the sensitivity in that section is $-179.7$ nm/RIU, so it is possible to detect changes in SRI of up to $\sim 3.1\cdot 10^{-4}$. Peaks with a positive concavity, such as $\#4$, are used as sensing points for the lower range of the sensing zone, due to their high sensitivity.

5.3 Temperature cross-sensitivity

Temperature cross-sensitivity has an important influence in terms of refractive index measurement. The characteristics of the diffraction grating are modified with temperature, as is the period ($\Lambda$). Firstly, it is important to note that the coefficient of linear thermal expansion (CLTE) of CYTOP is $\sim 7.4\cdot 10^{-5}\; ^\circ \textrm {C}^{-1}$, which is slightly lower than PMMA. Therefore, it is not a relevant factor here, even less so in a temperature range below $30\; ^\circ \textrm {C}$. To evaluate the behaviour of the SRI sensor, a temperature sweep between $23^\circ$C (room temperature) and $40^\circ$C is performed. In the measurement process, the device is immersed in water ($n_L=1.3153$), and the peaks/valleys indicated in Fig. 9(a) are tracked. The variation of the refractive index of the liquid in this temperature range is less than $5\cdot 10^{-4}$. The results obtained are depicted in Fig. 10. The most pronounced variation is observed at peak $\#1$, with $1.67$ pm$/^\circ$C. In any case, this is a value that can be considered negligible, since it gives rise to an error in the SRI measurement that is in the order of $10^{-5}$. It can therefore be concluded that the measurement of the surrounding refractive index is insensitive to the temperature of the fluid.

 figure: Fig. 10.

Fig. 10. Wavelength shift ($\Delta \lambda$) of each peak/valley as a function of temperature. The structure has been immersed in water ($n_L=1.3153$), and a temperature sweep between $23$ and $40$ $^\circ$C has been performed. A linear fit is made, with the slope of each fit. It should be noted that the wavelength step in the spectrum is $50$ pm.

Download Full Size | PDF

6. Discussion

To conclude the paper, it is interesting to briefly outline the main innovations and advantages compared to existing work on fibre-based refractometers. In a $2019$ review, Urrutia et al. present the state of the art associated with optical fibre refractometers [25], classifying them according to whether they are interferometers, grating-based structures, or resonance-based structures (coating). Although it can be concluded that there is not a great difference in performance between these three groups, it can be highlighted that tilted fibre Bragg gratings (TFBG) have been traditionally used as they are absolute refractometers in a wide wavelength range [26,27]. Moreover, they are relatively robust and simple, as they do not need specialty fibres, nor modifications of their geometry that degrade their practical use from a commercial point of view. With an accuracy of $\pm 5\cdot 10^{-5}$, the grating-based refractometer presented in [26] indicates that works are beginning to exist that obtain comparable performances to the Abbe refractometer.

TFBG-based refractometers have even been developed on polymer fibres [28]. These are custom PMMA fibers with a core diameter of $8.2\;\mathrm{\mu }$m and $\textrm {NA}=0.23$, which gives rise to a few-mode behaviour ($4$ guided modes). Although the performance is lower than silica single-mode fibres, sensitivities of up to $\sim 13$ nm/RIU are achieved. They also allow to increase their high sensitivity measurement range up to $\textrm {SRI}=1.49$, compared to silica fibers, with a range of up to about $\textrm {SRI}=1.45$. However, commercially available perfluorinated fibres have a highly multimode behaviour, which makes it impossible to use TFBGs to detect SRI with high accuracy. The great advantage of using polishing in combination with diffraction gratings is that, by eliminating cladding in a curved area, it is possible to induce most of the guided modes to interact with the external medium and condition their propagation as a function of SRI. Also, although these fibres do not allow a high sensitivity region as wide as PMMA fibres, they do have a maximum sensitivity region of around $\textrm {SRI}=1.34$. This makes them very suitable for applications that require a high sensitivity region. This makes them very suitable for applications requiring discrimination of aqueous solutions.

7. Conclusion

To the authors’ knowledge, this is the first work using perfluorinated polymer fibres for the sensing of the surrounding refractive index. It also presents the novelty of using an in-fibre diffraction grating for the detection. This optical structure, inscribed with a femtosecond laser on a previously polished curvature zone, makes it possible to modify the diffraction orders transmitted (losses) and reflected (propagated by the fibre) as a function of the SRI, thus obtaining a specific “fingerprint” of the external liquid. Due to the RI of the fibre used, this is a device designed for use in aqueous solutions, with a potential use in biomedical applications. Also, as can be deduced from its final presentation, depicted in several images of the work, it is easily encapsulated for a robust practical use. It has a maximum sensitivity of $-208.8$ nm/RIU, can detect changes in SRI of around $2.5\cdot 10^{-4}$, and is insensitive to temperature, thus avoiding cross-sensitivity issues.

Funding

Ministerio de Ciencia, Innovación y Universidades, Agencia Estatal de Investigación and FEDER funds of European Commission (PID2019-107270RB-C21 by MCIN/AEI/10.13039/501100011033); Ministerio de Educación , Cultura y Deporte of Spain (PhD grant FPU2018/02797).

Disclosures

The authors declare no conflicts of interest.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

References

1. D. Pallarés-Aldeiturriaga, P. Roldán-Varona, L. Rodríguez-Cobo, and J. M. López-Higuera, “Optical fiber sensors by direct laser processing: A review,” Sensors 20(23), 6971 (2020). [CrossRef]  

2. K. Peters, “Polymer optical fiber sensors—a review,” Smart Mater. Struct. 20(1), 013002 (2011). [CrossRef]  

3. L. Bilro, N. Alberto, J. L. Pinto, and R. Nogueira, “Optical sensors based on plastic fibers,” Sensors 12(9), 12184–12207 (2012). [CrossRef]  

4. J. Zubia and J. Arrue, “Plastic optical fibers: An introduction to their technological processes and applications,” Opt. Fiber Technol. 7(2), 101–140 (2001). [CrossRef]  

5. N. Yoshihara, “Low-loss, high-bandwidth fluorinated POF for visible to 1.3-μm wavelengths,” in OFC ’98. Optical Fiber Communication Conference and Exhibit. Technical Digest. Conference Edition. 1998 OSA Technical Digest Series Vol.2 (IEEE Cat. No.98CH36177), (Opt. Soc. America, 1998).

6. M. Naritomi, “CYTOP amorphous perfluoropolymer for low loss POF,” POF Asia Pacific Forum ’96 (1996).

7. A. Lacraz, M. Polis, A. Theodosiou, C. Koutsides, and K. Kalli, “Femtosecond laser inscribed Bragg gratings in low loss CYTOP polymer optical fiber,” IEEE Photonics Technol. Lett. 27(7), 693–696 (2015). [CrossRef]  

8. A. Theodosiou and K. Kalli, “Recent trends and advances of fibre Bragg grating sensors in CYTOP polymer optical fibres,” Opt. Fiber Technol. 54, 102079 (2020). [CrossRef]  

9. Y.-G. Nan, D. Kinet, K. Chah, I. Chapalo, C. Caucheteur, and P. Mégret, “Ultra-fast fiber Bragg grating inscription in CYTOP polymer optical fibers using phase mask and 400 nm femtosecond laser,” Opt. Express 29(16), 25824–25835 (2021). [CrossRef]  

10. C. Broadway, D. Kinet, A. Theodosiou, K. Kalli, A. Gusarov, C. Caucheteur, and P. Mégret, “CYTOP fibre Bragg grating sensors for harsh radiation environments,” Sensors 19(13), 2853 (2019). [CrossRef]  

11. A. G. Leal-Junior, A. Theodosiou, C. A. R. Diaz, L. M. Avellar, K. Kalli, C. Marques, and A. Frizera, “FPI-POFBG angular movement sensor inscribed in CYTOP fibers with dynamic angle compensator,” IEEE Sens. J. 20(11), 5962–5969 (2020). [CrossRef]  

12. A. G. Leal-Junior, H. R. O. Rocha, A. Theodosiou, A. Frizera, C. Marques, K. Kalli, and M. R. N. Ribeiro, “Optimizing linearity and sensitivity of 3D-printed diaphragms with chirped FBGs in CYTOP fibers,” IEEE Access 8, 31983–31991 (2020). [CrossRef]  

13. A. G. Leal-Junior, A. Theodosiou, R. Min, J. Casas, C. R. Diaz, W. M. D. Santos, M. J. Pontes, A. A. G. Siqueira, C. Marques, K. Kalli, and A. Frizera, “Quasi-distributed torque and displacement sensing on a series elastic actuator’s spring using FBG arrays inscribed in CYTOP fibers,” IEEE Sens. J. 19(11), 4054–4061 (2019). [CrossRef]  

14. A. Theodosiou, X. Hu, C. Caucheteur, and K. Kalli, “Bragg gratings and Fabry-Perot cavities in low-loss multimode CYTOP polymer fiber,” IEEE Photonics Technol. Lett. 30(9), 857–860 (2018). [CrossRef]  

15. A. Theodosiou, A. Lacraz, A. Stassis, C. Koutsides, M. Komodromos, and K. Kalli, “Plane-by-plane femtosecond laser inscription method for single-peak Bragg gratings in multimode CYTOP polymer optical fiber,” J. Lightwave Technol. 35(24), 5404–5410 (2017). [CrossRef]  

16. M. Ochoa, J. F. Algorri, P. Roldán-Varona, L. Rodríguez-Cobo, and J. M. López-Higuera, “Recent advances in biomedical photonic sensors: A focus on optical-fibre-based sensing,” Sensors 21(19), 6469 (2021). [CrossRef]  

17. M. Lomer, A. Quintela, M. López-Amo, J. Zubia, and J. López-Higuera, “A quasi-distributed level sensor based on a bent side-polished plastic optical fibre cable,” Meas. Sci. Technol. 18(7), 2261–2267 (2007). [CrossRef]  

18. C. Schulze, A. Lorenz, D. Flamm, A. Hartung, S. Schröter, H. Bartelt, and M. Duparré, “Mode resolved bend loss in few-mode optical fibers,” Opt. Express 21(3), 3170–3181 (2013). [CrossRef]  

19. D. Pallares-Aldeiturriaga, L. Rodriguez-Cobo, A. Quintela, and J. M. Lopez-Higuera, “Curvature sensor based on In-Fiber Mach–Zehnder interferometer inscribed with femtosecond laser,” J. Lightwave Technol. 35(21), 4624–4628 (2017). [CrossRef]  

20. J. E. Harvey and R. N. Pfisterer, “Understanding diffraction grating behavior: including conical diffraction and Rayleigh anomalies from transmission gratings,” Opt. Eng. 58(08), 1 (2019). [CrossRef]  

21. M. S. Amin, J. W. Yoon, and R. Magnusson, “Optical transmission filters with coexisting guided-mode resonance and Rayleigh anomaly,” Appl. Phys. Lett. 103(13), 131106 (2013). [CrossRef]  

22. J. E. Saunders, C. Sanders, H. Chen, and H.-P. Loock, “Refractive indices of common solvents and solutions at 1550 nm,” Appl. Opt. 55(4), 947–953 (2016). [CrossRef]  

23. G. M. Hale and M. R. Querry, “Optical constants of water in the 200-nm to 200-µm wavelength region,” Appl. Opt. 12(3), 555–563 (1973). [CrossRef]  

24. M. F. H. Arif and M. J. H. Biddut, “Enhancement of relative sensitivity of photonic crystal fiber with high birefringence and low confinement loss,” Optik 131, 697–704 (2017). [CrossRef]  

25. A. Urrutia, I. D. Villar, P. Zubiate, and C. R. Zamarreño, “A comprehensive review of optical fiber refractometers: Toward a standard comparative criterion,” Laser Photonics Rev. 13(11), 1900094 (2019). [CrossRef]  

26. W. Zhou, Y. Zhou, and J. Albert, “A true fiber optic refractometer,” Laser Photonics Rev. 11(1), 1600157 (2017). [CrossRef]  

27. Z. Liu, C. Shen, Y. Xiao, J. Gong, J. Wang, T. Lang, C. Zhao, C. Huang, Y. Jin, X. Dong, Y. Zhang, Z. Jing, W. Peng, and Y. Semenova, “Liquid surface tension and refractive index sensor based on a tilted fiber Bragg grating,” J. Opt. Soc. Am. B 35(6), 1282–1287 (2018). [CrossRef]  

28. X. Hu, C.-F. J. Pun, H.-Y. Tam, P. Mégret, and C. Caucheteur, “Tilted bragg gratings in step-index polymer optical fiber,” Opt. Lett. 39(24), 6835–6838 (2014). [CrossRef]  

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (10)

Fig. 1.
Fig. 1. Diagram of the CYTOP-based optical structure.
Fig. 2.
Fig. 2. Geometric parameters that define the (a) major axis and (b) minor axis of the ellipse generated in the XY plane (Fig. 1).
Fig. 3.
Fig. 3. Schematic of the fibre in the $Z$-plane (no transformation) and in the $W$-plane (conformal mapping applied) for (a) unpolished ($R-$sweep simulation), (b) and polished (at $\Psi =90^\circ$) fibre ($\epsilon -$sweep simulation). Refractive index profiles are depicted for each case (with and without transformation).
Fig. 4.
Fig. 4. RSoft simulation of the resulting field in the fibre cross-section at $\Psi =90^\circ$ for bending radii of (a) $2$, (b) $4$, (c) $5$, (d) $6$, (e) $10$ and (f) $20$ mm. The attenuation obtained at the end of the curvature ($\Psi =180^\circ$) for each radius is $31$, $16.3$ $12.7$, $9.7$, $5.1$ and $2.7$ dB, respectively. The boundaries corresponding to core (dotted line) and cladding (solid line) are also depicted.
Fig. 5.
Fig. 5. (a) RSoft simulation of fibre losses ($L$) due to curvature ($R=5$ mm) as a function of polishing depth ($\epsilon$) and surrounding refractive index. (b) $\partial L/\partial \textrm {SRI}$ of the values represented in (a). The SRI range between $1.318$ and $1.336$ is highlighted, due to its remarkable sensitivity.
Fig. 6.
Fig. 6. Diagram of 1-D diffraction grating. The diffraction orders in transmission and reflection are depicted. It is assumed that there is no cladding (removed in polishing).
Fig. 7.
Fig. 7. Femtosecond laser inscription setup.
Fig. 8.
Fig. 8. (a) $5\times$ microscope image of the ellipse generated as a result of polishing, showing both cladding and core, as well as overcladding. $20\times$ microscope images (b) before (c) and after diffraction grating inscription of the region indicated in (a) are also shown. (d) Image of the complete CYTOP-based optical structure.
Fig. 9.
Fig. 9. (a) Transmission spectra when the structure is immersed in liquids with different RI. The $6$ peaks and valleys on which the tracking in sensing is performed are highlighted. (b) Wavelength shift ($\Delta \lambda$) of each of the $6$ peaks/valleys for the liquids used, with their quadratic fit and error bars ($\textrm {CL}=95\%$). (c) Simulation and experimental transmission losses for $\lambda =1550$ nm, as well as experimental results for $\lambda =1500$ nm and $\lambda =1600$ nm, with the associated error bars ($\textrm {CL}=95\%$)
Fig. 10.
Fig. 10. Wavelength shift ($\Delta \lambda$) of each peak/valley as a function of temperature. The structure has been immersed in water ($n_L=1.3153$), and a temperature sweep between $23$ and $40$ $^\circ$C has been performed. A linear fit is made, with the slope of each fit. It should be noted that the wavelength step in the spectrum is $50$ pm.

Tables (4)

Tables Icon

Table 1. Works related to the manufacture of optical structures in CYTOP fibres by means of femtosecond laser.

Tables Icon

Table 2. Simulation parameters used in RSoft for both cases ( R sweep and ϵ sweep).a

Tables Icon

Table 3. Solute concentration, n D and n ( 1550  nm ) of the liquids used in the experiments.

Tables Icon

Table 4. Parameters of the quadratic and linear fit of the six curves depicted in Fig. 9(b).

Equations (8)

Equations on this page are rendered with MathJax. Learn more.

y 2 + ( R + d o c 2 ϵ ) 2 = ( R + d o c 2 ) 2 ,
x 2 + ( d o c 2 ϵ ) 2 = ( d o c 2 ) 2 .
W = u + i v = f ( Z ) = f ( x + i y ) = R ln ( Z R ) ,
| d Z d W | = exp ( u R ) .
n t i l t ( u , v ) = n r e a l ( u , v ) exp ( u R ) n r e a l ( u , v ) ( 1 + u R ) . d c l a d R
r ( θ ) = R + d c l a d 2 ϵ sin θ , θ [ π 2 α , π 2 + α ] ,
n L sin θ d = n c o r e sin θ i + λ Λ p , p Z .
n c o r e sin θ r = n c o r e sin θ i + λ Λ p , p Z .
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.