2.1. Design and Characterization of the Meta-Grating
Figure 1(a) schematically depicts the geometry of the meta-grating we designed, which is composed of an array of silver nanoshells and a continuous gold film separated by a photoresist (PR) grating mask that is not developed to the end. The gold film prohibits any transmission through the system so that only the reflection of the device needs to be considered, laying the basis for generating high-contrast RA diffraction [
26]. Near-filed coupling between two metallic layers (i.e. Ag nanoshells and the gold film) forms a series of LSPR[
27,
28] in different wavelengths with their features dictated by geometrical details as defined in
Figure 1(a). The structural parameters of the MIM meta-grating are defined as follow: the thickness of the Au film
d1; the thickness of the photoresist spacer
d2; the modulation depth of the PR grating
h; the full width at half maximum of the PR grating lines
W; the period of the nanocavity array
P.
Figure 1(b) and 1(c) depict, respectively, the scanning-electron-microscopic (SEM) pictures of part of the fabricated sample before and after the deposition of Ag nanoshells on surface of the PR grating. Refer to [
15,
29] for the detailed fabrication process of such MIM meta-grating. Through parameter sweeps on the geometrical parameters, such as the full width at half maximum of the PR grating lines, the thickness of photoresist spacer and Ag nanoshells, using finite-difference time-domain (FDTD) simulations, we optimized the geometric parameters of the MIM mete-grating in the experiments as follow:
P≈465 nm,
d1=70 nm,
d2= 50 nm,
h=260 nm,
W=165 nm.
We experimentally characterize the optical properties of the fabricated samples depicted in Figures 1(b) and 1(c) for the sake of comparison. Red circles in Figures 1(d) and 1(e) represent the measured reflectance spectra of the samples under the illuminations of normally incident lights for TM polarization (perpendicular to the extending direction of grating lines), which are in nice agreement with corresponding FDTD simulations as plotted by blue dotted curves in Figures 1(d) and 1(e), respectively. In our simulations, geometric parameters of the simulated nanostructures are derived from the experiments. Meanwhile, the intrinsic absorptions of metallic materials are also considered in the simulation. For the dielectric grating without Ag nanoshells, the RA resonances in two spectra (i. e. the measured and simulated reflectance spectra of the grating) as marked by black dotted circular rings in
Figure 1(d) are observed as two extremely shallow dips at 463 nm for their natural sub-radiation properties. The linewidth of RA mode in FWHM is about 8 nm. The other two dips appear at 645 nm and 740 nm in two spectra, respectively, corresponding to SPP modes bound to different interfaces (i.e., the interface between Au film and PR spacer layer, and the interface between Au film and the SiO2 substrate) (see Figure S1 in the supplementary information for the corresponding electric-intensity distribution). Those spectral features can also be confirmed by the corresponding angle-resolved reflectance spectrum in below.
As to the MIM meta-grating, the RA resonances in two spectra as highlighted by the black and red dotted circular ring, respectively in
Figure 1(e) turn to weak peaks due to the relatively weak Fano couplings between RA and the broadband plasmonic modes with their central wavelength at
Near-field interaction between Ag nanoshells and gold film forms such broadband plasmonic modes with stronger absorption and electric-field enhancement (~150 nm in FWHM and a reflectivity below 20% at central wavelength), which lays the basis for the realization of wide range tunable Fano resonances. Compared to the simple RA mode mentioned in
Figure 1(e), the linewidth of the hybrid RA modes in FWHM is increased to 12 nm in this case. To quantitatively evaluate the optical properties of RA, the plasmonic modes and the coupling behavers between them, we applied the coupled-mode theory (CMT) mentioned in [
15,
30] to help us precisely extract the basic parameters of the resonance modes, including the intrinsic damping, the radiation damping, and the coupling strength between such two resonance modes, where the system is simplified as a one-port (reflection geometry) two-mode (considering only the RA and LSPR modes in this work) model. The black solid curve in
Figure 1(e) shows the CMT fitting spectrum with parameters as follow: radiation damping of RA and the broadband LSPR:
ΓRAe=5,
ΓLSPRe=25; intrinsic damping of such two modes:
ΓRAi=0.5,
ΓLSPRi=55; coupling strength between such two resonance modes:
t=5. Good agreement is observed between the measurement and the CMT fitting. In this case (i.e.
), the radiation ability to far field and the linewidth of the hybrid RA mode are slightly influenced by the relatively weak coupling process.
While, it’s worth noting that the relatively low reflectance of the MIM meta-grating in the wavelength range of 400-550 nm actually results from the intrinsic absorption of the gold film, which can be verified by the measured reflectance spectra of gold film as depicted in
Figures 1(f) and plotted by blue circles. Therefore, the spectral superposition of the intrinsic absorption of gold film and the broadband LSPR at
can perfectly reproduce the relatively complex spectral lineshape of MIM meta-grating in
Figure 1(e), which can be seen from the corresponding CMT fitting results as depicted in Figures 1(f) and plotted by blue (the intrinsic absorption of gold film) and red (the broadband LSPR) solid curves, respectively. Meanwhile, the reflectance peaks appear at 720 nm, corresponding to the Fano resonances between the broadband LSPR and the SPP mode (see Figure S2 in the supplementary information for the corresponding electric-intensity distribution), which can also be demonstrated by the following angle-resolved reflectance spectra.
Figure 2(g) shows the electric-intensity distribution of the MIM meta-grating at
and
, corresponding to the RA and LSPR modes, respectively, where the enhanced electric fields are mainly distributed on surface of the MIM grating and the interspace between two Ag nanoshells, respectively. Own to the weak Fano coupling between LSPR and RA, the near-field distribution of RA is slightly modulated by the LSPR, inducing relatively weak electric field localized in the interspace between two Ag nanoshells.
Figures 2(a) and 2(b) show the angle-resolved reflectance spectra of the fabricated samples depicted in Figures 1(b) and 1(c), respectively, with the incident angles turned from 0° to 10°, where the spectral shift of RA resonances are highlighted by series of dotted circular rings. Three obvious features near the RA diffracted wavelengths can be observed: (1) The diffracted wavelengths of RA in experiments are basically consistent with the theoretical calculations as plotted by black dotted lines in Figures 2(a) and 2(b). The central wavelength of RA is determined by the diffraction formula of
[
24], where
is the refractive index of the environment and
equals 1 when the grating is located in air,
is the grating period, and λ is the diffracted wavelength of RA. Such an angle-resolved behaver further verify their identities of RA diffractions; (2) The RA resonances in dielectric grating are indeed reflectance dips (see
Figure 2(a)), while they are reflectance peaks in MIM meta-grating (see
Figure 2(b)), except for the resonances as marked by red dotted circular rings where the Fano resonances between RA modes and the broadband plasmonic resonances are disturbed to a certain extent by the intrinsic absorption of gold film; (3) The coupling strengths between RA and the broadband plasmonic mode can be tuned by the incident angles (see
Figure 2(b)). Disturbed by the intrinsic absorption of gold film, the coupling strengths have no evident changes at small incident angles (i. e.
), indicated by weak radiation ability to far field of RA resonances. While, as the incident angle increased to
, an enhanced reflectance peak near RA mode can be recognized from the spectrum, implying the enhancement of the coupling strengths between such two modes. Meanwhile, we also note other angle-resolved features in Figures 2(a) and 2(b) at the wavelength range of 600-750 nm, which are marked by blue arrows and shallow-blue dotted lines. The angle-resolved features for those resonance modes are basically consistent with that of SPPs, thus, further verify their identities.
2.2. High-Contrast Fano Resonance and Its Application in Reflective Filters
Continuing to increase the angles of incident light, we observe gradually enhanced Fano couplings between RA and the broadband LSPR.
Figure 3(a) shows the corresponding angle-resolved reflectance spectra of the MIM meta-grating depicted in
Figure 1 (c) with the incident angles turned from 10° to 26°, where the meta-grating is labeled as “Samp. A”. In the angle-resolved spectra, we observe a series of narrowband and high-contrast reflection peaks turned from 567 nm to 667 nm, and with relatively high reflectance changed from 75% to 88%, which imply the gradually enhanced coupling strength between RA and the broadband LSPR. The reflectance reaches to the peak as the RA mode shift to the central wavelength of LSPR (i.e. λ=610 nm ) at the incident angle of
.
Figure 3(b) shows the corresponding measured and simulated reflectance spectra of Samp. A at
, which are plotted by red circles and blue dotted curves, respectively. The reflectance of the hybrid RA mode is up to 88%. While, as to strong Fano coupling between RA and the broadband LSPR, the linewidth of the hybrid RA mode in FWHM is increased to about 20 nm. The black solid curve in
Figure 3(a) shows the CMT fitting spectrum with parameters as follow: radiation damping of RA and broadband LSPR:
ΓRAe=10,
ΓLSPRe=25; intrinsic damping of such two modes:
ΓRAi=0.5,
ΓLSPRi=55; coupling strength between such two resonance modes:
t=28. Good agreement can be observed between the measurement and the CMT fitting except for the region of intrinsic absorption of gold film. Obviously, stronger coupling between such two modes enhances the radiation ability of RA mode to a certain degree (
ΓRAe turn from 5 to10 as the incident angles increased from
to
), while, it also increases the linewidth of the hybrid mode correspondingly, which is actually not appreciated in biological/chemical sensing applications.
The peak reflectance location of the hybrid RA can be turned by changing the central wavelength of LSPR, which can be verified by the angle-resolved reflectance spectra of the MIM meta-grating labeled as “Samp. B” as shown in
Figure 3(c). The central wavelength of LSPR in “Samp. B” is about 630 nm. Therefore, the reflectance reaches to about 100% near the central wavelength of LSPR at the incident angles of
and
. Obviously, the hybrid RA in “Samp. B” exhibit higher reflectance (i.e., 76%-100% as the incident angles turned from 13° to 33°). Compared to “Samp. A”, the thickness of the photoresist spacer
d2 and the modulation depth of the PR grating
h of “Samp. B” are turned to 55 nm and 300 nm, respectively, while, other geometric parameters keep unchanged. Such superior spectral performances of such MIM meta-grating make it promising for the application of high-efficiency and narrowband reflective filters. The diffracted wavelengths of the hybrid RA mode in Figures 3(a) and 3(c) are basically consistent with the theoretical calculation, which are demonstrated in
Figure 3(d) and plotted by red (Samp. A) or blue (Samp. B) circles (measurement results) and black solid line (theoretical calculation), respectively.
2.3. Sensor Performances of the MIM Meta-Grating
Benefit both from extremely high RI sensitivity and sharply asymmetric lineshape of RA diffraction, the FOM of such hybrid sensor based on Fano resonance between RA and the broadband LSPR can be described specifically as:
where the subscript “RA” denotes the hybrid RA mode. Based on above results, we know that for the MIM metagrating with defined period
, larger incident angle would introduce two inverse effects. On the one hand, it can effectively increases the sensitivity
of the hybrid sensors, on the other hand, it broadens the linewidth of the RA mode due to the strong coupling between RA and the broadband LSPR, and thereby worsening the sensors performances. As for such contradiction, we should find a balance between the sensitivity and linewidth of the resonance mode in such hybrid sensor by carefully controlling the angles of incident light.
To demonstrate above analysis, we simulate the reflectance spectra of the hybrid sensor immersing in different RI environments for incident angles of
=
,
and
, and depict them in Figures 4(a)-(c), respectively. The environmental RI are changed from 1 to 1.02 with a small step size of 0.005. It should be mentioned that “Samp. A” is applied in this section, while, the intrinsic absorption of gold is ignored to avoid its influence on sensor performances. As the incident angle is
, relatively broadband hybrid RA mode due to the degeneracy of
order RA modes appears at about
in
Figure 4(a) [
18]. The sensitivity of the hybrid sensor is about
, which shows good agreement with the theoretical expectation. As the incident angles increase to 10
and 20
, the sensitivities of the hybrid sensor increase to about
and
, respectively. The slight deviation between simulation results and theoretical calculations (i. e.,
and
, respectively.) may be due to angular dispersion of the reflectance spectra at large incident angles with relatively broadband source in FDTD simulations.
Figure 4(d) summarizes the simulated results and re-plots the wavelength shifts of sensor signal as a function of the environmental RI at incident angles of 𝜃=
,
and
with different colour dotted lines and decorated with different colour circles. It is obvious that the sensitivity of the hybrid sensors at 𝜃=
is the highest. While, affected by the relatively broad linewidths of the hybrid RA, the FOM of the sensor at 𝜃=
is not the highest. The simulated linewidths of the hybrid RA at the incident angles of
,
and
are 22 nm, 9 nm, 14 nm, respectively. Correspondingly, the FOM of those sensors are 22, 60 and 44.5. Therefore, the hybrid sensor performs better overall or has higher sensory quality when working at incidence angles of around 10
, which can also be demonstrated in the following sensor performances.
Figure 5(a)-5(c) show the reflectance spectra of the hybrid sensor which is covered uniformly by different thickness of thin film with RI
n=1.5 for incident angles of 𝜃=
,
and
, where the thickness of the coating layer are changed from 0 to 8 nm with a step size of 2 nm. Such a scheme is applied to mimic the local RI changes caused by the binding of target molecules on surface of the nano-devices. Through parameter sweeping on the incident angles, we ultimately choose such three typical sensor spectra at 𝜃=
,
and
. Obviously, the hybrid sensor worked at 𝜃=
is really the best.
Figure 5(d) summarizes the simulated results and re-plots the wavelength shifts of sensor signal as a function of the thickness of coating film at incident angles of 𝜃=
,
and
with different colour dotted lines and decorated with different colour circles. Meanwhile, we compare the performance of the hybrid sensor with the simple RA based sensor worked at 𝜃=
and the broadband plasmonic sensor worked at 𝜃=
for the sake of fairnessas, which are plotted in
Figure 5(d) by different colour dotted line and decorated with yellow and purple pentagonal shapes, respectively. According to the figure of merit for thin layer detecting (i.e.
with
denoting the thickness of the coating film) introduced by J. Becker et al.[
5], the average
of the hybrid sensor worked at 𝜃=
, the simple RA worked at 𝜃=
and the LSPR based sensor are about 0.139, 0.024 and 0.0092, respectively. Obviously, the hybrid sensor worked at 𝜃=
does really perform the best.
Furthermore, simulated results in
Figure 5(d) help us understand some more intriguing and important phenomenon. First, the spectral shift now depends on the relative thickness of the molecules to the volume of the resonance modes field penetrates into the medium. Due to strong electric-field enhancement of the LSPR and hybrid RA mode in nanoscale, the sensitivity of those sensors decreases slightly with the increase of film thickness. Therefore, the spectral shift of those sensors as a function of the thickness of thin film are non-linear distributions. Second, LSPR based sensor is the most sensitive to the local RI changes caused by the bonding of thin film, while, affected by the broadband lineshape, the
of the LSPR sensor is the lowest. Simple RA based sensor has narrowband linshape, while, it is not sensitive to local RI changes due to RA relatively large mode volume and weak photon state density. Finally, own to both the relatively narrowband lineshape of hybrid RA and its ability to modulate the enhanced electric-intensity distribution near the nanostructure, the FOM of hybrid sensor in detection of local RI changes is enhanced to 6 and 15 times than the simple RA and LSPR based sensors, respectively.