**1. Introduction**

Over the past two decades, there has been extensive research on plasmonic devices due to the capability to break the optical diffraction limit and confine light in sizes that are much smaller than the diffraction limit, thus enhancing the electric field [1,2]. Because of the surface wave nature, this type of device is usually very sensitive to surrounding refractive index changes [3–8], making them potential candidates for sensing application. A recent review paper summarized the state of the art in this area [9]. However, the intrinsic loss from the materials limits the linewidth of the device resonance and hinders further development. For example, interacting between impinging light and lossy materials will generate a thermal issue that modifies the surrounding medium and changes the sensitivity. One possible way to address this problem is to use BICs to mitigate the issue [10]. A figure-of-merit (FOM), defined as S/FWHM (S denotes the sensitivity, and FWHM denotes the full width at half maximum), is often adopted to characterize a sensor performance. An ideal sensor would possess a larger FOM, namely high sensitivity and a small FWHM. Therefore, an approach to reduce the FWHM of plasmonic sensors while maintaining their sensitivity is highly desired.

BIC, which is referred to as embedded eigenvalues or embedded trapped modes, was first proposed by von Neumann and Wigner in 1929 [11]. It is a localized state of an open structure with access to radiation channels, yet it remains highly confined with, in theory, an infinite lifetime and Q-factor. Therefore, it is considered an important approach for designing high-Q optical resonators [12–15] and has been applied to lasers [16–19] and sensors [20–22]. Recently, Azzam et al. [23] proposed a hybrid plasmonic-photonic structure that possesses two types of BICs, namely symmetry-protected BICs (SP-BICs)

**Citation:** Tang, S.; Chang, C.; Zhou, P.; Zou, Y. Numerical Study on a Bound State in the Continuum Assisted Plasmonic Refractive Index Sensor. *Photonics* **2022**, *9*, 224. https://doi.org/10.3390/ photonics9040224

Received: 28 January 2022 Accepted: 2 March 2022 Published: 28 March 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Copyright:** © 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

and Friedrich–Wintgen BICs (FW-BICs), of which the former is caused by the symmetry incompatibility between the local modes and the continuum spectrum and the latter is induced by the destructive interference of two scattering channels. Friedrich–Wintgen BICs (FW-BICs), of which the former is caused by the symmetry incompatibility between the local modes and the continuum spectrum and the latter is induced by the destructive interference of two scattering channels.

sensors [20–22]. Recently, Azzam et al. [23] proposed a hybrid plasmonic-photonic structure that possesses two types of BICs, namely symmetry-protected BICs (SP-BICs) and

In this paper, we numerically study a BIC-assisted plasmonic sensor that leverages the BIC effect to suppress the linewidth of the resonance. A quasi-SP-BIC mode at 1559 nm can be initiated by slightly tilting the incident angle to break the symmetry, featuring a narrow FWHM of 1.8 nm. The device possesses high sensitivity (~938 nm/RIU) (RIU denotes refractive index unit) and high FOM (~521/RIU) characteristics in refractive index sensing. By tuning the operating wavelength to 1905 nm, our device exhibits a higher sensitivity of ~1264 nm/RIU, a narrower FWHM of 0.18 nm, and a better FOM of ~7022/RIU. The devices also feature a 0.1 nm/(g/L) sensitivity for glucose solution detection at 1559 nm and a 0.14 nm/(g/L) sensitivity at 1905 nm. In this paper, we numerically study a BIC-assisted plasmonic sensor that leverages the BIC effect to suppress the linewidth of the resonance. A quasi-SP-BIC mode at 1559 nm can be initiated by slightly tilting the incident angle to break the symmetry, featuring a narrow FWHM of 1.8 nm. The device possesses high sensitivity (~938 nm/RIU) (RIU denotes refractive index unit) and high FOM (~521/RIU) characteristics in refractive index sensing. By tuning the operating wavelength to 1905 nm, our device exhibits a higher sensitivity of ~1264 nm/RIU, a narrower FWHM of 0.18 nm, and a better FOM of ~7022/RIU. The devices also feature a 0.1 nm/(g/L) sensitivity for glucose solution detection at 1559 nm and a 0.14 nm/(g/L) sensitivity at 1905 nm.

#### **2. Design of a BIC-Assisted Plasmonic Sensor 2. Design of a BIC-Assisted Plasmonic Sensor**

*Photonics* **2022**, *9*, x FOR PEER REVIEW 2 of 11

The schematic of the proposed device is shown in Figure 1a. It consists of a onedimensional aluminum nitride (AlN) grating on a metal–glass substrate. The momentum matching condition between the excited resonance mode and the incident light is [24]: The schematic of the proposed device is shown in Figure 1a. It consists of a one-dimensional aluminum nitride (AlN) grating on a metal–glass substrate. The momentum matching condition between the excited resonance mode and the incident light is [24]:

$$k\_{sp} = \frac{2\pi}{\lambda} n\_i \sin \theta + m \frac{2\pi}{a} = \frac{2\pi}{\lambda} \sqrt{\frac{\varepsilon\_m n\_s^2}{\varepsilon\_m + n\_s^2}} \tag{1}$$

where *ksp* is the momentum of the surface plasmon, *λ* is the free-space wavelength, *n<sup>s</sup>* is the refractive index of the surrounding medium, *ε<sup>m</sup>* is the dielectric constant of the metal layer, *θ* is the angle of incidence, *a* is the grating pitch, and *m* is the grating diffraction order. the refractive index of the surrounding medium, is the dielectric constant of the metal layer, is the angle of incidence, *a* is the grating pitch, and *m* is the grating diffraction order.

**Figure 1.** (**a**) Schematic of the proposed device, consisting of one-dimensional aluminum nitride (AlN) grating on Ag-SiO2 substrate. (**b**) Cross-section of the device. Here, the width and height of the grating are w = 320 nm and d = 440 nm, respectively, the period is a = 1070 nm, and the thickness of the silver film is h = 100 nm. **Figure 1.** (**a**) Schematic of the proposed device, consisting of one-dimensional aluminum nitride (AlN) grating on Ag-SiO2 substrate. (**b**) Cross-section of the device. Here, the width and height of the grating are *w* = 320 nm and *d* = 440 nm, respectively, the period is *a* = 1070 nm, and the thickness of the silver film is *h* = 100 nm.

In our device, we choose silicon dioxide as the substrate and silver as the metal layer for its low-loss property in the near-infrared range among all the metals [25]. We pick AlN as the grating material for its transparency in the near-infrared, which can effectively reduce the device loss and lead to a narrower resonance peak. For our simulation, we fix the refractive index of silicon dioxide at 1.45 for its negligible variations. The refractive indices of silver and AlN are from Refs [26,27], respectively. We set water as the surrounding medium for sensor design and fix its refractive index at 1.33. Figure 1b shows the crosssection of the device. The AlN grating with the width, *w* = 320 nm, is periodically arranged in the x-direction with a period, *a* = 1070 nm. The thicknesses of the grating, *d*, and silver film, *h*, are 440 nm and 100 nm, respectively. The value of the geometric parameter shown in Figure 1b is listed in Table 1.

**Table 1.** Value of the Geometric Parameter Shown in Figure 1b.


We employ the finite difference time domain (FDTD) method to calculate the reflection of the proposed device. Due to the invariance in the z-direction of our device, for convenience, 2D-FDTD is chosen. We simulate the device using a non-uniform grid with a minimum grid size of 1 nm. The perfectly matched boundary condition and periodic Bloch boundary condition are selected for the y-direction and the x-direction, respectively. We use TM polarized light as the incident light to excite the plasmonic mode at the metal-dielectric interface.

The calculated reflected spectra for the wavelength range of 1350–1650 nm under different incident angles are plotted in Figure 2a. Two counter-propagating plasmonic modes, SP [+1] and SP [−1], form standing waves, resulting in the formation of symmetric and anti-symmetric modes with different energy bands as shown in the insets of Figure 2b. The bandgap between the two energy bands covers the wavelength region of 1475–1540 nm, and the width indicates the strength of the interaction between the two modes. It is worth noting that, in the upper band, there is a discontinuous region, circled by the red dash line, at the wavelength ~1559 nm when the light is normally incident, which does not exist in the lower band. To better illustrate this phenomenon, we plot the reflected spectral when the incident angles are 0◦ and 2◦ in Figure 2b. Compared with the normal incidence, a clear Fano resonance, marked by a red dashed square, is observed around 1559 nm when the device is illuminated with a 2◦ incident angle, indicating an SP-BIC exists in the device. The corresponding magnetic field distribution is shown in the inset of Figure 2b, where the magnetic field is antisymmetric to the central line. Since the incident plane wave is symmetric when the wavevector is normal to the plane. Additionally, the reflected wave and transmitted wave are supposed to be plane waves with no diffraction order other than the 0th order due to the small lattice constant. The different types of symmetry prevent coupling between the SP-BIC mode and the output wave, granting the strong localization of this mode. Since the incident wave cannot couple to the other modes at this wavelength, most of the energy is reflected, resulting in the high reflective coefficient, which explains the discontinuity of the upper band in Figure 2a. To excite a quasi-SP-BIC mode, we can slightly tilt the incident angle to break the symmetry and provide a certain amount of mode overlap with the quasi-SP-BIC mode.

**Figure 2.** (**a**) The reflected spectra for different incident angles. A discontinuity is observed at the center of the red circle, indicating the existence of an SP-BIC. (**b**) The reflection spectrum when the incident angles are 0° and 2°, respectively. The insets show the magnetic field distribution for the corresponding resonances circled by the red and blue lines in Figure 2a. (**c**) The reflection spectral for different incident angles. The inset shows the value of the CR/FWHM for different incident angles. **Figure 2.** (**a**) The reflected spectra for different incident angles. A discontinuity is observed at thecenter of the red circle, indicating the existence of an SP-BIC. (**b**) The reflection spectrum when theincident angles are 0◦ and 2◦ , respectively. The insets show the magnetic field distribution for the corresponding resonances circled by the red and blue lines in Figure 2a. (**c**) The reflection spectral for different incident angles. The inset shows the value of the CR/FWHM for different incident angles.

Compared with the anti-symmetric mode in the upper band, the lower band features symmetric mode distribution, as shown in the blue dashed square in Figure 2b. Due to the same symmetry between the incident light and the mode, even normal incidence can excite the resonance, which results in the continuity of the lower band in Figure 2a.

In Figure 2b, we also observe the different linewidths of the two modes. As seen from the insets, compared with the symmetric mode, the anti-symmetric mode is better confined in the device, which results in a reduced radiation loss. Therefore, the anti-symmetric mode (FWHM = 1.8 nm) has a narrower resonance peak than the symmetric one (FWHM = 5.7 nm). Although exciting the SP-BIC is hard, we can leverage the quasi-SP-BIC mode to suppress the resonance linewidth and significantly improve the device sensing performance.

We also study the property of the quasi-SP-BIC mode by gradually tilting the incident angle to deviate from zero-degree, as shown in Figure 2c. As the incident angle gradually increases, the resonance contrast ratio (CR), which refers to the reflection difference between the off-resonance and the on-resonance, becomes larger, indicating more energy couples into the quasi-SP-BIC mode due to increasing mode overlap. However, this process is associated with increased radiation loss due to the deviation from the SP-BIC condition, which, in turn, widens the FWHM. Therefore, we make a trade-off between the FWHM and the resonance CR by calculating the value of CR/FWHM, as shown in the inset of Figure 2c, and choose 2◦ as the final incident angle for its biggest value of CR/FWHM.

### **3. Results**

### *3.1. Sensitivity Characterization and Refractive Index Sensing*

To evaluate the device sensing performance, we simulate the response of the device under different environmental mediums with a refractive index range from 1.33 to 1.57. During the calculation, we fixed the grating period, height, and width to 1070 nm, 440 nm, and 320 nm, respectively, and the silver thickness to 100 nm. Here, we set the incident angle to 2◦ to slightly break its symmetry and enter the quasi-SP-BIC mode. Figure 3a shows multiple reflected spectra when the surrounding refractive index gradually increases. The quasi-SP-BIC mode red-shifts from 1559 nm to 1784 nm while the other resonance peak red-shifts from 1447 nm to 1647 nm. Linear fittings in Figure 3b indicate a linear resonance shift for both the two modes, of which the former has a sensitivity of ~938 nm/RIU and the latter's sensitivity is 833 nm/RIU. Considering the sharpness of the resonance, we obtain the FOM of the quasi-SP-BIC mode about 521/RIU, which is also the FOM of our device. In the following part, we will focus our discussion on the quasi-SP-BIC mode for its higher FOM.

Figure 3c presents the simulated resonance shift corresponding to different glucose concentrations in water with the incident angle at 2◦ . The refractive index of the glucose solution is taken from Ref. [28]. The resonance red-shifts linearly as the glucose concentration increases, showing 0.1 nm/(g/L).

### *3.2. Performance Change with the Thickness of the Silver Layer*

In our device, we choose silver to generate the plasmonic mode. Its thickness plays an important role in the device's performance. It has been proved that when the metal layer is thinner, parts of the energy are radiated into the substrate, resulting in a larger radiation loss, which gradually stabilizes when the thickness of the metal layer approaches the skin depth of the metal [29]. We calculate the reflected spectra of the quasi-SP-BIC under different silver layer thicknesses in Figure 4. During the calculation, all other parameters are fixed (*a* = 1070 nm, *w* = 320 nm, *d* = 440 nm, *n<sup>s</sup>* = 1.33, *θ* = 2◦ ). The thickness of the silver layer could slightly affect the position and the linewidth of the resonance. When the thickness changes from 60 nm to 120 nm with a step of 20 nm, the resonance wavelengths are 1558.64 nm, 1558.42 nm, 1558.56 nm, and 1558.47 nm, and the corresponding FWHM, shown in the inset of Figure 4, are 2 nm, 1.8 nm, 1.8 nm, and 1.8 nm. When the silver layer is thicker than 80 nm, the FWHM becomes stable, indicating no further substrate leakage occurs. In addition, the resonance CR stabilizes when the thickness is greater than 80 nm. Considering the FWHM and the resonance CR, a silver film with a thickness larger than

80 nm is preferred. We thus pick 100 nm as our silver thickness to provide a sufficient buffer for device implementation. symmetric mode distribution, as shown in the blue dashed square in Figure 2b. Due to the same symmetry between the incident light and the mode, even normal incidence can excite the resonance, which results in the continuity of the lower band in Figure 2a.

Compared with the anti-symmetric mode in the upper band, the lower band features

*Photonics* **2022**, *9*, x FOR PEER REVIEW 5 of 11

**Figure 3.** (**a**) The reflected spectra for different surrounding mediums. (**b**) The resonance wavelengths shift as a function of the surrounding index variation. The linear fittings show 833 nm/RIU and 938 nm/RIU for Mode1 and Mode2, respectively. (**c**) The reflected spectra for different concentrations of the glucose solution. The inset shows the resonance wavelength shifts with the glucose solution concentration variation. **Figure 3.** (**a**) The reflected spectra for different surrounding mediums. (**b**) The resonance wavelengths shift as a function of the surrounding index variation. The linear fittings show 833 nm/RIU and 938 nm/RIU for Mode1 and Mode2, respectively. (**c**) The reflected spectra for different concentrations of the glucose solution. The inset shows the resonance wavelength shifts with the glucose solution concentration variation.

Figure 3c presents the simulated resonance shift corresponding to different glucose concentrations in water with the incident angle at 2°. The refractive index of the glucose

ness to provide a sufficient buffer for device implementation.

tration increases, showing 0.1 nm/(g/L).

*3.2. Performance Change with the Thickness of the Silver Layer* 

**Figure 4.** The reflected spectra for different thicknesses of the silver layer. The inset shows the resonance CR and the FWHM change with the silver layer thickness. **Figure 4.** The reflected spectra for different thicknesses of the silver layer. The inset shows the resonance CR and the FWHM change with the silver layer thickness.

solution is taken from Ref. [28]. The resonance red-shifts linearly as the glucose concen-

In our device, we choose silver to generate the plasmonic mode. Its thickness plays an important role in the device's performance. It has been proved that when the metal layer is thinner, parts of the energy are radiated into the substrate, resulting in a larger radiation loss, which gradually stabilizes when the thickness of the metal layer approaches the skin depth of the metal [29]. We calculate the reflected spectra of the quasi-SP-BIC under different silver layer thicknesses in Figure 4. During the calculation, all other parameters are fixed (*a* = 1070 nm, *w* = 320 nm, *d* = 440 nm, = 1.33, = 2°). The thickness of the silver layer could slightly affect the position and the linewidth of the resonance. When the thickness changes from 60 nm to 120 nm with a step of 20 nm, the resonance wavelengths are 1558.64 nm, 1558.42 nm, 1558.56 nm, and 1558.47 nm, and the corresponding FWHM, shown in the inset of Figure 4, are 2 nm, 1.8 nm, 1.8 nm, and 1.8 nm. When the silver layer is thicker than 80 nm, the FWHM becomes stable, indicating no further substrate leakage occurs. In addition, the resonance CR stabilizes when the thickness is greater than 80 nm. Considering the FWHM and the resonance CR, a silver film with a thickness larger than 80 nm is preferred. We thus pick 100 nm as our silver thick-

### *3.3. Performance Change with the Grating Variables 3.3. Performance Change with the Grating Variables*

The AlN grating is also crucial for our device. There are multiple parameters in a grating, e.g., the grating height, width, and period, that would affect the device's performance. In this section, we provide a detailed analysis of these parameters on the perfor-The AlN grating is also crucial for our device. There are multiple parameters in a grating, e.g., the grating height, width, and period, that would affect the device's performance. In this section, we provide a detailed analysis of these parameters on the performance of the sensor. Note that in this section, the silver thickness is fixed at 100 nm.

mance of the sensor. Note that in this section, the silver thickness is fixed at 100 nm. First, we calculate the reflected spectra for different grating heights with fixed grating width (w = 320 nm), period (a = 1070 nm), and incident angle ( = 2°). Figure 5a shows that when the grating height gradually thickens from 400 nm to 480 nm, the resonance moves to the longer wavelength without significant distortion of its shape. The corresponding FWHM of the resonant peak is around 1.8 nm during the parameter scanning. We also do refractive index sensitivity characterization, showing the sensitivity is within the range of 920 nm/RIU to 960 nm/RIU. Therefore, a thinner grating possesses a higher sensitivity and thus a higher FOM, as shown in the inset of Figure 5a. However, from Figure 5a, we observe a slight decrease in the resonance CR for the grating height is smaller than 440 nm. Considering the overall performance, we pick 440 nm as our grating height to provide a First, we calculate the reflected spectra for different grating heights with fixed grating width (*w* = 320 nm), period (*a* = 1070 nm), and incident angle (*θ* = 2◦ ). Figure 5a shows that when the grating height gradually thickens from 400 nm to 480 nm, the resonance moves to the longer wavelength without significant distortion of its shape. The corresponding FWHM of the resonant peak is around 1.8 nm during the parameter scanning. We also do refractive index sensitivity characterization, showing the sensitivity is within the range of 920 nm/RIU to 960 nm/RIU. Therefore, a thinner grating possesses a higher sensitivity and thus a higher FOM, as shown in the inset of Figure 5a. However, from Figure 5a, we observe a slight decrease in the resonance CR for the grating height is smaller than 440 nm. Considering the overall performance, we pick 440 nm as our grating height to provide a higher FOM without sacrificing the CR of the resonance. This also provides a reference for our future design.

Next, we plot the reflected spectra for different grating widths with fixed grating height (*d* = 440 nm), period (*a* = 1070 nm), and incident angle (*θ* = 2◦ ) in Figure 5b. We find a red-shifted and slightly broadened resonant peak as we widen the grating width, i.e., increasing the duty cycle. The FWHM increases from 1.6 nm for the 280 nm wide grating to 2 nm for the 360 nm wide one, while the sensing sensitivity monotonically drops from 958 nm/RIU for the 280 nm wide grating to 923 nm/RIU for the 360 nm wide grating. The corresponding FOM also monotonically drops from 598/RIU to 461/RIU, as shown in the inset of Figure 5b, indicating a better sensing performance for the grating with a lower duty cycle. However, the inset of Figure 5b also shows an increasing resonance CR for wider grating teeth. Therefore, by taking a trade-off, we pick 320 nm as the grating width.

We also simulated the reflected spectra for different grating periods (from 970 nm to 1170 nm with a 50 nm step) and kept all other parameters fixed (*w* = 320 nm, *d* = 440 nm, and *θ* = 2◦ ). As shown in Figure 5c, when the period increases, the corresponding FWHM slightly drops from 1.9 nm to 1.7 nm, and the sensitivity increases from 826 nm/RIU to 1060 nm/RIU, providing a monotonically increasing FOM from 434/RIU to 623/RIU. In the meantime, we also notice a decreasing resonance CR as we increase the period. Therefore, we trade-off these factors and select the period *a* = 1070 nm for our design.

for our future design.

higher FOM without sacrificing the CR of the resonance. This also provides a reference

**Figure 5.** The reflected spectra for the grating with (**a**) different heights, (**b**) different widths, and (**c**) different grating periods. Insets show the corresponding FOM and the resonance CR versus the grating (**a**) height, (**b**) width, and (**c**) period. **Figure 5.** The reflected spectra for the grating with (**a**) different heights, (**b**) different widths, and (**c**) different grating periods. Insets show the corresponding FOM and the resonance CR versus the grating (**a**) height, (**b**) width, and (**c**) period.

### Next, we plot the reflected spectra for different grating widths with fixed grating *3.4. Performance Change with the Grating Shape Deviation*

height (*d* = 440 nm), period (*a* = 1070 nm), and incident angle ( = 2°) in Figure 5b. We find a red-shifted and slightly broadened resonant peak as we widen the grating width, i.e., increasing the duty cycle. The FWHM increases from 1.6 nm for the 280 nm wide grating to 2 nm for the 360 nm wide one, while the sensing sensitivity monotonically drops from We also study the device tolerance by considering the fabrication-induced sidewall angle deviations. In reality, it is common to have a trapezoid cross-section similar to the inset of Figure 6. In this section, we fix all the parameters same as in Section 2, namely *a* = 1070 nm, *d* = 440 nm, *w* = 320 nm, *h* = 100 nm, *θ* = 2◦ , and *n<sup>s</sup>* = 1.33, and scan the ∆w

958 nm/RIU for the 280 nm wide grating to 923 nm/RIU for the 360 nm wide grating. The

from −50 nm to 50 nm with a step of 25 nm. We plot the reflected spectra for different values of ∆w in Figure 6. As ∆w increases from −50 nm to 50 nm, similar to widening the grating width, the resonance shifts to the longer wavelength, and the CR gradually increases. The corresponding sensitivity decreases from 945 nm/RIU to 924 nm/RIU. During the whole process, the FWHM does not vary too much. It is fair to say that the device shows a sufficient tolerance for certain fabrication-induced geometric errors. ∆w in Figure.6. As ∆w increases from −50 nm to 50 nm, similar to widening the grating width, the resonance shifts to the longer wavelength, and the CR gradually increases. The corresponding sensitivity decreases from 945 nm/RIU to 924 nm/RIU. During the whole process, the FWHM does not vary too much. It is fair to say that the device shows a sufficient tolerance for certain fabrication-induced geometric errors.

We also study the device tolerance by considering the fabrication-induced sidewall angle deviations. In reality, it is common to have a trapezoid cross-section similar to the inset of Figure 6. In this section, we fix all the parameters same as in Sec II, namely *a* = 1070 nm, *d* = 440 nm, *w* = 320 nm, *h* = 100 nm, = 2°, and = 1.33, and scan the ∆w from −50 nm to 50 nm with a step of 25 nm. We plot the reflected spectra for different values of

we trade-off these factors and select the period *a* = 1070 nm for our design.

*3.4. Performance Change with the Grating Shape Deviation*

corresponding FOM also monotonically drops from 598/RIU to 461/RIU, as shown in the inset of Figure 5b, indicating a better sensing performance for the grating with a lower duty cycle. However, the inset of Figure 5b also shows an increasing resonance CR for wider grating teeth. Therefore, by taking a trade-off, we pick 320 nm as the grating width. We also simulated the reflected spectra for different grating periods (from 970 nm to 1170 nm with a 50 nm step) and kept all other parameters fixed (*w* = 320 nm, *d* = 440 nm, and = 2°). As shown in Figure 5c, when the period increases, the corresponding FWHM slightly drops from 1.9 nm to 1.7 nm, and the sensitivity increases from 826 nm/RIU to 1060 nm/RIU, providing a monotonically increasing FOM from 434/RIU to 623/RIU. In the meantime, we also notice a decreasing resonance CR as we increase the period. Therefore,

*Photonics* **2022**, *9*, x FOR PEER REVIEW 8 of 11

**Figure 6.** The reflected spectra for different ∆w. The inset shows the cross-section of the device. **Figure 6.** The reflected spectra for different ∆w. The inset shows the cross-section of the device.

### **4. Discussion 4. Discussion**

According to Equation (1), the resonance wavelength is proportional to the grating period. Further increasing the period may move the operating wavelength towards a longer wavelength. On the other hand, if we define the complex refractive index of silver According to Equation (1), the resonance wavelength is proportional to the grating period. Further increasing the period may move the operating wavelength towards a longer wavelength. On the other hand, if we define the complex refractive index of silver as *n<sup>m</sup>* = *n* + *ik*, for a plasmonic device, the quality factor is [30]:

$$Q\_{sp} = \frac{k^3}{2n} \tag{2}$$

(2)

A larger quality factor indicates a lower loss of a device. Take 1559 nm and 1905 nm as an example. Even though moving from 1559 nm to about 1905 nm increases the imaginary part of the complex refractive index, , from 11.439 to 13.862, it also raises the real part, *n*, from 0.1453 to 0.23119. Therefore, the indeed increases from about 5150 to about 5760, implying a lower device loss when working at a longer wavelength. In the A larger quality factor indicates a lower loss of a device. Take 1559 nm and 1905 nm as an example. Even though moving from 1559 nm to about 1905 nm increases the imaginary part of the complex refractive index, *k*, from 11.439 to 13.862, it also raises the real part, *n*, from 0.1453 to 0.23119. Therefore, the *Qsp* indeed increases from about 5150 to about 5760, implying a lower device loss when working at a longer wavelength. In the meantime, a resonance working at a longer wavelength shifts more than the one at a shorter wavelength.

meantime, a resonance working at a longer wavelength shifts more than the one at a shorter wavelength. We thus increase the period to shift the resonance peak of the quasi-SP-BIC from 1559 nm to a longer wavelength. Then, we optimize the device parameters as *a* = 1350 nm, *w* = 400 nm, *d* = 240 nm, and *h* = 100 nm. With incident angle *θ* = 2◦ , we apply glucose solution with different concentrations to the device, showing 1264 nm/RIU, corresponding to 0.14 nm/(g/L) for glucose solution, and 0.18 nm FWHM in Figure 7. The sensitivity is more than 30% higher than its counterpart at 1559 nm, while the FWHM is one order narrower than the 1559 nm one, leading to a much higher FOM of ~7022/RIU. We compare the performance of our devices with other recently reported plasmonic sensors in Table 2.

**Figure 7.** (**a**) The reflected spectra for different concentrations of the glucose solution. The inset shows resonance wavelength shifts with the glucose solution concentration variation. (**b**) The resonance wavelengths shift with the glucose solution concentration change-induced surrounding index variation. The linear fitting shows 1264 nm/RIU. **Figure 7.** (**a**) The reflected spectra for different concentrations of the glucose solution. The inset shows resonance wavelength shifts with the glucose solution concentration variation. (**b**) The resonance wavelengths shift with the glucose solution concentration change-induced surrounding index variation. The linear fitting shows 1264 nm/RIU.

We thus increase the period to shift the resonance peak of the quasi-SP-BIC from 1559 nm to a longer wavelength. Then, we optimize the device parameters as *a* = 1350 nm, *w* = 400 nm, *d* = 240 nm, and h = 100 nm. With incident angle = 2°, we apply glucose solution with different concentrations to the device, showing 1264 nm/RIU, corresponding to 0.14 nm/(g/L) for glucose solution, and 0.18 nm FWHM in Figure 7. The sensitivity is more than 30% higher than its counterpart at 1559 nm, while the FWHM is one order narrower than the 1559 nm one, leading to a much higher FOM of ~7022/RIU. We compare the performance of our devices with other recently reported plasmonic sensors in Table 2.


This paper @1905 nm 1264 0.18 7022


### **5. Conclusions**

> In summary, we propose a grating-assisted plasmonic refractive index sensor. Combining the quasi-SP-BIC effect and the surface wave nature, our device features a small FWHM and high sensitivity, resulting in FOM of 521/RIU at 1559 nm and 7022/RIU at 1905 nm. We also do glucose solution sensing, showing a 0.1 nm/(g/L) sensitivity at 1559 nm and a 0.14 nm/(g/L) sensitivity at 1905 nm. The performance of the proposed scheme is significantly better than the other grating-based plasmonic sensors reported earlier. Our approach provides a feasible way for designing a high-performance sensor with narrow FWHM and high sensitivity.

**Author Contributions:** Conceptualization, S.T. and P.Z.; methodology, S.T.; software, S.T. and C.C.; validation, S.T.; formal analysis, S.T.; investigation, S.T.; resources, S.T. and Y.Z.; data curation, S.T.; writing—original draft preparation, S.T.; writing—review and editing, Y.Z.; visualization, S.T.; supervision, Y.Z.; project administration, Y.Z.; funding acquisition, Y.Z. All authors have read and agreed to the published version of the manuscript.

**Funding:** The research was sponsored by the National Natural Science Foundation of China (NSFC) (61705099) and the Natural Science Foundation of Shanghai (21ZR1443100).

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Not applicable.

**Conflicts of Interest:** The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

### **References**

