2.1 Super-long range surface plasmon polaritons modes
Our SPP sensor consists of an Ag-Si-Ag slab and a flow cell filled with the detecting fluid on top of the waveguide, as shown in
Fig. 1
. The two Ag layers have the same thickness of 30 nm, with the waveguide Si layer sandwiched between them. The variable d
Si and n
D denote the Si layer thickness and the refractive index of the detecting fluid, respectively. The working wavelength is fixed at 1550 nm with the Si refractive index of 3.518 and the Ag permittivity of −86.6424 + 8.7422i [
13]. The SPP modes are propagating along the z axis with the device length of L. The sensor can be thought as an air-Ag-Si-Ag-fluid five-layer slab in x-z plane since it is thick enough in y axis.
Fig. 1 The geometry of the SPP sensor consisting of an Ag-Si-Ag slab and a flow cell filled with the detecting fluid.
Optical waveguide theory [
14D. Marcuse, Theory of dielectric optical waveguides (Academic, 1974).
] and transfer matrix method are applied to calculate the complex effective index of our five-layer slab, with its real part and imaginary part representing the mode propagation constant and losses, respectively. As SPP are TM waves which contain three electromagnetic component-H
y, E
x and E
z, we solve the scalar Helmholtz equation for TM mode in each layer as the following:
where,
k
0 = 2π/λ, it is the wave vector of vacuum and ε is the relative permittivity.
β represents the z-component propagation constant of the mode, and the effective index n
eff can be calculated by
β/
k
0. Once the H
y component is known, E
x and E
z can be calculated by the following equations:
In
formula (2), ω is the angular frequency. E
y, H
x and H
z components vanish for TM mode. By applying the continuity of the tangential components H
y and E
z in the boundary [
14D. Marcuse, Theory of dielectric optical waveguides (Academic, 1974).
], we obtain the relations between each adjacent layer. Then the total transfer matrix
M can be obtained via applying the transfer matrix method. By solving the equation |
M| = 0, the complex effective index is obtained.
We choose water first as our detecting fluid with refractive index of 1.33.
Figure 2
and
Fig. 3
illustrate the dispersion relations of our five-layer slab with different Si layer thicknesses ranging from 0 nm to 800 nm. The real part of the mode complex effective index, which we call the mode index below, is illustrated in
Fig. 2. This figure illustrates that the five-layer slab supports more TM modes as the Si layer thickness increases. As the Si layer thickness approaches 0 nm, TM
1 and TM
2 modes exhibit the mode index approaching water refractive index and air refractive index, respectively. These results suggest that an Ag/water interface SPP and an Ag/air interface SPP are formed, respectively. When the Si layer is thicker, higher order modes will form the Ag/water and the Ag/air interface SPP instead of TM
1 and TM
2 modes, with their mode indices approaching water refractive index and air refractive index, respectively. However, these SPP modes exhibit some differences between their fields in the Si layer.
Figures 4(a)
-
4(d) show the z direction normalized energy flows-S
z component field of four points marked in
Fig. 2 with a (TM
2 mode with d
Si = 300 nm), b (TM
4 mode with d
Si = 700 nm), c (TM
3 mode with d
Si = 300 nm) and d (TM
5 mode with d
Si = 700 nm), respectively. It is found from
Figs. 4(a) and
4(b) that the mode energy of points a and b are focused at the Ag/water interface, while some differences appear that the mode of point b exhibits more standing waves in the Si layer because of a thicker Si layer (point b with d
Si = 700 nm to point a with d
Si = 300 nm). The same phenomenon can also be found when comparing the field of point c to point d. These results show that the Ag/water interface SPP and Ag/air interface SPP formed by higher order modes exhibit more standing waves in the Si layer.
Fig. 2 The real part of the complex effective index with different Si layer thicknesses of the air-Ag-Si-Ag-water five-layer slab, both the Ag layer thicknesses are fixed at 30 nm. Three so-called interspace cut-off regions, which represent the discontinuous regions of the dispersion curves, are marked with A, B and C, respectively.
Fig. 3 The imaginary part of the complex effective index with different Si layer thicknesses of the air-Ag-Si-Ag-water five-layer slab, both the Ag layer thicknesses are fixed at 30 nm. Three so-called interspace cut-off regions, which represent the discontinuous regions of the dispersion curves, are marked with A, B and C, respectively.
Fig. 4 Normalized S
z field of (a) TM
2 mode with d
Si = 300 nm (b) TM
4 mode with d
Si = 700 nm (c) TM
3 mode with d
Si = 300 nm (d) TM
5 mode with d
Si = 700 nm in the air-Ag-Si-Ag-water five-layer slab. The four cases are also marked in
Fig. 2 as points a, b, c and d, respectively.
It is important to note that, as the Si layer thickness increases, an unique phenomenon appears in which the dispersion curves of higher order modes are non-continuous: TM
2 mode appears to be cut-off around the mode index close to the water refractive index where we name as the
interspace cut-off region (see region A in
Fig. 2 and
Fig. 3), while other higher order modes (e.g., TM
3 mode) have two interspace cut-off regions around the mode index close to the air refractive index and the water refractive index (see region B and C in
Fig. 2 and
Fig. 3), respectively. Here in these calculations, ‘interspace’ means the space where the dispersion curves are non-continuous, and we regard ‘cut-off’ as the region where the mode do not have solutions with positive imaginary effective index. But the mode may have solutions with negative imaginary part of effective index which require growing waves and we will discuss this more thoroughly later. To the best of our knowledge, the so-called interspace cut-off region has not been reported before. The possible origins of the interspace cut-off region lie in the thinness of the Ag layer, which we will also discuss later and we first put emphasis on the properties of the mode near the interspace cut-off region, which can be found in
Fig. 3. This figure shows the imaginary part of the mode effective index of our five-layer slab. Since TM
4 mode and other higher order modes have the same properties as TM
3 mode, they are not given in
Fig. 3, in which we find that the mode losses become infinitely small when the mode operates near the interspace cut-off region.
From the field’s point of view, all the modes near their interspace cut-off regions exhibit radiative fields into water since all their mode indices are lower than the water refractive index.
Figure 5(a)
shows the H
y component field of TM
3 mode with Si layer thickness of 401 nm, where the mode is near the interspace cut-off region C. Inset of
Fig. 5(a) shows the local field divided by three regions (air, Ag-Si-Ag and water) with two red dashed lines. The effective index of the mode is calculated to be 1.2485 + 0.0052i with a damping attenuation field of about 17.9 μm into water due to the losses of the mode. As TM
3 mode approaches the interspace cut-off region C, the mode losses approaches zero since more and more energy propagates in the lossless medium water, which is also the reason why we call this the super-long range SPP mode. The so-called super-long range SPP mode can exhibit infinitely small losses with finite Ag layer thickness, which is the same as the case of the LRSPP mode in an asymmetric IMI slab, and both the cases of lossless limit require an infinite plane wave in the higher index medium that compensates for dissipation in the metal [
11
J. J. Burke, G. I. Stegeman, and T. Tamir, “Surface-polariton-like waves guided by thin, lossy metal films,” Phys. Rev. B Condens. Matter
33(8), 5186–5201 (1986). [CrossRef]
[PubMed]
].
Fig. 5 (a) Normalized H
y field of TM
3 mode with Si layer thickness of 401 nm and Ag layer thickness of 30 nm in the air-Ag-Si-Ag-water five-layer slab. Inset shows the local field divided by three regions (air, Ag-Si-Ag and water) with two red dashed lines. (b) The real part of the complex effective index with different Si layer thicknesses, where the Ag layer thickness is changed to 70 nm with other parameters remained the same. (c) (d) Normalized S
z field of TM
3 mode with Si layer thickness of (c) 401.998 nm (d) 412.745 nm and Ag layer thickness of 30 nm. The two cases are marked in
Fig. 2 as point e and f, respectively.
Here we make another simulation to compare the super-long range SPP mode in our structure with the LRSPP mode in the asymmetric IMI slab calculated by Burke
et al. [
11
J. J. Burke, G. I. Stegeman, and T. Tamir, “Surface-polariton-like waves guided by thin, lossy metal films,” Phys. Rev. B Condens. Matter
33(8), 5186–5201 (1986). [CrossRef]
[PubMed]
]. In their work, they called the LRSPP mode as the Symmetric bound (S
b) mode. They found that there is a critical metal layer thickness h
C, above which the S
b mode corresponds to a bound wave. But when the metal layer thickness is below h
C, the bound wave no longer exists and instead of which, a growing wave solution appears which require both the amplitude grow with propagation distance and the fields in the higher index medium grow exponentially away from the metal layer. The growing wave solution exhibits a negative imaginary effective index which is only logical in mathematics. So we regard the region with growing wave solution as a cut-off case because it is not a real solution in physics. To compare with Burke
et al.’s work, we calculate the growing wave solutions around the interspace cut-off region C (see
Fig. 2) and the result is shown in
Fig. 6(a)
and
6(b). This Figure shows the real (
Fig. 6(a)) and imaginary (
Fig. 6(b)) part of effective index with Si layer thicknesses around interspace cut-off region C. Other than
Fig. 2 and
Fig. 3, the growing wave solutions which exhibit negative imaginary effective index are included. The red lines correspond to TM
3 mode, while the blue lines correspond to the growing wave solutions that take up the space of interspace cut-off region C to the moment. Two green points named as h
C1 and h
C2 represent the critical Si layer thicknesses in the boundary between TM
3 mode and the mode with growing wave solutions. Both our points h
C1, h
C2 and Burke
et al.’ point h
C represent the lossless limit that requires a plane wave in the higher index medium. So when comparing with Burke
et al.’ result, our design effectively squeezes the region with growing wave solutions to very small ranges of Si layer thicknesses. Again, since the growing wave solutions are not real solutions in physics, they are not shown in
Fig. 2 and
Fig. 3.
Fig. 6 (a) The real and (b) imaginary part of complex effective index with different Si layer thicknesses around interspace cut-off region C in
Fig. 2, both the Ag layer thicknesses are fixed at 30 nm. The growing wave solutions are also included. The red lines correspond to TM
3 mode, while the blue lines correspond to the growing wave solutions that take up the space of interspace cut-off region C to the moment. Two green points named as h
C1 and h
C2 represent the critical Si layer thicknesses in the boundary between TM
3 mode and the mode with growing wave solutions.
Then we discuss the possible cause of the interspace cut-off region. In
Fig. 2 and
Fig. 3, as the Si layer thickness increases near any interspace cut-off region, the focused energy will be transferred from the inner Si layer to the outer water. To better demonstrate this, we also choose TM
3 mode and interspace cut-off region C for instance. We first calculate the two critical interspace cut-off Si layer thicknesses both from the left side and the right side of interspace cut-off region C, and they are found to be 401.998 nm and 412.745 nm (see points e and f in
Fig. 2, also points h
C1 and h
C2 in
Fig. 6, precision of 0.001 nm), respectively. It means that TM
3 mode is interspace cut-off with the Si layer thickness ranging from 401.999 nm to 412.744 nm. The S
z field of points e and f are given in
Fig. 5(c) and
5(d), respectively. By comparing
Fig. 5(c) with 5(d), we find that more energy are focused in the inner Si layer in the field of point e than that of point f, and from point e to point f the energy are transferred from the inner Si layer to the outer water. From the field’s point of view, there should be a transition from point e to point f with the Si layer thickness increasing from 401.998 nm to 412.745 nm, but the mode is interspace cut-off.
So we make some other simulations to analyze the forming of the interspace cut-off region and we find that the interspace cut-off region will disappear if we increase the Ag layer thickness to some extent (e.g., 70 nm) with other parameters remained the same. The result is shown in
Fig. 5(b), which illustrates the real part of effective index with different Si layer thicknesses where the Ag layer thickness is 70 nm. It can be found from
Fig. 5(b) that the dispersion curves are continuous, which indicates that the possible cause of the existence of the interspace cut-off region in
Fig. 2 may lie in the thinness of the Ag layer. In
Fig. 2, the Ag layer is not thick enough to supply enough losses thus most of the mode energy flows propagate in water near interspace cut-off region (see S
z field of both point e and point f in
Fig. 5(c) and
5(d)), and the fields approach those of a plane wave in water together with the interspace cut-off of the mode. However, with a thicker Ag layer, less energy is supplied by the plane wave in water than is dissipated in the Ag layer, and the energy can be transferred from the inner Si layer to the outer water continuously with the Si layer thickness increasing. The disappearance of the interspace cut-off region with thicker Ag layer (e.g., 70 nm) strongly supports our explanations.
2.2 Sensor design
Here we are very interested in the mode properties around the interspace cut-off region, where the mode losses are very sensitive to the Si layer thickness variation. By taking TM3 mode for instance, as the Si layer thickness increases from 401 nm to 401.9 nm, the imaginary part of effective index falls from 0.0052 to 0.0004, which means that the mode can propagate 13 times longer in z direction. Since most of the mode energy propagates in water, the mode losses will also become sensitive to fluid refractive index change, which is the basic theory foundation of our SPP sensor.
Figure 7
shows the attenuation of TM
3 mode with four different fluid refractive index n
D around water. The attenuation is calculated by:
. We choose the Si layer thickness around the second interspace cut-off region of TM
3 mode (see region C in
Fig. 2 and
Fig. 3) for calculation. The four different detecting fluid refractive indices are 1.33, 1.3305, 1.333 and 1.335 with red, green, blue and purple lines illustrated in
Fig. 7, respectively. It can be seen clearly that with a properly chosen thickness of Si layer and device length, we can distinguish two fluids whose refractive index are very close to each other. This is caused by the mode characteristics that the attenuation curve is extremely sharp near the interspace cut-off region. For instance, with a Si layer thickness of 401.99 nm, the attenuation of n
D = 1.33 and n
D = 1.3305 are α
1 = 0.0005 dB/μm and α
2 = 0.0035 dB/μm respectively (see the black double arrow in
Fig. 7). As a result, with a device length of only 3.3 μm (0.01 dB/(α
2-α
1)), this sensor can distinguish the two fluids. Where we assume that the resolution of the power measurement is 0.01 dB.
Fig. 7 Relation between the attenuation of TM3 mode and Si layer thickness with four different fluid refractive index: nD = 1.33 (red line), nD = 1.3305 (green line), nD = 1.333 (blue line) and nD = 1.335 (purple line).
Moreover, our target lies in high resolution and large scale of refractive index detection, but there are always conflictions between resolution and detection scale. With a longer device, we can distinguish smaller difference of fluid index, but the larger attenuation induced by the longer device will decrease the scale of refractive detection, since the output power is too tiny to detect. So in consideration of resolution, compact in size and scale of detection, the device length is fixed at 200 μm. We choose TM
3 mode for sensing, as it exhibits a symmetric field distribution, which is easy for coupling in (see
Fig. 5(a), TM
3 mode exhibits a symmetric field distribution in the Si layer). For simplicity, we ignore the input mode coupling losses, which means the output power only depends on the losses of TM
3 mode in our SPP sensor. The output power is calculated as: P
out = P
inexp(−2αL), where P
in = 1 mW, L = 200 μm, α = Imag(n
eff) × 2π/λ.
Figure 8
shows the relation between output optical power of TM
3 mode and fluid refractive index n
D with three different Si layer thicknesses: d
Si = 402 nm (blue line), d
Si = 405 nm (red line) and d
Si = 410 nm (green line). Three linear dashed curves are given to compare with the curves of the detection so as to show the linear properties of the detection. According to the blue line and blue dashed line, the average sensitivity in the range of n
D = 1.33 to n
D = 1.41 can be as high as 1262 dB/RIU, which is equal to the slope of the blue dashed line. The sensitivity is calculated by: SA(RIU
−1) = Δα/Δn
DL [
15
J. Homola, S. S. Yee, and G. Gauglitz, “Surface plasmon resonance sensors: review,” Sens. Actuators B Chem.
54(1-2), 3–15 (1999). [CrossRef]
]. Since our SPP sensor is based on the intensity measurement, if the resolution of the power measurement is 0.01 dB, the minimum detectable refractive index change can be as small as 7.9 × 10
−6 RIU. In the case of n
D < 1.33, TM
3 mode is interspace cut-off, while in the case of n
D > 1.41, the linear properties of the detection gradually vanish. In spite of this, the result shows that in comparison with other sensors based on intensity measurement for refractive index detection [
7
A. Hassani and M. Skorobogatiy, “Design criteria for microstructured-optical-fiber-based surface-plasmon-resonance sensors,” J. Opt. Soc. Am. B
24(6), 1423–1429 (2007). [CrossRef]
,
9
R. Y. Wan, F. Liu, and Y. D. Huang, “Ultrathin layer sensing based on hybrid coupler with short-range surface plasmon polariton and dielectric waveguide,” Opt. Lett.
35(2), 244–246 (2010). [CrossRef]
[PubMed]
,
15
J. Homola, S. S. Yee, and G. Gauglitz, “Surface plasmon resonance sensors: review,” Sens. Actuators B Chem.
54(1-2), 3–15 (1999). [CrossRef]
], our designs exhibit larger scale of detection (e.g., 0.08, for 1550 nm ~1.33 to 1.41) and higher resolution (e.g., 7.9 × 10
−6 RIU) with a compact device size (length of 200 μm).
Fig. 8 Relation between output optical power of TM3 mode and fluid refractive index nD with three different Si layer thicknesses: dSi = 402 nm (blue line), dSi = 405 nm (red line) and dSi = 410 nm (green line). The device length is fixed at 200 μm. Three linear dashed curves are given to compare with the curves of the detection so as to show the linear properties of the detection which is, more specifically, scale of 1.33-1.41 for dSi = 402 nm, scale of 1.415~1.485 for dSi = 405 nm and scale of 1.552-1.592 for dSi = 410 nm.
Furthermore, the sensing range can be adjusted by changing the Si layer thickness, which can also be concluded from
Fig. 8. The red line corresponds to the Si layer thickness of 405 nm and as seen from it, good linear properties are obtained with a detecting range of n
D = 1.415 to n
D = 1.485 but the sensitivity falls to 1089 dB/RIU. So we find that although the mode losses are very sensitive to the Si layer thickness change (see
Fig. 7) which make the device difficult to be fabricated, small tolerance of the Si layer thickness will not affect the sensing range and the sensitivity too much (see
Fig. 8). The sensitivity can be compensated by increasing the device length, and the sensing range may be adjusted by changing the Ag layer thickness. This will be the topic of our next calculation.