Introduction
In the past few decades, the development of surface plasmon based sensing (SPR) techniques has been emphasized as an effective optical detection phenomenon owing to high sensitivity, label free, and real time monitoring [
1]. Particularly, SPR is highly attractive to refractive index (RI) deviations, which performs as a suitable sensor for measuring the chemical, biological, and physical parameters [
2,
3]. The working principle of SPR sensing method can be well explained among the interface of metal and dielectric materials when the light signals and the free electrons jointly overlap at the same frequency [
4]. Mostly commonly functional SPR-sensing platform rely on prism structure and optical fiber technology. However, prism structure demands mechanical and optical modules, these devices have massive arrangements that cause high cost and are not valuable for the integration of practical applications [
5,
6].
On the other hand, conventional optical fiber based SPR sensors have progressed a lot due to their high resolution, a miniaturized structure, and reliability in remote sensing applications [
7,
8]. But these types of fibers are single mode or multimode fibers which consist of core and cladding section; therefore, the external part of the fiber is chemically or physically isolated to reveal the core region to enhance the linking among the surface plasmon polariton (SPP) mode and core mode [
9,
10].
The invention of photonic crystal fiber (PCF) was a prominent breakthrough in optical fiber technology which is utilized in various fields such as spectroscopy, biomedical, and industrial machinery [
11]. PCF has fascinated attraction over the conventional optical fiber due to their remarkable performance such as large effective area, nonlinearity, low loss, multi-parameter sensors, endlessly single mode, tunable dispersion, high index contrast and extensive range of environmental monitoring [
12,
13]. Therefore, the combination of PCF-based SPR sensing mechanism pays much attention to many researchers due to their flexibility structures and controllable light [
14]. The air holes in the core and cladding region moderate the effective refractive index of the core which enhances the phase matching between the SPP and core modes [
15]. The selection of plasmonic elements is also an important factor for PCF-SPR sensor, commonly there are gold, graphene, silver, copper, aluminum, and titanium nitride [
16,
17]. Among these materials, mostly silver and gold have created extreme and developing concern in the research study of PCF-SPR sensor to develop the precious optical features [
18]. Silver approaches a high sensitivity, but it becomes oxidized in the presence of humid environment that can be improved by coating a layer of other elements like graphene, gold, or titanium oxide. Contrary, gold is chemically stable, biocompatible, and more appropriate for microfluidic liquids sensing [
19].
Recently, many groups of researchers have proposed microstructure fibers based SPR sensor with various configuration of PCF where analytes or layer of plasmonic materials placed inside or outside the fiber to measure the different parameters [
20]. In 2017, Weng et al. introduced the D-shape hole with dual cladding photonic crystal fiber based on surface plasmon resonance for sensing the temperature parameters [
21]. In 2018, Ahmed et al. proposed the plasmon resonance PCF biosensor based on Titanium Nitride (TiN) to achieve the RI sensitivity as 7700
nm/RIU [
22]. Furthermore, a dual channel photonic crystal fiber based on gold plasmonic material has been numerically investigated with the sensitivity of 3750
nm/
RIU for RI measurement [
23]. In 2020, a multi-channel PCF sensor based on surface plasmon resonance with multi-analytes channel has been studied with the sensitivity of 3083
nm/RIU [
24]. Additionally, for the detection and quantification of coronavirus a photonic crystal fiber based on the plasmonic material is theoretically examined in 2022. The obtained sensitivity was 2009
nm/RIU, 2745
nm/RIU, 1984
nm/RIU, for the analytes RBD, RNA, and IgG respectively [
25].
From the above literature review it has been seen that mostly researchers work with the complex structure of PCF’s which are not available commercially and need specific arrangements to fabricate which could be lengthy and costly. Further, the single or multiple plasmonic materials are utilized at the same structure of PCF to measure the sensitivity whereas the comparative study still is missing among the two materials at the same structure of PCF.
In this study, we introduced a broad analysis of commercially available solid core photonic crystal fiber based on plasmonic materials for liquid refractive index sensing. The Finite element method (FEM) has been applied to compute the sensing performance and the guiding properties of light spectrum. The most highly sensitive plasmonic materials silver and gold are utilized with the analytes refractive index range of 1.32-1.40. The proposed sensor is explored by finding the different optical properties of confinement loss, wavelength sensitivity, amplitude sensitivity, power spectrum and transmission intensity which creates the shift of peaks wavelength. The maximum sensitivities achieved as 1.15 μm/RIU, 1.10 μm/RIU for gold and silver materials respectively. On the other hand, by analyzing the power spectrum the obtained sensitivity was 513μm/RIU. Moreover, the comparative study of plasmonic particles and the effect of size of air holes and metal layers also has been investigated.
Materials and Methodology
The PCF-SPR design sensor involves cladding air holes with a refractive index of 1, and the background materials is silica whose refractive index can be analyzed by the sellimeier equation. The light dispersion of the sellimeier method can be calculated by the following equation [
28].
where
n is the effective refractive index of the pure silica glass, the values for the constant variables are
A1 = 0.6961663,
A2 = 0.4079426,
A3 = 0.8974794,
B1 = 0.0684043,
B2 = 0.1162414,
B3 = 98.96161.
λ is the operating wavelength in (
μm).
Moreover, silver/gold has been selected as one of the promising candidates for plasmonic material which improves the detection capacity of the sensor. The advantages of these materials are no interband shift at visible wavelength, narrow resonance peak, and low optical damping. The material dispersion of plasmonic materials to realize the metal-dielectric behavior is illustrated by the Drude–Lorentz model [
29].
where
ɛAu is the permittivity of gold,
ɛ∞ = 5.9673 is the permittivity of gold at high frequency,
ω =
2πc/
λ is known as the angular frequency where
c is the velocity of light,
ωD = 2113.6 THz ⨯ 2π is known as plasma frequency,
γD = 15.92 THz is known as damping frequency and
Δɛ = 1.09 is denoted by the weighting factor.
ΩL= 650.07 ⨯
2π THz and
ГL =104.86 ⨯
2π THz are known as Lorentz oscillator strength and spectral width, respectively.
In the simulation study, the finite element method (FEM) based on the COMSOL Multiphysics Software is applied to investigate and evaluate the light transmission of the proposed PCF-SPR sensor. In computation, the physics-controlled way of fine mesh analysis is studied which consists of the number of 499097 solved elements for both domain and boundary elements. The distribution of electric fields in the core mode for
X and
Y polarization has been shown in
Figure 2a,b respectively which illustrates the maximum light confining into the core region. On the other hand,
Figure 2c,d shows the light spectrum for the SPP modes. The arrows indicate the direction of light in
X and
Y polarization respectively.
A broadband wavelength ranges from 1.30
μm to 1.70
μm is carried out for the whole simulation analysis. By using the silver material
Figure 3a demonstrates the distribution of the light transmission among the SPP and fundamental core modes in which the RI of analytes, size of air holes, and the thickness has been fixed at RI=1.33,
d=1.0
μm, and
ts=26
nm respectively. In this graph, the confinement loss (CL) and the real part of the effective refractive index (
neff) are taken from the right and left side respectively. The red line exhibits the loss peak for core mode whereas the black and blue solid lines represent the effective refractive index of fundamental modes and SPP modes respectively. It has been noticed that at the wavelength of 1.44
μm, the same refractive index of SPP and core mode intersect with each other which creates the resonance wavelength. At the resonance wavelength, the maximum light is transferred from the fundamental mode to SPP mode, where the sharp loss peak is monitored which is called phase matching condition. The unidentified samples can be detected by varying this peak to a longer or a shorter wavelength with the refractive index of different analytes. As the wavelength is going to increase the effective refractive index for the core and SPP mode is steadily decreased. So, the maximum refractive index is achieved at the minimum wavelength. Likewise, the short dash line in
Figure 3b illustrates the loss spectra with respect to wavelength for gold material. Here the same parameters have been fixed like silver analysis but the resonance peak shifts towards the higher wavelength with higher refractive index and minimum loss. Usually, the core and SPP modes work at visible and infrared radiation respectively. The propagation of light has been evaluated for core-
X and core-
Y polarization spectrum. Both polarization spectrums have the same mechanism of light transmission. Therefore, core-
X-polarized mode is considered to explore the sensor schemes for further study.
Results and Discussions
Moreover, the change of direction of light coupling shifts between the core and SPP mode happens at the resonance wavelength. Basically, the computation and examinations of the proposed PCF-SPR sensor have been explored by CL, one of the most important aspects in evaluating the function of the sensor. The CL of the fundamental mode can be calculated by the following equation [
30].
Where is the free space of wave number, λ is the operating wavelength (μm), and is the imaginary part of effective refractive index.
To realize the sensing mechanism of the proposed PCF-SPR structure, the properties of the analytes RI have been examined. The effective refractive index of the SPP mode is highly responsive to the surrounding environment. The small variation of RI of analytes makes a significant shift of the real part of core mode. The comparatively confinement loss of silver and gold as a function of wavelength due to the alteration of analyte RI from 1.32 to 1.40 is shown in
Figure 4. This graph depicts the solid lines to indicate the silver resonance wavelength. At the minimum RI (1.32) with sliver material, the small loss peak (162
dB/cm) is observed. By increasing the RI of analytes, the loss peak is going to increase at the higher RI (1.40), the maximum loss peak (340
dB/cm) is obtained. The resonance wavelength shifts move towards the shorter wavelength with the increment of analytes RI. Similarly, the same phenomena of loss spectra with short dash lines of gold material have been shown in
Figure 4. The resonance wavelength shift attained at the higher wavelength as compared to silver analysis also the higher loss peak (386
dB/cm) is achieved at the RI of 1.40. During the analysis of loss spectrum with the variation of analytes RI, the rest of parameters like metal layer thickness (t
g=26
nm) and air holes (d=1.0
μm) were kept constant.
Furthermore, the effect of air holes and the thickness of plasmonic layers has been evaluated.
Figure 5a describes the solid and short dash lines of the variation of silver and gold layer (22
nm to 34
nm) with the change of each step 4nm respectively. It seems that by increasing the thickness of the silver and gold layer the resonance wavelength is moving towards the shorter wavelength and comparatively both have higher resonance wavelength. In both cases (silver and gold) at the thickness of (t
s,t
g=26
nm) the highest loss peak (269
dB/cm) appears, and at the maximum thickness (t
s,t
g=34
nm) the minimum loss (119
dB/cm) and (135
dB/cm) is achieved respectively. As a result, the alteration of plasmonic layers has a great impact on the sensing phenomena.
Correspondingly, in
Figure 5b different sizes of air holes (0.96
μm to 1.0
μm) with the interval step of 0.02
μm have been considered for both silver and gold material. Here the maximum loss peak (252
dB/cm) is attained at the size of (1.0
μm) for gold metal.
Sensitivity is the most important factor to validate the actual performance of the sensor that depends on the peak wavelength shift and the variation of the analytes RI. The following equation can be employed to determine the sensitivity (
S) of analytes [
31].
Where, is the resonance wavelength shift, and is the variation of analytes RI.
Another meaningful feature is the Resolution (
R) of the sensor which is applied to identify the minor variation of analyte RI. The resolution of a sensor can be determined by the following equation [
32].
Where, represents the difference between two resonance peak wavelength, is the alteration of analytes RI, and is the maximum peak of resonance wavelength.
Furthermore, for practical execution of sensors, the inquiry of amplitude approach is modest and efficient. Amplitude sensitivity can be estimated by operating the following equation [
33].
Where
is the confinement loss as a function of the wavelength and RI of analytes.
is the difference between two loss peaks under two different refractive indexes at the same wavelength. The values of amplitude sensitivity have been calculated by tracing
Figure 4.
Figure 6a implies the amplitude sensitivity of the refractive index of different analytes for silver layer. At the refractive index of (1.36) the maximum amplitude sensitivity (-46
RIU-1) is achieved. On the other hand,
Figure 6b indicates the amplitude sensitivity for the gold layer. Here the maximum amplitude sensitivity (-32
RIU-1) is obtained at the resonance wavelength (1.46
μm) with the RI of (1.38).
Figure of merit (FOM) is the significant attribute to value the function of any type of sensors. The following equation is employed to find the FOM [
34].
Where
is the wavelength sensitivity measured in
μm/RIU, and FWHM is the full width at half-maximum calculated in nanometers from the resonance wavelength.
Figure 7a shows the function of FOM with the variation of analytes RI for silver layer. The maximum FOM attained as 25
RIU-1 at the analyte RI of 1.34. Next in
Figure 7b the variation of analyte RI for gold layer has been depicted. Here the maximum FOM 31
RIU-1 obtained at the RI of 1.32 which is comparatively higher than the silver analysis.
On the basis of the above equations (4,5,6,7) the wavelength sensitivity, amplitude sensitivity, resolution, and FOM of proposed PCF-SPR sensor has been calculated.
Table 1 explains the comparative results among the silver and gold in which the performance of the gold metal secured higher sensitivity than silver layer.
Additionally, the birefringence, coupling length, power spectrum, and transmission intensity also has a great influence to validate the working function of the sensors. Therefore, these parameters have been calculated based on the difference of effective refractive index of the X and Y polarization mode. Here the explanation of the most important aspects such as power and transmission intensity have been considered.
The optical power variations operate from one mode to another, along with the proposed sensor. The expression of birefringence and coupling length define the power spectrum. The power spectrum ensures the sensor ability, which is measured by using the following equation [
35].
Where B is the birefringence, L is the coupling length, and λ is the operating wavelength.
Figure 8a shows the power spectrum as a function of wavelength with the variation of analytes RI (1.32 to 1.40). By increasing the RI of analytes, the power peak spectrum moves towards the shorter wavelength. The maximum sensitivity obtained is 1.13
μm/RIU at the refractive index of 1.32.
The transmission intensity plays a significant role in calculating the sensitivity of the proposed PCF-SPR sensor. The change of intensity propagation has been manipulated. The whole plasmonic impression extensively fluctuates on the transmittance. However, the wavelength sensitivity of the designed sensor can be figured out through the alteration of the sharp peak of the transmittance curve with the change of the refractive index. The transmittance spectrum can be evaluated by using the following equation [
36].
where
Tr is the transmission of light,
Pin is the maximum power as a fraction of
Pout.
The transmission spectrum with respect to wavelength by varying the analyte RI has been depicted in
Figure 8b. This graph explains the variation of resonance wavelength to shorter wavelength as the refractive index increases. With the increases of analyte refractive index, the transmission is going to higher wavelength. For the transmission spectrum the attained sensitivity is 1.21
μm/RIU.
A linear reaction is a valuable mechanism for sensor execution as it allocates to resolve the unidentified quantities by regular data extrapolation. Therefore, by considering the good results of
Figure 4 and the
Figure 8a (power spectrum) with gold material has been explored for linear regressions.
Figure 9a illustrates the linear fitting curve of loss function with respect to the resonance wavelength and the variation of analytes RI. In this graph the slope values consider the sensitivity of the sensor which is obtained as 513
μm/RIU with the adjacent R
2 of 0.997. Similarly,
Figure 9b exhibits the linear regression of power spectrum as a function of wavelength shift and the variation of analyte RI. The achieved sensitivity was 0.706
μm/RIU with an R-square value of 0.979.
A comparison study of the proposed structure and the previous work done based on RI sensing has been drawn in
Table 2, which consists of the sensitivity results. In the overall analysis, the PCF-SPR sensor has a remarkable impact in refractive index sensing, particularly gold has higher sensitivity than silver layer. Moreover, this designed structure enhances solidity and stability for sensing applications. By tuning the structural parameters and sensing materials, the results can be further improved.