1. Introduction
Surface plasmon resonance (SPR) is a kind of oscillation which gathers electrons from a conduction band and interacts with the periodical varying electric field generated from incident light [
1]. In this optical effect, incident light is introduced with an appropriate incidence angle, which brings photons to the surface of plasmonic metal materials. These photons interact with surface-located electrons and further cause resonance where these original surface-located electrons appear to propagate parallel to the surface of plasmonic metal materials and create energetic plasmon electrons through non-radiation excitation [
2,
3]. When SPR occurs, a minuscule change in the refractive index (RI) in analyte or another medium in the environment in simple or complicated systems such as single-object or multi-physics can alter the phase-matching condition. Change of condition further influences the existence of the SPR effect [
4,
5]. This renders the measurement of analytes and related sensing technologies feasible for various systems, particularly in applications such as biomedical detection services in the face of global public health crises like the Coronavirus [
6].
Photonic Crystals can be utilized as integrated optical devices which enhance interactions between light and matters [
7]. Many applications have been demonstrated, such as all-optical logic gates [
8], filters [
9], highly efficient modulators [
10], and sensors. For photonic-crystal-based sensors, there are two types. One type is based on an array of different number ring-slot and input–output line defect coupling waveguides, which places light sources vertically on the rods to connect waveguides from both input and output paths [
11]. This type has a small sensitivity value but has a narrow curve width on the nanometer scale. The second type is utilized in optical fibers. In recent decades, with a large quantity of research concentrating on photonic crystal fibers (PCFs) [
12], it has been well concluded that PCFs can be effectively applied with the SPR effect for a reliable reason, as PCFs perfect optical properties along with different geometries and materials to choose for design [
13,
14]. A large amount of work has been carried out to design SPR-based PCF sensors to achieve good performances in the fields of measurement and detection. Some typical applications for single-object or multi-physics systems are known as liquid RI, salinity [
15], temperature [
16], magnetic field [
17], etc.
There are a few types of SPR-based PCF sensors. The realization of the first type depends solely on various topological geometries and different plasmonic metal materials chosen for single-parameter measurement. In 2018, Chen et al. proposed a D-shaped SPR-based PCF applied as an RI sensor with numerical studies. Unlike normal D-shaped structures, this study employed an open-ring channel coated with a gold film to excite the SPR effect. Simulation results indicated that the sensor achieved a sensing range of RI from 1.20 to 1.29. When the analyte’s RI exceeded 1.25, the anti-crossing effect started to appear and the peak value of the loss curve remained nearly constant with a change in RI. The proposed design achieved a maximum spectral sensitivity of 11,055 nm/RIU and a high resolution of about
RIU when RI was selected as 1.29 [
18]. Some unusual shapes are also proposed. In 2023, Chao et al. proposed a novel trapezoidal-shaped temperature sensor. The fiber structure conducts two layers of elliptical air holes with different parameters and a polished trapezoidal-shape surface where thin silver and silica dioxide are deposited. Simulation results indicated that the sensor obtained a maximum temperature sensitivity and resolution up to 5200 pm/
C and
, respectively. Moreover, the temperature sensing range was from 10 to 60
C [
19].
The second type is focused on measurement in multi-channels to reach the scenario of multi-parameter sensing without considering cross influence. In 2020, Bing et al. proposed a dual-channel SPR-based PCF sensor that possessed the advantage of improving the detection efficiency of samples on different channels. In this design, a horizontal polished structure was proposed. The fiber structure was designed to be immersed directly in the liquid analyte to reduce coating difficulty and prevent microfluidic filling. Through the optimization of the dielectric layer, the ability to identify channels was efficiently improved. The wavelength sensitivities of two channels achieved 11,600 nm/RIU and 10,600 nm/RIU, respectively. The proposed structure was capable of simplifying the process complexity in multi-sample detection [
20].
The third type focuses on a simultaneous performance in multi-physics systems for a multi-parameter sensing scenario. In 2019, Han et al. proposed a PCF structure with dual surface plasmon resonance to realize the measurement of temperature and strain while addressing the non-cross-sensitivity problem simultaneously. In this design, two holes in a symmetric relation of the fiber geometric were coated with goad coatings. One of those two holes was infiltrated with temperature-sensitive liquid in order to measure temperature. Simulation and calculation results indicate that the proposed design achieved simultaneous measurement of temperature and strain, with temperature sensitivies of −6.83 nm/
C and a strain sensitivity of 1.30 pm/
with a scenario when diethylene glycol was chosen as the temperature-sensitive liquid [
21].
In this article, an optimized novel SPR-based PCF sensor structure with a polished-D-shape is proposed for measurement of RI with high values of maximum and average sensitivity. Sensing performance of the fiber sensor structure has been analyzed through numerical simulations with approaches based on the finite element method (FEM). In the correspondence, different capillaries are gathered and processed to a polished-D-shape silica structure with nano-scale gold material coated on flattened surface. By applying a thin gold film and a solid silica-background core area to excite and enhance the SPR effect, the loss curve shifts towards longer wavelengths when the RI value increases. By optimizing structural parameters, the values of maximum and average sensitivity reach 32,000 and 12,167 nm/RIU, respectively, along with a sensing coverage of RI from 1.26 to 1.32, in steps of 0.01 RIU. Moreover, the resonant wavelength of the proposed design ranges from 1810 to 2540 nm.
Label-free sensors are systems that are capable of detecting physical, chemical and biological matter in micrometer and nanometer scales [
22]. These devices have extraordinary sensitivity along with a fast time applied for measurement in numerous scenarios [
23]. An SPR-based sensor is a kind of label-free sensor. Benefiting from these outstanding performances along with an intuitive geometrical structure, the proposed fiber sensor will have an advantage in the field of organic and biological detection and related applications.
This article is organized as follows:
Section 1 is a brief introduction to the research background and the relationships between these abbreviations.
Section 2 provides the geometry of proposed model and some parameters utilized in the study.
Section 3 shows a comparison and corresponding analysis of related models concentrating on different sensing performances.
Section 4 offers a detailed description of the effect of geometrical parameters on sensing performances.
Section 5 conducts a simulation on the final optimized model and demonstrates a baseline comparison with similar studies. The summarized conclusions are presented in
Section 6.
2. Model and Method of the Design
Figure 1a,b illustrate a 3D schematic and a 2D sketch of the proposed fiber sensor, respectively. Panel (i) demonstrates mode field distributions of the core region area A in
-polarization, which is utilized to measure RI.
The proposed cylinder structure is primarily composed of an area of an air hole array with a triangle lattice arrangement. The lattice spacing, which is interpreted as the distance between two nearby geometrical centers of the air holes in the array, is represented by . In the initial design, the structure was simplified, with all air holes having same radius as . To excite the SPR effect, part of yjr silica structure is polished to form a flattened plain and nano-scale gold material is coated on the surface. The thickness and width of the coating are presented by and . For convenience, these two parameters are abbreviated as t and L.
In further improvement, six air holes at the tips of the array are removed and four air holes at the edge of the upper and lower array are narrowed to a smaller radius denoted by
, which aims to decrease the silica-background area between the air hole array and external part in order to concentrate more energy on the solid-silica-background core area and further improve sensing performances. The solid-silica-background core area is named core area A, which generates the core mode. Area A is formed by removing two air holes, which is illustrated in the sky-blue rectangular area in
Figure 1b.
In related works, it has been found that elliptical holes demonstrates a superior performance to circular holes as they are capable of remarkably affecting the coupling energy between the core mode and the surface plasmon polariton mode (SPP) [
24]. Thus, the air hole on the left side of area A is replaced with a horizontally placed elliptical hole with semi-major and semi-minor axes, denoted by
and
. For convenience, these two parameters are abbreviated as
a and
b. A tiny air hole with a radius denoted by
is added on the right side of area A in order to further strengthen the SPR effect. The adjacent two air holes of the elliptical hole are magnified to a larger radius denoted as
in order to extrude the core region and to make sure most of the energy from area A interacts with the thin gold film on the right side, thereby improving sensing performances. On the external side of the silica structure, liquid analyte is introduced.
Figure 2a illustrates the 2D schematic of the proposed fiber sensor in the process of fabrication. For the proposed structure, it can be practically realized with the well-known stack-and-draw fabrication technique [
25]. In manufacturing procedures, an external container is applied to gather a variety of different rods and capillaries at the start-up stage. In the correspondence, thick-wall, thin-wall, elliptical-core and square circular-core fiber capillaries are utilized as basic units to form initial holes, large holes, elliptical core and tiny air hole, respectively. Some simple auxiliary solid rods and two little-bit-complicated auxiliary polished solid rods are applied as the removed air holes in order to fix the position inside the container and create core region area A during the assembly process. An elliptical-core capillary is used to create the elliptical hole on the left side of area A. In addition, a special circular-core square capillary is embedded in these two auxiliary polished solid rods in order to form the tiny air hole and fix its position in the fabrication process. After gathering these rods and capillaries together in a symmetry-shaped hexagonal external container, liquid silica is infiltrated into the space among these micro structures. At this moment, the start-up hexagonal structure is made. The structure is then embedded into a hollow cylinder external container and filled with liquid silica once again to form a cylinder-shaped fiber structure. Then, the structure is shaped into a polished D-shape fiber, and plasmonic metal material, namely nano-scale gold, is coated onto the flattened surface. Subsequently, the final structure of the proposed design is fully realized.
Figure 2b demonstrates a cross-sectional FEM mesh of the proposed fiber sensor in simulations. In the correspondence, the finite element method is used for solving eigenmode with a quantity of square and triangle grids in numerical simulations. Furthermore, a perfectly matched layer (PML) is designed at the exterior part of the analyte area to cut off regions beyond the analyte area and limit the computational area for numerical simulation by absorbing all radiation energy at the edge of the PML region. The thickness of PML is designed as
2.0
m. The numerical studies and results are simulated and calculated with the help of Wave Optics Module in the software COMSOL Multiphysics 6.0.
Fused silica is chosen for its superior properties and physical performances. It is designed as the background base material of the sensor structure. The material refractive index of fused silica is obtained by a third-order Sellmeier Equation (1),
where
is the wavelength in free-space, measured in nanometers. The value of coefficients
are given in
Table 1 in a non-dimensional form. These coefficients are suitable for the wavelength range in the this work [
26,
27].
Nano-scale gold is a kind of plasmonic metal material. It is coated at the polished surface near the core region to excite the SPR effect. The dispersion of nano-scale gold is obtained through the Lorentz–Drude model expressed by Equation (2) [
28,
29].
where
stands for the dielectric constant in a scenario where frequency is ultra-high, and
is the angular frequency of the light propagating through the PCF system.
is the weight factor.
and
represent the plasma and damping frequency, respectively.
and
are the frequency and spectral width of the Lorentz oscillator, respectively.
and
are
rad/s and
rad/s in angular frequency.
It is widely known that PCFs have multiple types, including high-birefringence (Hi-Bi) PCFs and SPR-based PCFs. Most Hi-Bi PCFs have cladding layers without plasmonic metal materials. This type has high index contrast and can be applied as polarization-maintaining fibers, which are capable of eliminating the effect of polarization mode dispersion or of stabilizing the operation of optical devices [
30]. SPR-based PCFs are coated with plasmonic metal materials. It works on the basis of interactions between electrons on the surface of the plasmonic metal materials along with the evanescent field [
31,
32]. In SPR-PCF sensor systems, when light with an appropriate incident angle propagates through the structure, some part of the light penetrates through the air hole array and undergoes interactions with plasmonic metal materials coated on the polished surface. The interactions generate free electrons excited from surface. SPR-PCF sensor systems use different-order SPP modes and core modes to produce the SPR effect [
33]. When the evanescent field of the incident light reaches a frequency that matches with the frequency of free electrons generated from the surface of plasmonic metal materials, the resonance effect is thus fully generated along with the surface plasmon wave [
34]. This phenomenon is characterized by confinement loss (CL), which is a good metric to represent the loss when the light propagates through the system. The
value is obtained by calculation calculation Equation (3) [
35],
where
represents the effective RI of the coupled mode. The value of
is measured in dB/cm. When the phase-matching condition is achieved, a loss peak occurs in the loss curve. The resonance peak will witness a red or blue shift with a change in the RI value of the dielectric medium, such as a sample or analyte. which is adjacent to the surface. The change of peak value reflects the concentration of the medium [
36].
Figure 3a demonstrates the
curve of the coupled mode in both
x- and
y-polarization along with the real part of the effective RI of the coupled mode in both the core mode and the SPP mode with RI selected as 1.28 with a list of geometrical parameters in the proposed structure in the initial stage:
a = 0.9
m,
b = 0.4
m,
= 0.2
m,
= 0.6
m,
= 1.1
m,
= 1.3
m,
t = 60 nm,
l = 5.6
m,
= 2.8
m. Furthermore,
Figure 3b–d illustrate the mode field distributions of the coupled mode in
x-polarization when the wavelength is selected as 1760, 1920, and 2300 nm, respectively. From
Figure 3a, the core mode in area A is completely coupled with the SPP mode on the surface of the polished nano-scale Au-coated plane when the wavelength is selected as 1920 nm. In some parts of the wavelength regions, the value of confinement loss is close to zero. Thus, it can be approximated that resonance almost does not exist. A reasonable explanation for this phenomenon is considered to be that in those wavelength regions, most of the energy is bound in the core region, as depicted in
Figure 3b and d. When it reaches an appropriate resonant wavelength denoted by
, the SPR effect is fully excited, with the core mode completely coupled with the SPP mode. During this time, the mode field distribution and energy distribution transition from the core area to the polished surface of the nano-scale Au coating, accumulating the largest amount of energy in a small space with a high energy density, as illustrated in
Figure 3c. This density reaches 160, which is higher than the values of 90 and 100 in
Figure 3b and d, respectively. At this point, the phase-matching condition is fully achieved, observed at the intersection between the yellow solid curve and the green solid curve in
Figure 3a. At this intersection, the
curve peaks. Furthermore,
Figure 3a shows that modes coupled in different polarization states exhibit a significant difference. In the
x-polarization, the curve has a large value in terms of dB/cm-scale. In the
y-polarization, the CL value is close to zero. This further indicates that the proposed design achieves well-performing single-polarization sensing. In the subsequent parts of this work, we will select
x-polarization for research and analysis.
The variation in RI has an influence on the effective RI of the coupled mode in the SPR-excited region in the
x-polarization. This, in turn, alters the mode field distributions in the area. The corresponding value of the resonance wavelength also changes. Eventually, the
curve will exhibit either a red-shift or a blue-shift, depending on the variation in the peak value of these curves. Therefore, it is significant to seek evaluation schemes to measure the combined impact of the varying structural parameters on the affected sensing performances. Under this scenario, it is appropriate to characterize the relationship between changes in the
curves and variations in RI. Therefore, the RI sensitivity, denoted as
, can be a good metric, and it is obtained through the following equation [
37],
where
is the change of resonant wavelength value and
denotes the varying RI with a step of 0.01 RIU. The unit of
is nm/RIU.
Also, amplitude sensitivity, which is denoted as
or
, is introduced as a another useful metric. It can reflect the relationship between the curve intensity of the
spectra along with the variation sequence of RI. It is defined as Equation (5), the equation below,
where
is the change of confinement loss value and
denotes the initial confinement loss. Metric
has a unit of
.
In addition, another important metric can combine sensitivity performances along with the shape of the
curve in the proposed sensor. This is the figure of merit, which is denoted as
. This metric focuses on the full width at half of the maximum (FWHM) value. The calculation of this metric is presented by Equation (6),
where the unit of
is nanometers.
Figure 4 presents a possible implementation of the experimental setup of the SPR sensor for detecting analyte RI. Initially, a broadband light source is employed to direct initial light through a polarizer. In this work, the light exits the polarizer as
x-polarized light. Subsequently, the light passes through the light controller device and enters a single-mode fiber (SMF). The SMF is capable of receiving light and transitioning it into an environment suitable for PCF. Molecules in the liquid analyte are utilized for measurement and sensing. These molecules pass through a mass flow controller, serving as the inlet to enter the sensor, and then as the outlet to flow out of the sensor. These interactions lead to the SPR effect and cause a red or blue shift in the loss curve, which can be observed using an optical spectrum analyzer (OSA). The relevant data can be visualized via a computer [
38].
3. Geometrical Comparisons
In this section, a variety of related geometrical designs are compared and analyzed.
Figure 5 shows the comparison of four 2D schematics of the design. In
Figure 5a, a symmetrical construction design is applied to form a dual-core sensor structure with an elliptical core set at the central location to isolate these two core areas.
Figure 5b is the initial design. In
Figure 5c, a small air hole is added on the left side of the air hole array to concentrate the mode field distribution on core region area A. In
Figure 5d, two circular air holes are replaced with elliptical air holes pointing to the core area in order to alter the region for the SPR effect and change the sensing performance.
Figure 6 demonstrates
curves of the coupled mode in
x-polarization of four compared structures with different ranges of analyte RI.
In
Figure 6a, with an increase in the sequence of RI with a changing range from 1.30 to 1.36, the
curve of the coupled mode in the core area in
x-polarization undergoes a long-distance red-shift. Additionally, the corresponding peak value of the
curve shows a trend from increase to decline. With the red shift,
follows a changing sequence: 1970, 2025, 2100, 2175, 2260, 2370, 2580 nm, with RI respectively selected as 1.30, 1.31, 1.32, 1.33, 1.34, 1.35, 1.36. According to calculation using Equation (4), the corresponding
values are obtained as 5500, 7500, 7500, 8500, 11,000, and 21,000 nm/RIU, respectively. The sensing range of RI in this scenario is from 1.30 to 1.36. Additionally, the average value of
is calculated as 10,167 nm/RIU.
In
Figure 6b, with an increase in the of RI with a changing range from 1.26 to 1.33, the
curve of the coupled mode of the core area in
x-polarization undergoes a significant red-shift. In addition, the corresponding peak value of the
curve demonstrates a trend from increase to decrease. With the red shift,
has a changing sequence: 1820, 1880, 1940, 2010, 2090, 2220, 2360, 2600 nm with RI set as 1.26, 1.27, 1.28, 1.29, 1.30, 1.31, 1.32, 1.33, respectively. According to Equation (4), the corresponding
values are obtained as 6000, 6000, 7000, 8000, 11,000, 14,000 and 24,000 nm/RIU, respectively. The sensing range of RI in this model is from 1.26 to 1.33. In addition, the average value of
is calculated as 11,000 nm/RIU.
In
Figure 6c, with an increase in the sequence of RI within a changing range from 1.26 to 1.32, the
curve of the coupled mode in the core area in
x-polarization witnesses a significant red-shift. Additionally, the corresponding peak value of the
curve shows a trend from rise to decrease. With the red shift,
follows a changing sequence: 1850, 1900, 1970, 2050, 2150, 2300, 2540 nm, when RI is 1.26, 1.27, 1.28, 1.29, 1.30, 1.31, 1.32, respectively. According to calculations using Equation (4), the corresponding
values are obtained as 5000, 7000, 8000, 10,000, 15,000, and 24,000 nm/RIU, respectively. The sensing range of RI in this scenario ranges from 1.26 to 1.32. Additionally, the average value of
is calculated as 11,600 nm/RIU.
In
Figure 6d, with an increase in the sequence of RI within a changing range from 1.28 to 1.34, the
curve of the coupled mode in the core area in
x-polarization witnesses a significant red-shift. Additionally, the corresponding peak value of the
curve shows a trend from increase to decline. With the red shift,
follows a changing sequence: 1900, 1960, 2020, 2100, 2220, 2400, 2620 nm, with RI respectively chosen as 1.28, 1.29, 1.30, 1.31, 1.32, 1.33, 1.34. According to calculations using Equation (4), the corresponding
values are obtained as 6000, 6000, 8000, 12,000, 18,000, and 22,000 nm/RIU, respectively. The sensing range of RI in this scenario ranges from 1.28 to 1.34. Additionally, the average value of
is calculated as 12,000 nm/RIU.
Figure 7 shows
curves of coupled mode in
-polarization of four compared structures with different range of analyte RI.
In
Figure 7a, with an increase in the sequence of RI within a variation range from 1.31 to 1.35, the
curve of the coupled mode in the core area in
x-polarization witnesses a long-distance red-shift. The
value, which represents the corresponding wavelength of the peak of the
curve in the negative semi-axis, follows a changing sequence: 2120, 2180, 2030, 2260, 2380, 2600 nm, when RI is 1.31, 1.32, 1.33, 1.32, 1.35, respectively. According to calculations using Equation (5), the corresponding absolute values of
are obtained as 99.69, 126.8, 176.9, 89.79, and 80.88 RIU
, respectively.
In
Figure 7b, with an increase in the sequence of RI within a variation range from 1.26 to 1.32, the
curve of the coupled mode in the core area in
x-polarization witnesses a dramatically red-shift. The
value follows a changing sequence: 1920, 1960, 2020, 2100, 2220, 2380, 2600 nm, with RI respectively set as 1.26, 1.27, 1.28, 1.29, 1.30, 1.31, 1.32. According to calculations using Equation (5), the corresponding absolute values of
are obtained as 75.92, 85.46, 112.1, 134.1, 102.6, 67.27, and 3.213 RIU
, respectively.
In
Figure 7c, with an increase in the sequence of RI within a variation range from 1.26 to 1.31, the
curve of the coupled mode in the core area in
x-polarization witnesses a drastically red-shift. The
value shows a changing sequence: 1920, 1980, 2040, 2160, 2300, 2400 nm, with RI respectively selected as 1.26, 1.27, 1.28, 1.29, 1.30, 1.31. According to calculations using Equation (5), the corresponding absolute values of
are obtained as 79.75, 107.8, 145.6, 99.71, 70.7, and 9.624
, respectively.
In
Figure 7d, with an increase in the sequence of RI within a variation range from 1.28 to 1.33, the
curve of the coupled mode in the core area in
x-polarization witnesses a significant red-shift. The
value has a changing sequence: 1980, 2040, 2100, 2220, 2400, 2500 nm, with RI respectively chosen as 1.28, 1.29, 1.30, 1.31, 1.32, 1.33. According to calculations using Equation (5), the corresponding absolute values of
are obtained as 86.39, 107, 147, 130.6, 73.49, and 15.25
, respectively.
For the model illustrated in
Figure 5a, this design is a dual-core type of SPR-PCF sensor, which is significantly different from the other three single-core sensing models. As a result, it exhibits different loss curve distributions, including curve width and curve intensity, as well as the sensing range for refractive index (RI). In the models shown in
Figure 5b and c, the most notable disparity lies in the existence of a small air hole on the left side. Thus, the loss curve distribution is almost the same, with the exception of an additional RI sensing range of 1.33, along with numerically different sensing performances.
For the models illustrated in
Figure 5c,d, the most significant difference is the variation in the core area geometry, specifically the number of elliptical air holes used. Consequently, the loss curve distributions for these two models exhibit visible differences, particularly in terms of the sensing range of RI and curve width. Taking into consideration both the fabrication feasibility and the pursuit of improved performance, including RI sensing range and sensitivities, the schematic in
Figure 5b is chosen as the final proposed geometrical structure for analysis and optimization.
4. Simulation Results and Formal Analysis
Figure 8 illustrates the
curve of the coupled mode in
x-polarization when
is chosen as 2.6, 2.8, 3.0
m, and RI is changing from 1.26 to 1.33 with a step of 0.01 RIU.
From
Figure 8, as the value of
rises from 2.6 to 3.0
m, the
spectra begin to red-shift along with a decrease in curve intensity. The sensing coverage of RI widens to 1.33 when
is set at 2.8
m, and it narrows to 1.31 when
is chosen as 3.0
m. Moreover, with the increase in the value of
, the envelope of
curves corresponding to varying RI shifts to the red and the curve intensity decreases.
The reason for these changes is considered to be that as the value of increases, the core area A expands, and the arrangement of the air hole array becomes looser. This alteration affects the mode field distributions of the coupled mode in x-polarization within area A. Consequently, it modifies the phase-matching condition and the range of RI sensing, resulting in the emergence of new detected RIs such as 1.32 and 1.33 when increases. Additionally, with the increase in the value of , the expanding core area leads to a portion of the energy not being effectively concentrated in the area. As a result, the intensity of the envelope visibly decreases.
According to calculation using Equation (4), the maximum values of are achieved as 22,000, 24,000, and 22,000 nm/RIU when is respectively set to 2.6, 2.8, and 3.0 m. The calculated results clearly indicate that the maximum value of in the proposed design exhibits a trend of rise followed by a decline as increases from 2.6 to 3.0 m, with a peak value attained at of 2.8 m. Furthermore, the average values achieved are 10,713, 11,000, and 13,400 nm/RIU with set to 2.6, 2.8, and 3.0 m, respectively. Taking into consideration both the maximum values and the coverage of measurable RIs, it is evident that the proposed sensor structure achieves the best sensing performance when is set to 2.8 m.
Figure 9 demonstrates the
curve of the coupled mode in
x-polarization when
is selected as 1.0, 1.1, 1.2
m, and RI is selected as 1.26, 1.31, 1.32 and 1.33.
From
Figure 9, as the value of
rises from 1.0 to 1.2
m, the
curve of the coupled mode in
x-polarization starts to blue-shift and the curve intensity has a noticeable fall when RI is set as 1.26. The
curve of the
x-polarized mode with RI set to 1.31 and 1.32 has a noticeable red-shift and the intensity of spectra shows a trend from rise to decline. Sensing coverage of RI is extended to 1.33 when
is set to 1.0
m. In addition, when
increases, the envelope of
curves associated with varying RI appears to shift to the red along with the occurrence of a trend from increase toecrease in envelope intensity. The reasons for these changes are considered to be that when the value of
becomes larger, core area A is gradually compressed, which concentrates the energy. However, the channel between the core area and the coating surface narrows due to the increasing
value, making the arrangement of the air hole array compact, which causes the intensity of the envelope to trend from ascending to descending. Also, the change of mode field distributions influences the phase-matching condition of the
x-polarized coupled mode on a very small scale, which further causes the sensing scope to be extended to 1.33 when
is chosen as 1.0
m. According to calculation in Equation (4), maximum values of
are achieved as 24,000, 31,000, 29,000 nm/RIU when
is set as 1.0, 1.1, 1.2
m, respectively. It is obvious from the calculation results that the maximum value of
in the proposed design has a trend from rise to decrease when
has an increase sequence from 1.0 to 1.2
m and achieves a peak value with
selected as 1.1
m. Moreover, the average
values achieved are 11,000, 12,000, 12,833 nm/RIU when
is selected as 1.0, 1.1, 1.2
m, respectively. When the maximum
value is given priority to be taken into consideration, the value of
is finally set as 1.1
m when the proposed sensor structure gains an optimal sensing performance.
Figure 10 shows the
curve of the coupled mode in
x-polarization when
is set as 1.1, 1.2, 1.3
m, and RI is selected as 1.26, 1.31, 1.32.
From
Figure 10, with a rising value of
from 1.0 to 1.2
m, the
curves of the coupled mode in
x-polarization with all chosen RI values almost stay still. The curve intensity when RI is chosen as 1.26 has a subtle increase and a shift of the curve does not occur. The curve intensity with RI selected as 1.32 first increases and then decreases. In addition, with the increase of
, the envelope of
curves associated with varying RI becomes narrow and curve intensity shows a decrease with minor changes. The reasons for these changes are considered to be that when the value of
becomes larger, the core area A in
x-polarization is compressed to a certain degree, ensuring more energy interaction with the right-side plasmonic metal material. However, the change of the
value has a small effect due to the limited range of variation in the fabrication process. Also, there is a long distance from
-radius air holes to the core area, which causes the
curve and peak to show a small change. Moreover, the change of the mode field distribution of
x-polarized area A influences the phase-matching condition and causes the
curve and peak to shift slightly. According to calculation in Equation (4), the maximum values of
are achieved as 31,200, 31,000, 31,500 nm/RIU when
is selected as 1.1, 1.2, 1.3
m. It is obvious from the calculation results that the maximum value of
in the proposed design shows a trend from climbing to falling when
has a varying sequence from 1.1 to 1.3
m. Moreover, the calculated average value of
in the proposed design rises when
has a changing sequence from 1.1 to 1.3
m and achieves an ideal large value when
is selected as 1.3
m. Additionally, the average
values achieved are 11,916, 12,000, 12,000 nm/RIU when
is selected as 1.1, 1.2, 1.3
m, respectively. To make sure most of the energy can interact with the right side of the thin gold film with high sensitivity values, the
is finally set to 1.3
m when the proposed sensor structure achieves an optimal sensing performance.
Figure 11 illustrates the
curve of the coupled mode in
x-polarization when
is selected as 0.6, 0.7, 0.8
m, and RI is chosen as 1.26, 1.31, 1.32. From
Figure 11, with a rising value of
from 0.6 to 0.8
m, the
curves with all chosen RI values almost stay still. The intensity of the
curve when RI is selected as 1.26 has a subtle decrease. The curve intensity when RI is selected as 1.32 first increases and then decreases. In addition, with
increasing, the envelope of
curves associated with varying RI become wider and the curve intensity decreases with minor changes. Reasonable explanations for these changes are considered to be that when the value of
becomes larger, core area A does not change, but the mode field distributions are altered by the solid-silica area between the air hole array and the exterior analyte. However, variation of
values has little effect due to the limited range of variation in fabrication along with the far distance between
-radius air holes and the core area, which causes a very small shift of the
peak and curve. Moreover, the changes have little effect on the phase-matching condition of the
x-polarized coupled mode of area A and further causes the
curve and peak to show a slight shift. According to calculation in Equation (4), the maximum values of
are achieved as 32,200, 31,500, 30,000 nm/RIU when
is selected as 0.6, 0.7, 0.8
m, respectively. It is obvious from the calculation results that the maximum
value of the proposed structure gradually drops when
has a variation sequence from 0.6 to 0.8
m. Moreover, the average
values achieved are 12,000, 12,000, 11,067 nm/RIU when
is selected as 0.6, 0.7, 0.8
m, respectively. In this work, efforts are made to ensure most of the energy is concentrated in core area in order to further improve performances. Taking these issues, especially sensing sensitivities, into consideration, the value of
is finally set as 0.6
m when the proposed sensor structure reaches an optimal sensing performance.
Figure 12 demonstrates the
curve of the coupled mode in
x-polarization when
L is selected as 5.2, 5.6, 6.0
m,
t is chosen as 50, 60, 70 nm, and there is a variation of RI from 1.26 to 1.34. From
Figure 12, with an increasing value of
L from 5.2 to 6.0
m, sensing coverage of RI is gradually narrowed. As seen in
Figure 12a,d,g,b,e,h, sensing scope is from 1.26 to 1.33 when
L is set as 5.2
m. When
L is set as 5.6 and 6.0
m, the proposed sensor structure achieves a narrowed sensing scope from 1.26 to 1.32. As seen in
Figure 12c,f,i, sensing scope is from 1.26 to 1.34 when
L is set as 5.2
m. When
L is set as 5.6
m, the proposed sensor structure achieves a narrowed sensing coverage from 1.26 to 1.33. When
L is set as 6.0
m, the sensing scope of proposed sensor structure is further narrowed from 1.26 to 1.32. Reasons for these changes are considered that when value of
L becomes larger, the amount of coating increases and area A is significantly compressed. Thus the intensity of the envelope demonstrates a trend from increase to decrease. Also, the increasing value of
t changes the amount of Au coating too. Moreover, the variation of
L influences mode field distribution at the right side of area A in proposed sensor structure, which has an effect on the phase-matching condition of the coupled mode of area A in
x-polarization and causes the
curve and peak to demonstrate a visible change. According to calculation in Equation (4), the maximum values of
are achieved as 18,500, 29,000, 21,000, 22,000, 32,000, 25,000, 26,000, 18,000, 31,000 nm/RIU for scenarios in
Figure 12a–i, respectively. Moreover, the average
values achieved are 10,642, 10,714, 10,125, 12,500, 12,000, 10,857, 14,600, 12,833, 12,667 nm/RIU for scenarios in
Figure 12a–i, respectively. When the maximum
value is given priority to be taken into consideration, values of
L and
t are set as 5.6
m and 60 nm when the proposed sensor structure achieves an optimal sensing performance.
Figure 13 shows the
curve of the coupled mode in
x-polarization when
is selected as 0.1, 0.2, 0.3
m, and RI is chosen as 1.26, 1.31, 1.32 and 1.33. From
Figure 13, with an increasing value of
from 0.1 to 0.3
m, the
curve of the mode with RI selected as 1.26 almost does not shift and the curve intensity demonstrates a dramatic decrease. The
curve of the mode demonstrates a small red shift and its intensity first increases and then decreases when RI is chosen as the maximum value of RI sensing coverage for all chosen
values. In addition, with
increases, the envelope of the
curves associated with varying RI demonstrates a visible red shift and becomes narrow along with a decrease in envelope intensity. Additionally, the sensing scope is extended to 1.33 when
increases to 0.3
m. The reason for these changes is considered to be that when value of
becomes larger, the channel between core area A and the coating surface is narrowed, which influences the mode field distribution. Also, the variation of
changes the area for the SPR effect, which has a visible influence on the phase-matching condition of the coupled mode of
x-polarized area A and further causes an obvious change in the
curve along with the appearance of a new sensing RI 1.33. Moreover, with the channel narrowed, the core area is further isolated, which affects the mode couple. Thus, the intensity of the envelope illustrates a trend from ascending to descending. According to calculation in Equation (4), the maximum
values are achieved as 26,200, 32,000, 27,000 nm/RIU when
is selected as 0.1, 0.2, 0.3
m, respectively. It is clear from the calculation results that the maximum value of
in the proposed design has a trend from rise to decline with a variation of
from 0.1 to 0.3
m. Moreover, the average value of
in the proposed design decreases with a changing sequence of
. In addition, the average
values achieved are 12,167, 12,000, 11,143 nm/RIU when
is selected as 0.1, 0.2, 0.3
m, respectively. When the maximum
value and fabrication difficulty are considered priorities, the value of
is finally set as 0.2
m, resulting in the proposed sensor structure achieving an optimal sensing performance.
Figure 14 illustrates the
curve of the coupled mode in
x-polarization when
a is set as 0.7, 0.9, 1.1
m, b is fixed at 0.4
m, and RI is chosen as 1.26, 1.31, 1.32. From
Figure 14, with an increasing value of
a from 0.7 to 1.1
m when b is fixed at 0.4
m, the
curve starts to red-shift and the related curve intensity demonstrates a slight increase when RI is set as 1.26. The
curve when RI is set as 1.31 has a noticeable red-shift and its intensity first rises and then drops. The
curve with RI chosen as 1.32 demonstrates a visible red-shift and its intensity demonstrates a dramatic decrease. In addition, with the value of
a increasing, the envelope of the
curves associated with varying RI demonstrates a red shift and appears to become broader. Moreover, the intensity of the envelope shows a trend from climbing to falling when
a increases. The reason for these changes are considered to be that when the value of
a becomes larger, the core region area A is gradually horizontally compressed, influencing the mode field distribution. This further alters the phase-matching condition of the
x-polarized coupled mode of area A and further causes the slight shift of the
curve and peak. Additionally, the variations in the core area alter the concentration of energy, causing a change in the envelope intensity. According to calculation in Equation (4), the maximum
values are achieved as 30,500, 32,000, 30,500 nm/RIU when
a is chosen as 0.7, 0.9, 1.1
m, respectively. It is obvious from the calculation results that the maximum value of
in the proposed design shows a trend from increase to decrease with a variation of
a from 0.7 to 1.1
m. Additionally, the calculated average
value of the proposed design increases when
a changes from 0.7 to 1.1
m. Moreover, the average
values achieved are 11,500, 12,000, 12,667 nm/RIU when
a is selected as 0.7, 0.9, 1.1
m, respectively. Taking both
values and fabrication difficulties into account, the value of
a is finally set as 0.9
m when the proposed sensor structure achieves an optimal sensing performance.
Figure 15 demonstrates the
curve of the coupled mode in
x-polarization when
b is set as 0.2, 0.4, 0.6
m, a is fixed at 0.9
m, and RI is chosen as 1.26, 1.31, 1.32. From
Figure 15, when
b increases from 0.2 to 0.6
m when a is fixed at 0.9
m, the
curve starts to red-shift and the related curve intensity demonstrates a slight increase when RI has a variation sequence from 1.26 to 1.31. Also, when RI changes from 1.31 to 1.32, the
curve of this design has a significant red-shift and the intensity decreases. In addition, with an increasing value of
b, the envelope of
curves associated with varying RI values demonstrates a red shift and becomes broader. The reason for these changes are considered to be that when value of
b becomes larger, core region area A is vertically compressed to a certain degree, influencing the mode field distributions. This also alters the phase-matching condition of the
x-polarized coupled mode of area A and further causes the slight shift in the
curve and peak. Meanwhile, the change of the core area in the vertical direction has little influence on the concentration of the energy. Thus, the change of envelope intensity is small. According to calculation in Equation (4), maximum
values are achieved as 31,000, 32,000, 31,000 nm/RIU when
b is chosen as 0.2, 0.4, 0.6
m, respectively. It is clear from the calculation results that the maximum value of
in the proposed design illustrates a trend from rise to decline with a variation of
b from 0.2 to 0.6
m. Additionally, the calculated average
value of the proposed structure increases with
b changing from 0.2 to 0.6
m. Moreover, the average
values achieved are 11,333, 12,000, 12,167 nm/RIU when
b is chosen as 0.2, 0.4, 0.6
m, respectively. When maximum
value is considered a priority, the value of
b is finally chosen as 0.4
m when the proposed sensor structure achieves an optimal sensing performance.
For the design in this correspondence, a sensitivity analysis of the single structural parameter to sensing performances of the proposed design has been explored and summarized in
Table 2, as below.