s13538-022-01088-6
s13538-022-01088-6
https://doi.org/10.1007/s13538-022-01088-6
Received: 28 September 2021 / Accepted: 9 March 2022 / Published online: 17 March 2022
© The Author(s) under exclusive licence to Sociedade Brasileira de Física 2022
Abstract
This article exhibits a fiber plasmonic sensor based on a silver (Ag) layer and a silicon nitride (Si3N4) layer with a wide
dynamic range of refractive index sensing. The Taguchi (L9) orthogonal array method is applied to optimize the design
parameters such as fiber core, sensing region length, and the thickness of Ag/Si3N4 layers. The performance of the structure
is investigated for the full-width half maxima (FWHM) as the smaller the better (STB). The smaller FWHM favors accurate
detection, high-quality factor, and better sensitive detection of biomolecules. With the use of analysis of variance (ANOVA),
it is evident that for the normalized transmitted power, the Ag layer thickness contributes 47.36%, while the S i3N4 layer
thickness contributes 0.06% only. Furthermore, the length of the sensing region has the highest dominating factor effect of
53.26% whereas the thickness of the Si3N4 layer has the least dominating factor effect of 6.90% on FWHM. This work has
shown the highest sensitivity of 6287 nm/RIU whereas the quality factor and detection accuracy are 873.19 RIU−1 and 87.31
respectively. Hence, Taguchi’s optimization approach is suitable in multilevel optimization of different control factors that
lead to the robust design of the fiber SPR sensor.
Keywords Taguchi optimization · Surface plasmon resonance · Silver · Silicon nitride · Sensor
13
Vol.:(0123456789)
80 Page 2 of 10 Brazilian Journal of Physics (2022) 52: 80
unclad fiber core’s surface. In comparison to these materials, fiber [24]. The RI of the fused silica core (ncore) may be
gold is extremely stable. Because silver and copper have calculated using the given dispersion equation [25] :
weak chemical stability and these are supposed to oxidize )1 ∕2
fast. This makes them lesser stable for the SPR sensor. As
(
a1 𝜆2 a2 𝜆2 a3 𝜆2
a result, applying a protective layer on the surface of these ncore (𝜆) = 1+ + + (1)
𝜆2 − b21 𝜆2 − b22 𝜆2 − b23
layers is a feasible option for making an SPR sensor more
durable [13–15]. In recent times, the use of transition metal Here, a1 = 0.6961663, b1 = 0.0684043, a2 = 0.4079426,
dichalcogenide (TMDC) materials such as WSe2, MoS2, b2 = 0.1162414, a3 = 0.8974794, and b3 = 9.896161 are the
MoSe2, and WS2 is seen in SPR sensing. These materi- Sellmeier coefficients and “λ” is the wavelength in nm [26].
als serve as supporting materials for enhancing plasmons in Ag metal layer of the thickness of 50 nm is grown over the
SPR sensing [16–21]. unclad portion of the fiber. Liquid-phase deposition tech-
In this work, a protective layer of silicon nitride (Si3N4) nique can be used to deposit silver layer over the unclad
has been deposited over the surface of the Ag layer that makes surface of optical fiber [27]. The dielectric constant of this
this structure more stable. Si3N4 has a huge bandgap (~ 5 eV) layer (𝜖 Ag ) is computed using the given Drude relation [28] :
and received interest due to its excellent chemical and thermal
stability [22, 23]. Here, we have used the Taguchi method to 𝜆2 𝜆c
optimize the designing parameters employed in the proposed
𝜖Ag (𝜆) = 1 − (2)
𝜆2p (𝜆c + i𝜆)
fiber SPR sensor. In the traditional approach, one factor at a
time is optimized such that when the thickness of one layer is here, 𝜆c stands for collision (damping) wavelength and is
optimized, the thickness of the other layer is kept fixed. This defined as the wavelength corresponding to the damping of
process makes the high numbers of optimization tests. Due to the oscillations of the electron density because of the colli-
this reason, the parameters generated with this approach are sion among the electrons. The value of 𝜆c is 1.7614 × 10−5 m.
not sufficient for the designing of such sensors. The Taguchi Furthermore, 𝜆p stands for plasma wavelength. It is the
technique aids in the creation of a robust design by taking wavelength that corresponds to the frequency of the elec-
into account all of the sensor’s control elements. As a result, tron density oscillations in the metal. The value of 𝜆p is
a high-performance SPR biosensor can be made with a small 1.4541 × 10−7 m [29]. Next, to protect the silver layer from
number of tests. oxidation, a thin layer of Si3N4 of the thickness of 5 nm is
grown over the covered region of the Ag layer. Magnetron
sputtering technique can be used to deposit S i3N4 film on the
surface of silver layer [30]. The dielectric constant for this
2 Sensor Design Parameters can be computed using following relation [31] :
( )
Taguchi’s optimized design of the proposed fiber-based SPR 3.0249𝜆2 40314𝜆2
𝜖Si3N4 (𝜆) = 1 + +
sensor is represented in Fig. 1. Here, the fused silica fiber 𝜆2 − 0.13534062 𝜆2 − 1239.8422
core having a diameter of 600 µm is utilized to carry light (3)
waves from one side of the fiber to another side, in which At the transmitting end of the fiber, light is incident and
a 10-mm length of fiber cladding is removed. Fiber optic propagates through the fiber. The evanescent wave is gen-
probe with 600 µm of core diameter can be fabricated by erated at the metal surface that excites the SPs wave when
removing 10-mm cladding portion from the middle region the light wave is completely reflected in the detecting area.
of a fiber using a sharp blade. After that, acetone and de- Lastly, we chose water as the detecting medium, which has a
ionized water can be used to clean the unclad area of the RI of 1.33 and measured the change in RI from 1.33 to 1.43
13
Brazilian Journal of Physics (2022) 52: 80 Page 3 of 10 80
13
80 Page 4 of 10 Brazilian Journal of Physics (2022) 52: 80
Table 1 Studied control factors and level representation the help of Fresnel multilayer reflection theory and transfer
Control factors Level 1 Level 2 Level 3
matrix method. Furthermore, Taguchi L9 orthogonal array
is used to calculate SNR for the desired control factors that
Core diameter (µm) 400 500 600 we can consider smaller the better for FWHM.
Length of sensing region (mm) 10 15 20 Mean and SNR are calculated for FWHM and normal-
Thickness of silver layer (nm) 40 50 60 ized transmitted power for all nine experiments. In Taguchi
Thickness of Si3N4 layer (nm) 3 5 7 design analysis, SNR can be categorized as larger the better,
nominal the best, and smaller the best. In designing of fiber
SPR sensor, small FWHM and large normalized transmitted
structure. These control factors are listed in Table 1. With these power are the prerequisites. Normalized transmitted power
control factors, the appropriate orthogonal array is L9. FWHM with SNR (smaller the better) for different experiments is
and normalized transmitted power for all nine experiments can shown in Table 4. Here, SNR is selected for smaller the bet-
be calculated from Fig. 2. With the Taguchi array structure, ter for best performance of the proposed sensor that can be
the number of simulation trials has been cut from 81 (34) to 9 calculated by the following equation [40] :
only. This saves cost, time, and the number of attempts while
1 ∑n 2
optimizing the SPR sensor’s design parameters. SNRSmaller the better = −10 log10 ( Y ) (12)
n i=1 i
FWHM is the important parameter in designing of sen-
sor that mainly gives the spectral width of the SPR curve. If here, n represents the number of trials and Yi denotes the
the FWHM of the SPR curve is large, the quality factor and response for a specific factor-level combination.
detection accuracy decrease as expressed in Eqs. 10 and 11
respectively. Here, we have optimized the performance for 5.1 Interpretation of Response Table for FWHM
smaller the better SNR. Furthermore, SNR for the FWHM
with the Si3N4 layer and without the Si3N4 layer is shown The response table is utilized to determine the appropriate
in Tables 2 and 3 respectively. It is seen that adding the level for each control factor. Response characteristics of
Si3N4 layer reduces the FWHM which is the primary pre- different control factors with SNR for FWHM are depicted
requisite for designing for such an SPR sensor. In Table 2, in Table 5. In this, the factors with the greatest effect on
FWHM is calculated from the SPR curve that can plot with FWHM are identified using rank values and delta. The rank
13
Brazilian Journal of Physics (2022) 52: 80 Page 5 of 10 80
Table 2 FWHM with their SNR Experi- Control factors FWHM (nm) SNR
for different experiments ment
no. Core Length of Thickness of Thickness of
diameter sensing region silver layer Si3N4 layer
(µm) (mm) (nm) (nm)
Table 3 FWHM with their Experiment Control factors FWHM (nm) SNR
SNR without the Si3N4 layer for no.
different experiment Core Length of sensing Thickness of silver
diameter (µm) region (mm) layer (nm)
Table 4 Normalized transmitted power with their SNR for different experiments
Experiment Control factors Normalized SNR (smaller
no. transmitted power the better)
Core diameter Length of sensing Thickness of silver Thickness of Si3N4 (STB)
(µm) region (mm) layer (nm) layer (nm)
of sensing region length is one while the rank of silver 5.2 Interpretation of Response Table
layer thickness, core diameter, and Si3N4 layer is 2, 3, and for Normalized Transmitted Power
4 respectively. This means that the length of the sensing
region has the highest contributing effect of smaller SNR The response table for different control factors for the nor-
for FWHM. The overall mean value of SNR for FWHM malized transmitted power is shown in Table 6. In this, the
from Fig. 3 is obtained at − 23.365. thickness of the Ag layer with a delta value of 19.80 and
13
80 Page 6 of 10 Brazilian Journal of Physics (2022) 52: 80
Table 5 Response table for SNR for different control factors for Table 6 Response table for different control factors for normalized
FWHM; smaller the better transmitted power; smaller the better
Level Core Length of Thickness of Thickness of Level Core Length of Thickness of Thickness of
diameter sensing region silver layer Si3N4 layer diameter sensing region silver layer Si3N4 layer
(µm) (mm) (nm) (nm) (µm) (mm) (nm) (nm)
13
Brazilian Journal of Physics (2022) 52: 80 Page 7 of 10 80
Table 7 Analysis of variance Control factors Degree of Sum of square Mean square Contribution
(ANOVA) for FWHM with freedom
different control factors
Core diameter (µm) 2 32.54 16.27 10.00%
Length of sensing region (mm) 2 173.34 86.67 53.27%
Thickness of silver layer (nm) 2 97.08 48.54 29.83%
Thickness of Si3N4 layer (nm) 2 22.47 11.23 6.90%
Table 8 Analysis of variance Control factors Degree of Sum of square Mean square Contribution
(ANOVA) for normalized freedom
transmitted power with different
control factors Core diameter (µm) 2 0.034161 0.017080 21.00%
Length of sensing region (mm) 2 0.051377 0.025688 31.58%
Thickness of silver layer (nm) 2 0.077075 0.038537 47.36%
Thickness of Si3N4 layer (nm) 2 0.000093 0.000047 0.06%
parameters, we have a designed a fiber SPR sensor for the RI this, when we are increasing the RI of the sensing medium,
ranging from 1.33 to 1.43. The plot for normalized transmit- the resonance wavelength shifts to a higher value. The
ted power with RI of sensing medium is shown in Fig. 5a. In value of resonance wavelength for respective RI of sensing
medium is plotted in Fig. 5b and also depicted in Table 10.
The values of sensitivity, quality factor, and detection accu-
Table 9 Best selection of control parameter with Taguchi method racy which are given in Table 10 can be calculated by using
Eqs. 9, 10, and 11 respectively. Here, the maximum sensi-
S. no. Control parameter Parameter
value tivity of 6287 nm/RIU is achieved for the varying RI from
1.33 to 1.43 while the maximum quality factor and detection
1 Core diameter (µm) 600 accuracy are 873.19 RIU−1 and 87.31 respectively. As a
2 Length of sensing region (mm) 10 result, the Taguchi approach can be utilized to optimize the
3 Thickness of silver layer (nm) 50 designing parameters and can assist in designing the fiber
4 Thickness of Si3N4 layer (nm) 5 SPR sensor.
13
80 Page 8 of 10 Brazilian Journal of Physics (2022) 52: 80
Fig. 5 (a) Plot of the SPR spectra and (b) the shift of resonance wavelength for the core-Ag-Si3N4 configuration with different RI of the sensing
medium
13
Brazilian Journal of Physics (2022) 52: 80 Page 9 of 10 80
13
80 Page 10 of 10 Brazilian Journal of Physics (2022) 52: 80
35. N. Mudgal, A. Saharia, A. Agarwal, J. Ali, P. Yupapin, G. Singh, 39. N.R. Mohamad, M.F.M.R. Wee, M.A. Mohamed, A.A. Hamzah,
Modeling of highly sensitive surface plasmon resonance (SPR) P.S. Menon, Multi-response optimization of chromium/gold-
sensor for urine glucose detection. Opt. Quant. Electron. 52, 307 based nanofilm Kretschmann-based surface plasmon resonance
(2020) glucose sensor using finite-difference time-domain and Taguchi
36. L. Wu, J. Guo, Q. Wang, S. Lu, X. Dai, Y. Xiang, D. Fan, Sensitiv- method. Nanomater. Nanotechnol. (2020). https://d oi.o rg/1 0.1 177/
ity enhancement by using few-layer black phosphorus-graphene/ 1847980420982119
TMDCs heterostructure in surface plasmon resonance biochemi- 40. S. Sahu, J. Ali, P. Yupapin, G. Singh, Effectiveness of Taguchi
cal sensor. Sens. Actuat. B Chem. 249, 542–548 (2017) method for the optimization of narrowband optical filters based
37. Y. Saad, M. Selmi, M.H. Gazzah, A. Bajahzar, H. Belmabrouk, on grating waveguides. MicrosystTechnol 25, 789–795 (2019)
Performance enhancement of a copper-based optical fiber SPR
sensor by the addition of an oxide layer. Optik 190, 1–9 (2019) Publisher's Note Springer Nature remains neutral with regard to
38. A. Sharma, B.D. Gupta, Fiber optic sensor based on surface plas- jurisdictional claims in published maps and institutional affiliations.
mon resonance with Ag–Au alloy nanoparticule films. Photon.
Nanostruct.-Fundam. Appl. 3, 30–37 (2005)
13