Mid-infrared water pollutant sensor based on SPR-PCF

In this paper, a highly sensitive water pollutant optical sensor is proposed and analyzed. The suggested sensor consists of photonic crystal fiber with a core surrounded by four elliptical holes infiltrated with the studied analyte (pure/polluted water sample). In addition, two gold nanorods are mounted horizontally at the inner surfaces of two horizontal elliptical holes. The proposed sensor can efficiently detect dissolved pollutants in water such as nitric acid (HNO3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathrm{HNO}}_{3}$$\end{document}) with concentrations of 14, 23 and 35% and H2O2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathrm{H}}_{2}{\mathrm{O}}_{2}$$\end{document} with concentrations of 7, 15 and 30%. The dissolved pollutants in water (analyte) have refractive indices (RIs) in the range of 1.350–1.355 in the mid infrared regime from λ = 2200 nm to λ = 3500 nm. A fully vectorial finite element method (FVFEM) is employed for the modal analysis of the reported structure. The geometrical parameters are studied to maximize the RI sensitivity where a high sensitivity of 36,000 nm/RIU is achieved between the studied RIs of 1.350 and 1.355. The obtained RI sensitivity is higher than those of the recent reported sensors in the literature especially those operating in the mid infrared wavelengths.


Introduction
Water is the most crucial commodity for life.Therefore, the consumption of water by humans should be safe, easily accessible, adequate, and free from any kind of contamination (Madhav et al. 2020).Pollutants in water bodies pose a severe threat to human health as well as aquatic ecosystem.Thus, the development of simple and accurate devices for the detection of water pollutants is necessary to have a healthy society.Optical sensors, via which the polluted water samples can be distinguished from the normal water samples, have their uniqueness in water quality assessment due to their compact size, high sensitivity, and immunity to stray fields.One of the major optical sensing mechanisms is surface plasmon resonance (SPR).Such technique relies on the coupling between surface plasmon modes (SPMs) constructed at the metal/dielectric interfaces and the guided core modes (Younis et al. 2017(Younis et al. , 2018)).At resonance, the effective indices of the SPM and core guided mode are equal.In this situation, phase matching occurs where maximum power transfer Extended author information available on the last page of the article 966 Page 2 of 18 takes place from the guided core mode to the SPM.A shift in the resonance wavelength is detected when the analyte refractive index is changed, and thus the sensitivity can be measured.Many forms of optical devices including optical absorbers (Wang et al. 2020b), sensors (Hameed et al. 2017;Abdelghaffar et al. 2022), filters (Gamal et al. 2021a), polarization rotators (Hameed et al. 2010(Hameed et al. , 2013b;;Hameed and Obayya 2011a, b), polarization converters (Hameed et al. 2011a(Hameed et al. , 2013a)), and optical instruments (Hameed et al. 2009a(Hameed et al. , b, 2011b) ) have been also designed based on the SPR mechanism.The general idea of an optical sensor is to track the sensitivity of electromagnetic coupling phenomenon to the change in the refractive index of the analyte.
In recent times, plasmonic photonic crystal fibers (PCFs) based sensors have gathered a considerable attention due to its exceptional characteristics (Abdelghaffar et al. 2023).PCFs have many advantages including large degrees of design freedom due to the flexibility of air holes arrangement and diversity of filling materials (Poli et al. 2007).Moreover, PCFs have many outstanding properties in comparison to the standard optical fibers, such as high birefringence, endless single-mode propagation over a wide wavelength range, low-loss guidance, high nonlinearity, transfer of high power and low bending loss (De et al. 2019).Therefore, PCFs have been widely employed in refractive index sensing applications.In this context, Dash and Jha (2016) have reported a simple Indium Tin Oxide (ITO) based D-shaped PCF sensor.Their sensor showed a relatively low wavelength sensitivity of 5200 nm/RIU and a resolution of 1.92 × 10 −5 RIU in the wavelength range from 1600 to 2100 nm (Dash and Jha 2016).Another D-shaped PCF sensor with gold as a plasmonic material has been numerically investigated in Wu et al. (2018).A maximum sensitivity of 31,000 nm/RIU has been obtained in the near-infrared band (Wu et al. 2018).Moreover, Chen et al. (2018) have demonstrated a different structure of D-shaped PCF-SPR sensor for low RI detection where a maximum spectral sensitivity of 11,055 nm/RIU and RI resolution of 9.05 × 10 −6 RIU have been achieved in the wavelength range from 2.35 to 2.65 µm.Additionally, Gao et al. (2014) have presented an optical fiber based SPR sensor with bi-layered plasmonic films of gold and TiO 2 .The RI resolution has been estimated to be 2.7 × 10 −5 RIU and the resonance wavelength has been extended to the near-infrared region.Further, an extra-broad PCF RI sensor has been investigated based on a triangular lattice PCF and four large-size channels (An et al. 2017b).The detection range of RI was 1.30-1.79and the resonance wavelength has been increased from 1680 to 1870 nm.Turduev et al. (2017) have presented a Mid IR T-shaped photonic crystal waveguide for RI sensing with sensitivities of 1040 nm/RIU, 500 nm/RIU, 530 nm/ RIU and 390 nm/RIU within the analyte RI ranges of 1. 05-1.10, 1.10-1.30, 1.10-1.25 and 1.00-1.30, respectively. Additionally, An et al. (2017a) have introduced a Mid IR RI sensor based on D-shaped PCF to achieve a sensitivity of 10,493 nm/RIU with a very high resolution of 9.53 × 10 −6 RIU.Wang et al. (2020a) have presented a near IR dual core SPR-PCF sensor based on dual gold nanowires with sensitivity of 17,500 nm/RIU.Further, Gayraud et al. (2007) have demonstrated a methane gas sensor based on hollow core PCF.Such sensor has been utilized to measure the transmission spectra through the methane absorption region of 3.15-3.35µm (Gayraud et al. 2007).Another methane sensor based on Fourier transform infrared spectroscopy has been developed by Gayraud et al. (2008).This sensor (Gayraud et al. 2008) is fabricated to measure qualitative methane concentrations to 1000 ppm (parts in 10 6 ).Furthermore, Liu et al. (2017) have reported a low RI sensor with two open-ring chan- nels based on SPR-PCF.The sensor achieves sensitivity of 5500 nm/RIU with maximum resolution of 7.69 × 10 −6 RIU in the Mid IR region between 2550 and 2900 nm.Recently, Gamal et al. (2022b) have introduced a multifunctional RI and temperature sensor with sensitivities of 34,600 nm/RIU and 21.05 nm/°C, respectively, in the near IR region.Moreover, Gamal et al. (2022a) have reported an environmental and low RI sensor in near IR with sensitivity of 19,100 nm/RIU.Herein, a highly sensitive plasmonic PCF based water pollutant sensor is introduced.The suggested PCF structure consists of two successive rings of air holes.The first ring of air holes is composed of four ellipses arranged in a square lattice and infiltrated with the targeted water samples.The outer ring consists of twenty circular air holes organized in a circular lattice to create the cladding region in a silica (SiO 2 ) background.In addition, two gold nanorods are attached to the inner surfaces of the two horizontal ellipses to allow the SP coupling with the core modes.COMSOL Multiphysics software package (COMSOL 5.1 2020) based on finite element method (FEM) (Obayya et al. 2002) is employed to perform the modal analysis of the reported sensor.The suggested optical sensor is designed and optimized to efficiently detect the dissolved pollutants in water such as nitric acid (HNO 3 ) with concentrations of 14, 23 and 35% and H 2 O 2 with concentrations of 7, 15 and 30%.In this regard, a very high RI sensitivity of 36,000 nm/RIU is obtained which is much higher than those reported by similar devices in the literature (De et al. 2019;Dash and Jha 2016;Wu et al. 2018;Chen et al. 2018;Gao et al. 2014;An et al. 2017b).Moreover, the structure is shown to be robust against fabrication tolerances.
The rest of the manuscript is organized into 4 subsections.Section 2 presents the design considerations, the utilized materials, and structural parameters of the proposed device.In Sect.3, the simulation results are shown and discussed.In addition, the optimization of the geometrical parameters to sustain sensitivity at maximum value is also proposed.Further, Sect. 3 also includes the study of different pollutants with different concentrations in water with the obtained sensitivities.The conclusions and summary of the proposed results are drawn in Sect. 4.

Design considerations
The two-dimensional cross-section of the suggested PCF water pollutants sensor is depicted in Fig. 1a.The proposed PCF has 4 elliptical holes; two large elliptical holes with radii of a x = 3 µm and a y = 5 µm, and two smaller ones with radii of b x = 2.1 µm and b y = 3 µm as may be seen in Fig. 1a.In addition, an outer ring of air holes arranged in a circular lattice with a hole diameter (d) = 1.1 μm and a hole pitch (Λ) = 2.2 μm is added to construct the cladding region.Moreover, two gold nano-rods are attached to the inner surfaces of the two large elliptical holes with diameter (d g ) of 300 nm as depicted in Fig. 1a.It is worth noting that the four large elliptical holes are infiltrated with the studied water sample.In the same context, the outer layer of air holes is used to reduce the cladding average refractive index.Thus, the index contrast between the solid core region and the cladding increases which ensures a good confinement of the core mode in the core region.Therefore, the light matter interaction with the analyte sample and hence sensor sensitivity will be increased.
The sensing system can be built by coupling the laser source with the PCF sensor via single mode fiber (SMF) (Leon-Saval et al. 2005) as depicted in Fig. 1b.A standard mid infrared single mode fiber (SMF) (Thorlabs) with core RI of 1.49229, cladding RI of 1.47, core diameter of 9 μm, and cladding diameter of 125 μm can be used for the excitation of the proposed PCF.To ensure the coupling between the fundamental modes of the SMF and the suggested PCF, the overlap integral and power coupling coefficient between them are investigated.The numerical results show that the overlap integral is close to unity (~ 0.99).So, it can be concluded that the coupling losses and reflections at the junctions between the SMF and the PCF have relatively little effect on the efficiency of the proposed structure.The transmitted light will be coupled again to the spectrometric system via SMF to measure the transmission and hence the losses through the proposed design.Due to the plasmonic effect, maximum loss occurs at a resonance wavelength depending on the RI of the water sample.The shift in the resonance wavelength can be used to detect the RI variation.
In this work, the suggested PCF is analyzed using the full vectorial finite element method (FVFEM) (Koshiba and Tsuji 2000) via COMSOL Multiphysics software package (COMSOL 5.1 2020) with computational domain of 328.6 µm 2 where the number of degrees of freedom solved for is 64,269 with a maximum element size of 0.5 µm 2 .The cross section of the proposed structure is discretized into very small triangular elements with nonuniform meshing capabilities.In the simulation model, the outer coating (polymer or plastic) of the proposed optical fiber is replaced by a perfectly matched layer (PML) (Koshiba and Tsuji 2000) is applied to the transverse directions to truncate the simulation domain.Initially, extensive convergence tests are performed to ensure the numerical accuracy and stability of the results for the proposed sensor.The numerical results show that convergence occurs when the PML thickness is ≥ 0.3 µm and the maximum element size is ≤ 3.0 µm is used.Therefore, a PML with thickness of 2.0 µm is applied to ensure a maximum element size of 0.5 µm.However, an element size of 1 nm is utilized in the thin layers, i.e., Au and silica layers.
Different modes supported by the suggested structure are characterized by two main parameters: mode effective index ( n ef f ) and confinement loss (L).In this study, L is calculated from the imaginary part of the mode effective index according to Eq. ( 3) (Crutchfield et al. 2004): The sensitivity of the proposed SPR sensor can be calculated for the RI range 1.35-1.355using the wavelength interrogation method depicted in Eq. ( 4) (Rifat et al. 2015): where Δλ r is the shift occurs to the resonance wavelength due to the polluted water sample and Δn is the analyte RI change that is equal to 0.001.The resolution is another critical parameter for sensing process as it represents the smallest change that can be detected in the RI and is determined through the following equation (Gamal et al. 2021b): where the wavelength resolution is represented by Δλ min .In this work, it is assumed that Δλ min is 0.1 nm as suggested in Wang et al. (2020a).Moreover, the figure of merit (FOM) parameter is also studied for the proposed sensor to show the accuracy and high degree of selectivity.In this investigation, FOM is defined as the ratio of the sensitivity to the full width half maximum as depicted in Eq. ( 6) (Gamal et al. 2022b): where FWHM is the full width at half-maximum of the corresponding resonance peak. (3) 966 Page 6 of 18 Figure 2 depicts the wavelength dependent refractive index (n) and the extinction coefficient (k) of HNO 3 and H 2 O 2 with different concentrations in water.It is revealed from Fig. 2 that there are slight changes in both n and k of the two materials in water by changing the concentration.This may reflect on the behavior of the proposed sensor.However, comparing each case with pure water is the most important thing in calculating the sensitivity of the proposed sensor.In addition, this small difference in refractive indices may be compensated by the structure of the sensor in order to keep the sensitivity as high as possible.In the proposed structure, the wavelength sensitive coupling process between the fundamental mode and SPM is significantly affected by very small changes in the analyte RI.Thus, the shift in λ r would be large and the sensitivity becomes sufficient.

Simulation results
The calculated effective indices of the studied x-and y-polarized core modes supported by the PCF core are 1.3676 + 7.0269 × 10 −4 i (Fig. 3a, x-polarized core mode), 1.3671 + 7.0293 × 10 −4 i (Fig. 3b, x-polarized SP mode), and 1.3664 + 1.23899 × 10 −7 i (Fig. 3c, y-polarized core mode) at λ = 3116 nm and analyte RI of 1.355.It is worth noting that the imaginary part of the effective index of the x-polarized core mode is large due to the strong coupling with the x-polarized SP mode shown in Fig. 3b where both modes have approximately the same field distribution and effective indices.However, the y-polarized mode (Fig. 3c) is well confined in the core region and is not coupled to the SP modes around the plasmonic materials.This is evident by the very small imaginary part of the effective index of the y-polarized core mode (Fig. 3c).
The dispersion characteristics and the loss behavior of the x-polarized and y-polarized core modes and SPM for n a = 1.355 are depicted in Fig. 4. The real parts of the n ef f of the x-polarized core mode and the corresponding SPM are shown as red and blue lines, respectively.It may be seen that throughout the whole studied wavelength range, R e (n ef f ) of both x-polarized mode and SPM decreases with the increase of the wavelength.The confinement loss of the x-polarized and y-polarized core modes are represented by the black and green lines, respectively.The confinement loss of the x-polarized core mode has a peak at λ = 3.116 μm, while the confinement loss of the y-polarized core mode (represented by the green line) is very small throughout the whole wavelength range with no loss peaks as shown in Fig. 4. The insets (a), (b) and (c) in Fig. 4 show the electric field norm of the SPM, at the phase matching point and the x-polarized core mode, respectively.At the resonance wavelength (3.116 μm), there is a strong coupling between the SPM and the x-polarized core mode where these two modes have approximately the same electric field distribution with nearly identical effective indices.This results in a very strong coupling and accordingly a peak in the loss spectrum x-polarized mode at λ = 3.116 μm as may be seen in Fig. 4.
According to the world health organization (WHO) statistics (WHO 2022), the accepted limits of the pollutants in water are in the range of milligram per 1 l of water.Therefore, these small, dissolved amounts will not have a large effect on the RI change of the water and the RI of the pure and polluted water are very close to each other.Accordingly, our analysis of water sample RI depends on a very small RI variation of 0.001 RIU throughout the rest of the work.This small variation will be taken above the MIR average water RI value (n = 1.35) according to Hale and Querry (1973).
Initially, the geometrical parameters of the proposed PCF structure are taken as a x = 3 µm, a y = 5 µm, b x = 2.1 µm, b y = 3 µm, d g = 300 nm, and d = 1.1 µm at which the calculated sensitivity is equal to 11,100 nm/RIU.Figure 5 shows the effective indices of the coupled x-polarized core and SP modes and the loss spectra of the x-polarized mode at the two analyte RIs of 1.35 and 1.355.It may be seen that by increasing the analyte refractive index, the index contrast between the core and the analyte region decreases.Thus, strong leakage from the core region to the analyte occurs leading to a strong interaction between the core mode and the SPM (that has larger n eff in this case).Due to the increase in n eff of the SPM, the coupling will occur at a shorter wavelength.Thus the variation in the analyte RI induces a shift in the resonance wavelength of 55.5 nm.It can also be observed that, for the two analyte samples, there is a strong coupling between the core and SP modes with the same effective indices and the phase matching occurs at the same wavelength of the confinement loss's peak.
The geometrical parameters of the proposed sensor are the studied to maximize the sensor sensitivity at n a = 1.35.It is expected that, the radius of gold nano-rods, (d g ) would be a very important parameter in controlling the sensitivity of the SPR proposed sensor.To prove this concept, a separate study on the effect of d g on the sensing performance is performed and the results are depicted in Fig. 6. Figure 6a shows the loss spectra of the x-polarized core mode at different values of d g .In this study, the parameters, a x , a y , b x , b y , Ʌ, and d are fixed to 3 µm, 5 µm, 2.1 µm, 3 µm, 2.2 µm and 1.1 µm, respectively.As may be noted from Fig. 6a, the loss peak has a blue shift from λ = 2807 nm to λ = 3150 nm by increasing d g from 300 to 340 nm.The radius of the metallic nano-rods influences the decay rate of the evanescent wave out from the core region and accordingly affects the shape of the SPM.As may be seen from Fig. 6a, thicker gold nano-rod is responsible for increasing the sensitivity, decreasing the loss value, and increasing the FWHM.With continuous increase in the gold nano-rod diameter, the energy of evanescent field decreases, and the confinement loss decreases dramatically.The RI sensitivity as a function of d g is also studied when the analyte RI is changed to n a = 1.355 as shown in Fig. 6b.It is evident from this figure that the RI sensitivity reaches its maximum of 33,200 nm/RIU at d g = 330 nm which will be used in the subsequent simulations.
As the guiding core is a small region that is surrounded by four nonidentical elliptical holes, it is expected that the geometrical parameters of the elliptical holes i.e., a x , a y , b x , and b y would influence the guided modes and consequently the sensitivity of the proposed water pollutant sensor.The variation of a x , a y , b x , and b y is classified into four cases as depicted in Table 1.The effect of these dimensions on the attenuation loss of the x-polarized core mode is studied and the obtained results are depicted in Fig. 7.In this investigation, d and d g are fixed to 1.1 µm and 330 nm, respectively, where the analyte RI is equal to 1.355.It may be seen from Fig. 7a that enlarging the two elliptical holes (i.e., increasing a x , a y , b x , and b y ) by 0.1 µm (Case 2), λ r increases from 2864 to 3234 nm.It is worth men- tioning that as the size of the two elliptical holes increases, the index contrast between the core and cladding regions increases.Thus, the core mode will be more confined in the core region with higher n ef f value.Further, increasing the operating wavelength decreases the n ef f of the core mode to approach the n ef f of the SPM.Thus, the phase matching between the two modes i.e., core mode and SPM can be obtained at longer wavelengths.The RI sensitivity as a function of a x , a y , b x , and b y is also calculated and depicted in Fig. 7b.It may be seen from this figure that the sensitivity has a peak value at a x , a y , b x , and b y equals 2.9 µm, 4.9 µm, 1.7 µm and 2.9 µm, respectively (Case 3).Therefore, the diameters of the two elliptical air holes ( a x , a y , b x , and b y ) are set to 2.9 µm, 4.9 µm, 1.7 µm and 2.9 µm, respectively, to keep the RI sensitivity at its highest value of 36,000 nm/RIU.Since the outer ring of air holes is positioned relatively far away from the core and analyte regions, the variation in its dimensions has a negligible effect on the optical characteristics of the core and SP modes.Thus, d is kept at its original value of 1.1 µm throughout the subsequent simulations.
The variations in the analyte RI have potential applications such as water pollutant detection therefore the RI sensitivity of the proposed structure is also studied.In this investigation, the RI sensitivity of the suggested sensor is analyzed at the optimized geometrical parameters; a x = 2.8 µm, a y = 4.9 µm, b x = 1.7 µm, b y = 2.9 µm, d g = 330 nm, and d = 1.1 µm as in Table 2.
Figure 8a presents the wavelength-dependent loss spectra at various analyte RIs of 1. 350, 1.351, 1.352, 1.353, 1.354 and 1.355 with very small RI variation of 0.001 RIU.Therefore, the proposed sensor can detect dissolved pollutants in water like nitric acid ( HNO 3 ) with percentages of 14, 23 and 35% in water and H 2 O 2 with percentages of 7, 15 and 30% in water and corresponding RIs that overlap with the studied refractive index range (1.350-1.355)(Myers et al. 2018).It is evident from Fig. 8a that generally the resonance wavelength decreases by increasing the analyte RI. Figure 8b shows the resonance wavelength as function of the analyte RI and the fitted line represents the linear fitting over the studied RI range from n = 1.350 to 1.355.The fitted line has a correlation coefficient (r) and coefficient of determination of 0.9999 and 0.9998, respectively, while the sensitivity slope is equal to 36,000 nm/RIU with 36 nm wavelength difference for analyte RI change of 0.001.Therefore, the suggested sensor has a linear-like behavior with high sensitivity through the RI range from 1.35 to 1.355 and is suitable for water pollutant screening.
Figure 9 depicts the FOM of the proposed sensor calculated in the RI range from 1.35 to 1.355.Best FOM value of 250 RIU −1 is obtained at RI of 1.355 .It can be noticed from Fig. 9 that the FOM increases with increasing the RI.Additionally, the obtained FOM in  In order to test the performance of the suggested sensor in detecting different types of pollutants dissolved in water, a study is performed to investigate the effect of adding different concentrations of HNO 3 and H 2 O 2 (Hale and Querry 1973) in water.In this study, the geometrical parameters are fixed to their optimal values based on the previous parametric analysis.Figure 10a, b show the wavelength dependent confinement loss of the core mode for pure water with different concentrations of HNO 3 and H 2 O 2 in water,  HNO 3 (35%) and b pure water, H 2 O 2 (7%), H 2 O 2 (15%) and H 2 O 2 (30%) 966 Page 12 of 18 respectively.It can be noticed from Fig. 10a that there is a large resonance shift between the case of pure water and that of HNO 3 (14%), HNO 3 (23%) and HNO 3 (35%) in water.In addition, Fig. 10b shows that there is also a resonance shift exists between the case of pure water and that of H 2 O 2 (7%), H 2 O 2 (15%) and H 2 O 2 (30%) in water.Table 3 summarizes the resonance wavelengths and the difference in Δλ r for each concentration.From this table, it is revealed that the proposed sensor can be used to detect different pollutants in water with relatively high sensitivity.
Table 4 summarizes a comparison between the proposed sensor with similar devices that have been reported in literature.It is evident that the resolution and RI sensitivity of the proposed water pollutants sensor are better than those introduced with similar devices reported in Dash and Jha (2016), Wu et al. (2018), Chen et al. (2018), Gao et al. (2014), An et al. (2017a, b), Turduev et al. (2017), Wang et al. (2020a), Gayraud et al. (2007Gayraud et al. ( , 2008)), Liu et al. (2017), Gamal et al. (2022a, b).
The suggested device can be fabricated using the well-known stack and draw method (Russell 2003;Pysz et al. 2014), while the elliptical air holes can be fabricated in different ways.Issa et al. successfully fabricated a consistently oriented elliptical-hole fiber (Issa et al. 2004).Guiyao et al. fabricated an SF6 glass made of PCF using the die casting and suggested that this method can also be used for making various PCF structures, including elliptical holes (Guiyao et al. 2006).Zhang et al. reported the fabrication of a methyl methacrylate-based elliptical-hole PCF using the in-situ chemical polymerization method (Zhang et al. 2006).Zhang et al. reported the fabrication of a PCF with elliptical holes using extrusion, canning, and drawing fabrication methods (Zhang et al. 2011).Further, 3D printing technology can be also an effective method for fabricating elliptical-hole PCFs (Ebendorff-Heidepriem et al. 2014).Then, two gold nanorods can be grown and self-aligned to the central air hole's inner surface using stepwise seeded growth method as reported in Feng et al. (2019).In addition, silicone glue can be utilized to cover the cladding circular air holes in order to selectively fill the four elliptical air holes with the analyte (Huang et al. 2004).The suggested PCF is then immersed in the analyte solution, and the fluid will fill the uncovered four elliptical air holes due to the capillary forces that keep the analyte inside them and avoid its drop out.To push out the analyte form the infiltrated holes, a pressure is applied on the syringe pump (Huang et al. 2004).Furthermore, the PCF is spliced to a standard mid infrared SMF (Thorlabs 2023) to ensure strong coupling to the proposed design.Thus, a laser source (Leon-Saval et al. 2005) is used to launch light to the SMF directly.The suggested design has a loss peak of 14,379 dB/cm at the resonance wavelength of 2.795 μm.The specified device length can also be adjusted to meet the requirement of acceptable total leakage loss.For example, a device length of 1 µm can be used to have a total loss of less than 1.5 dB.Therefore, it is believed that the reported device can be experimentally fulfilled.

Fabrication tolerance study
The fabrication tolerance of the proposed water pollutant sensor is also investigated.As of current fabrication technology of silica-based optical fibers, the fabrication tolerance of ± 5% is considered for the proposed designed parameters (Chang et al. 2019).It is found that the diameter of small circular air holes (d) has no effect on the loss spectra of the guided x-polarized mode.However, d g , a x , a y , b x , and b y are considered in this study and the corresponding sensitivities are summarized in Table 5.It may be seen that the sensitivity is better than 32,000 nm/RIU with a variation of ± 5% in the considered geometrical parameters that may occur during the fabrication process.Thus, the tolerance study shows that with the available technology (Chang et al. 2019), the proposed sensor can still be resistant to fabrication inconsistencies.

Conclusion
A highly sensitive RI water pollutant sensor over the range of (1.35:1.355)based on PCF with SPR mechanism is reported and analyzed.The suggested device can detect pollutants like nitric acid ( HNO 3 ) with percentages of 14, 23 and 35% in water and H 2 O 2 with per- centages of 7, 15 and 30% in water which have RIs within the sensing range (1.35-1.355).
The proposed sensor has four elliptical holes infiltrated with the analyte with two gold nanorods attached to the inner surface of the two horizontal elliptical holes.An ultra-high RI sensitivity of 36,000 nm/RIU over the range of (1.35:1.355) is achieved.Therefore, the sensitivity of the suggested sensor surpasses all the recent designs operating in the same ranges.In addition, the reported sensor has advantages in terms of linear sensitivity, simplicity of fabrication, and high robustness to fabrication errors.

Fig. 1 a
Fig. 1 a 2D cross section view of the proposed water sensor, and b The schematic diagram of the proposed sensing system

Fig. 2 Fig. 3
Fig. 2 a, b Wavelength dependent refractive index and extinction coefficient of HNO 3 , respectively, c, d wavelength dependent refractive index and extinction coefficient of H 2 O 2 , respectively with different concentrations in water(Hale and Querry 1973)

Fig. 4
Fig. 4 Wavelength dependent n ef f of the x-polarized mode core (red) X-polarized SPM (blue), y-polarized core mode (green), and the confinement losses (black) at n a = 1.355.(Color figure online)

Fig. 5 Fig. 6 a
Fig. 5 Corresponding effective mode index and loss spectra of the x-polarized core mode at analyte refractive indices of 1.35 and 1.355

Fig. 8 a
Fig. 8 a Confinement loss spectra of the multi-functional sensor for analyte RI of 1.35, 1.351, 1.352, 1.353, 1.354 and 1.355, and b the effect of the analyte refractive index on the resonance wavelength

Fig. 9
Fig.9The FOM of the suggested water pollutant sensor as a function analyte RI

Table 1
Elliptical air holes diameters of the PCF structure depicted in Fig.1a Fig. 7 a Attenuation loss spectra as a function of ellipse dimensions for the four cases depicted in Table 1, where d g is taken as 330 nm, and b the RI sensitivity of the studied four cases 966 Page 10 of 18

Table 2
The optimized geometrical parameters of the proposed structure

Table 3
The resonance wavelengths and difference shift for different tested pollutants

Table 5
Fabrication tolerance study of the different geometrical parameters of the proposed sensor: diameter of gold nanorods d G and ellipses dimensions a X , a Y , b X and b YEthical approvalThe authors declare that there are no conflicts of interest related to this article.Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material.If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.To view a copy of this licence, visit http:// creat iveco mmons.org/ licen ses/ by/4.0/.