University of Birmingham Terahertz particle-in-liquid sensing with spoof surface plasmon polariton waveguides

We present a highly sensitive microﬂuidic sensing technique for the terahertz (THz) region of the electromagnetic spectrum based on spoof surface plasmon polaritons (SPPs). By integrating a microﬂuidic channel in a spoof SPP waveguide, we take advantage of these highly conﬁned electromagnetic modes to create a platform for dielectric sensing of liquids. Our design consists of a domino waveguide, that is, a series of periodically arranged rectangular metal blocks on top of a metal surface that supports the propagation of spoof SPPs. Through numerical simulations, we demonstrate that the transmission of spoof SPPs along the waveguide is extremely sensitive to the refractive index of a liquid ﬂowing through a microﬂuidic channel crossing the waveguide to give an interaction volume on the nanoliter scale. Furthermore, by taking advantage of the insensitivity of the domino waveguide’s fundamental spoof SPP mode to the lateral width of the metal blocks, we design a tapered waveguide able to achieve further conﬁnement of the electromagnetic ﬁeld. Using this approach, we demonstrate the highly sensitive detection of individual subwavelength micro-particles ﬂowing in the liquid. These results are promising for the creation of spoof SPP based THz lab-on-a-chip microﬂuidic devices that are suit-able for the analysis of biological liquids such as proteins and circulating tumour cells in buffer solution. © 2017 Author(s).

We present a highly sensitive microfluidic sensing technique for the terahertz (THz) region of the electromagnetic spectrum based on spoof surface plasmon polaritons (SPPs). By integrating a microfluidic channel in a spoof SPP waveguide, we take advantage of these highly confined electromagnetic modes to create a platform for dielectric sensing of liquids. Our design consists of a domino waveguide, that is, a series of periodically arranged rectangular metal blocks on top of a metal surface that supports the propagation of spoof SPPs. Through numerical simulations, we demonstrate that the transmission of spoof SPPs along the waveguide is extremely sensitive to the refractive index of a liquid flowing through a microfluidic channel crossing the waveguide to give an interaction volume on the nanoliter scale. Furthermore, by taking advantage of the insensitivity of the domino waveguide's fundamental spoof SPP mode to the lateral width of the metal blocks, we design a tapered waveguide able to achieve further confinement of the electromagnetic field. Using this approach, we demonstrate the highly sensitive detection of individual subwavelength micro-particles flowing in the liquid. These results are promising for the creation of spoof SPP based THz lab-on-a-chip microfluidic devices that are suitable for the analysis of biological liquids such as proteins and circulating tumour cells in buffer solution. © 2017 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/). https://doi.org/10.1063/1.4998566

I. INTRODUCTION
The terahertz (THz) region of the electromagnetic spectrum, lying between microwaves and infrared light, is an important band for sensing due to its non-ionizing nature and the many intermolecular vibrational and rotational modes of polar molecules present. One area where THz sensing is predicted to have a significant impact is in biological sensing, 1,2 particularly, bio-liquids such as aqueous solutions of proteins, DNA, and cells. The use of THz waves in these applications is hindered by its relatively long wavelength compared with the size of the sample to be characterised. Therefore, techniques which enhance light-matter interactions in the THz regime through strong confinement of the electromagnetic field are highly desirable. 3 One approach to achieve this is to use metamaterials consisting of periodic metallic or dielectric resonators in which the localized electromagnetic fields are strongly enhanced. 4-8 A successful alternative has been to borrow the ideas of surface plasmon polaritons (SPPs) from the visible and near-infrared regions. The past few years have witnessed the rapid growth of research in plasmonics for its capability to provide strong field confinement and its high sensitivity to environmental a Electronic mail: s.maier@imperial.ac.uk 2378-0967/2017/2(11)/116102/9 2, 116102-1 © Author(s) 2017 refractive index changes. 9 However, in the terahertz frequency band, the high conductivity of metals only allows for poorly confined Sommerfeld-Zenneck waves to propagate on their surfaces. To overcome this, Pendry et al. proposed the use of structured metal surfaces with holes in them to mimic the SPPs that appear in the visible regime but at much lower frequencies, 10,11 as was experimentally verified by Hibbins et al. [12][13][14] These surface waves are termed spoof or "designer" surface plasmon polaritons. Since then, spoof SPPs have been extensively studied for the confined waveguiding of THz electromagnetic waves in various situations. [15][16][17][18][19][20][21][22][23] Spoof SPPs can be excited by matching their momentum to free-space radiation using the edge coupling method 13 or through a tapered mode converter to couple with conventional waveguides or transmission lines. 17,[24][25][26] In addition, the localized version of spoof surface plasmons can also be excited in closed geometries, and subwavelength electric and magnetic dipolar resonances have been demonstrated. 27,28 Recent researches have also shown that, in the mid-infrared and terahertz frequencies, surface plasmon can be excited and propagated on two-dimensional materials such as graphene and single-layer black phosphorus, in which the free carrier density can be tuned by chemical doping or electrical gating. [29][30][31] The most commonly employed spoof SPP platforms consist of corrugated metal surfaces or waveguide geometries, 11,32-34 including ultrathin waveguides. 16,17,21,35 In these designs, the optical properties of the spoof surface or waveguide can be readily controlled by the geometrical parameters. Spoof SPP waves can confine electromagnetic energy on subwavelength scales, greatly enhancing light-matter interaction with the local dielectric environment. In the THz regime, spoof SPPs provide the advantage of shrinking the relatively long wavelength of THz waves to a subwavelength scale comparable with the dimension of biological analyte targets. All these properties, together with the capability of guiding electromagnetic energy with spoof SPP waveguides, provide a new operating principle for optofluidic sensors, which combine the advantages of compact and rapid handling of microfluidics and the high detection sensitivity of the optical analysis. So far, a variety of waveguide-based optofluidic sensors have been demonstrated, including resonant microcavities, 36,37 photonic crystal resonators, 38 interferometer waveguide, 39 and others. 6,[40][41][42][43] Spoof surfaces as THz sensing platforms have been experimentally demonstrated for sensing both the real and imaginary parts of the refractive index of certain materials. [44][45][46] However, these sensing schemes operate with all of the grooves filled with the sample analyte, which may not be feasible for many situations where only a small volume of sample liquid is available. Moreover, the liquids are also not well guided, making integration with microfluidics for lab-on-chip sensing devices difficult.
In this paper, we demonstrate a microfluidic sensing platform consisting of a spoof SPP domino waveguide 34,47,48 and a micro-capillary through which a liquid sample can be introduced into the spoof SPP's field. We show through numerical simulations that high sensitivity detection of the refractive index of liquids can be achieved in the THz range. When the refractive index of the liquid inside the capillary is changed, a significant shift in the transmission spectrum of the spoof SPP waves occurs. The small interaction area between the micro-capillary and the spoof SPP mode, which is, in particular, determined by the lateral dimension of the waveguide, implies that a liquid volume only in the nanoliter scale is required. Furthermore, by using a tapered design with a gradually narrowing lateral block width, spoof SPP waves can be further confined in the lateral direction, making it possible to detect single micro-particles with diameters significantly smaller than the radiation wavelength passing through the capillary. The proposed sensing platform has the advantage that the liquid samples are completely contained in the microfluidic system, which can be useful for integrated lab-on-chip sensing devices in the future.

II. THEORY AND MODEL
The schematics of the two proposed spoof SPP waveguide optofluidic sensors are shown in has an outer diameter of 130 µm, a wall thickness of 10 µm, and a constant permittivity of 3.8. The homogeneous domino waveguide in Fig. 1(a) has periodic blocks of the same lateral width W. The tapered domino, sketched in Fig. 1(b), features a domino width that decreases slowly from the neighbour block, and the capillary is positioned at the center of the waveguide where W is narrowest. Experimentally, the domino blocks may be fabricated by deep reactive ion etching (DRIE) of silicon 4 or 3D printing of polymers 49 followed by metal deposition. In the THz regime, the skin depth of gold is very small compared with the wavelength of radiation, and the ohmic loss is very small. For simplicity, the metal is assumed to be a perfect electric conductor in all the simulations. We start by considering the dispersion relation of the spoof SPP modes that propagate on a domino waveguide. The dispersion relations are calculated using the finite element method (FEM) in COMSOL Multiphysics. In Fig. 2(a), the calculated dispersion curves are plotted for different values of W and the two-dimensional case which corresponds to W → ∞, with other  Fig. 1(a)]. k denotes the wavenumber of spoof SPP waves propagating in the x-direction, and it is normalized to the maximum value π/p. The modes that are below the light-line are confined non-radiative modes. For larger k values, the deviation from the light-line is greater and the mode is more confined to the surface, propagating at a slower group velocity. 19 The dispersion relation for an infinitely wide domino array, i.e., a spoof surface, can be analytically approximated by 15,34,50 where k 0 is the wavenumber of free space radiation. From the equation, one can see that the asymptotic frequency of the spoof SPP mode is mainly determined by the groove depth h. Deep grooves allow for the emergence of higher order modes in the vertical direction inside the grooves, 50 and in this paper, we concentrate only on the first order mode. In the three-dimensional case with finite lateral width W, it was previously demonstrated that the dispersion relation of spoof SPPs on domino blocks is relatively insensitive to the lateral width W. 34 In Fig. 2(a), we see that the asymptotic frequencies only change slightly with significantly different lateral widths. For W = 400 µm, the asymptotic frequency equals to 0.23 THz, and it is indicated with a dashed line parallel to the x axis. With smaller lateral width, the asymptotic frequency increases, and the overall spoof SPP mode appears to be less confined from the dispersion curve; however, the electric fields are more concentrated in the lateral y-direction due to the reduction in that dimension. In Fig. 2(b), the normalized electric (E)-field intensity and the E-field vectors at the asymptotic frequency are plotted in the transverse cross section (xz-plane). From the electric near-field distribution, we can see that the field intensity is most concentrated at the top of the domino, and the field decays exponentially away from that interface, similar to a conventional surface plasmon mode. In order to maximize the interaction between the liquid and the spoof SPP wave, the capillary is aligned with the top of the domino where the field enhancement is the greatest.
In Figs. 2(c) and 2(d), the normalized E-field distributions are plotted on the yz-plane at the middle of the grooves for widths W = 400 µm and W = 80 µm, showing more concentrated electric field for narrower domino width. Figures 2(c) and 2(e) both show that, in the lateral y-direction, the field intensity is relatively constant except for the hotspots at the block corners. It should be noted that in the domino structures with large lateral width, higher order transverse modes with varying charge distribution in the lateral direction are supported. These modes can have anti-nodes with zero field intensity in the center, 47 which is not beneficial for our purpose and we do not discuss them further in this paper.
For the simulation of the transmission behavior of the domino waveguide, an electric dipole oscillating in the x-direction is positioned at the beginning of the dominoes to excite the spoof SPP mode. Two-plane field monitors integrating the Poynting vector of the spoof SPP waves are set apart with a distance of 40 periods of dominoes. The transmission rate is calculated by normalizing the power integrated on the farther plane from the dipole to that from the nearer plane. In Fig. 2(e), the normalized E-field intensity on the homogeneous domino waveguide (W = 400 µm) without the capillary is plotted. This figure shows that the E-field amplitude remains constant over the dominoes, suggesting that there is no radiation or reflection loss along propagation. In contrast, the E-field amplitude for the tapered waveguide is plotted in Fig. 2(f), showing field enhancement in the center. In the study, the lateral width of the tapered domino waveguide starts at W = 400 µm and decreases in steps of 40 µm to the narrowest width of W = 80 µm, with an array of 4 domino blocks with constant W = 80 µm at the center. From Fig. 2(a), we observe that the domino with smaller lateral width has slightly higher cutoff frequency than the large width domino. Therefore, a tapered domino waveguide with slowly decreasing lateral width can support gradually confined spoof SPP waves without significant propagation loss. We can observe in Fig. 2(f) that in this tapered waveguide, the fields are more confined in the center region with smaller cross section, and it is expected that light-matter interaction between the SPP wave and the analyte will be enhanced in this region.

III. RESULTS AND DISCUSSION
In this section, we demonstrate how a change in the refractive index of liquid flowing through the micro-capillary or a single micro-particle passing through it affects the transmission of the spoof SPP waves on the domino waveguide, which can be used for sensing applications. The basic structures of the waveguide sensor are as shown in Fig. 1. First, we consider sensing of homogeneous liquids flowing through the micro-capillary channel [see Fig. 1(a)]. Liquids with different refractive indices are pumped through the capillary, and the transmission spectrum is measured for each case.
In the simulation, the tested materials correspond to nitrogen (n = 1), gasoline (n = 1.4), liquid paraffin (n = 1.5), glycerin (n = 1.82), and glucose (n = 2.1), 6,44 neglecting losses. In Fig. 3(a), the transmission (T ) spectra for a domino waveguide (W = 400 µm) with different liquids in the capillary are plotted. It is shown that when no capillary is inserted, the spoof SPP wave can propagate on the waveguide with nearly no radiative loss over the propagation band, maintaining a high transmission rate over 90%, until it meets the cutoff frequency of 0.23 THz where no spoof SPP mode can propagate. In the presence of a liquid-filled capillary, the transmissions start to drop at frequencies lower than the cut-off frequency of the domino waveguide, and the higher the index of the liquid, the more redshifted the transmission cutoff frequency is. This way, the transmission spectrum reflects changes in the refractive index value. These results demonstrate liquid sensing at THz frequencies with spoof SPPs using only a single groove, with greatly reduced interaction volume on the nanoliter scale. Such a small interaction volume also provides the advantage of preventing the transmitted signal from being completely attenuated by the high absorption of some lossy liquids in THz such as water.
To explore why the domino waveguide sensor has such high sensitivity to the index change in one of the grooves, we performed an eigenfrequency study for the unit domino structure with the capillary filled with materials of different refractive index. The dispersion curves of the spoof SPP mode are  plotted in Fig. 3(b). The dashed lines connecting Figs. 3(a) and 3(b) serve as a guide for the eye to show the connection between the cut-off frequencies of the spoof SPP modes and the corresponding reduction in the transmission of the domino waveguide with filled capillary. As shown in the figure, the higher the refractive index inside the capillary, the lower the cut-off frequency of the spoof SPP mode is because the presence of the capillary and liquids increases the average refractive index in the groove, effectively increasing the groove cavity length. From Figs. 3(a) and 3(b), one can see that the cutoff frequency of a groove with capillary is closely related to the frequency point where the transmission of the domino waveguide reduces to 1/e. As the liquid refractive index increases, the cut-off frequency of the groove decreases and the transmission spectrum displays a low-pass behavior with a lower cut-off frequency. In the lower frequency region of the transmission spectrum, the mode is less confined, with the spoof SPP wavelength being much larger than the capillary size. Therefore, the spoof SPP waves can pass around the capillary without too much scattering or reflection and the transmission rate is kept close to one, independent of the liquid refractive index. However, it should be pointed out that because the properties of only one groove are changed by the presence of the capillary, its effect is similar to scattering from a defect in a periodic structure. A complete transmission-forbidden bandgap is not formed, and spoof SPP waves can tunnel through the defect barrier and propagate along the waveguide. This explains why the transmission drop is not sharp as in the case of a domino waveguide without the capillary [see dark blue line in Fig. 3(a)]. To observe the effect of the capillary on the spoof SPP waves, we plot in Fig. 3(c) the near-field of the instantaneous x-component of the electric field on the homogeneous waveguide with capillary, with k = 0.5, and f = 0.21 THz, in the xz viewing plane. It can be seen that part of the spoof SPP field is scattered by the capillary and radiated into free space, and the transmitted spoof SPP wave has a greatly reduced amplitude.
The sensitivity of the system can be characterized by S = ∂T ∂n at a given frequency and liquid refractive index n, where T is the power transmission rate. In the homogeneous waveguide sensor, the maximum S has a value of 1.89/RIU at 0.19 THz. In Fig. 3(d), we show the transmission spectra for the tapered domino waveguide with the same liquids as in homogeneous waveguide. Overall, it shows the same trend as in the homogeneous domino waveguide, the transmission cut-off redshifts with increasing liquid refractive index. However, the shift is less significant compared to that in the homogeneous waveguide, and the sensitivity is also weaker with a calculated value of 1.54/RIU. The reduced performance in the tapered waveguide sensor for the liquid index change can be explained by the smaller interaction volume with the liquids due to the reduction in the lateral y-direction, and the fact that the slight rise of the cutoff frequency with smaller W also results in lower confinement of the spoof SPP waves in the z-direction. FIG. 4. Sensing a single micro-particle passing through the capillary. (a) Transmission spectra of a homogeneous domino waveguide with a single spherical particle in a liquid inside the capillary. The particle, with diameter 90 µm, is positioned at the center of the groove, and its refractive index takes values of n = 1, 1.73, 2.24. The background liquid in the capillary n = 1.73. (b) Transmission spectra for the tapered domino waveguide with capillary and particle, with particle parameters being the same as in (a). Now we discuss the performance of these two domino waveguide sensors for the detection of a sub-wavelength micron-sized particle passing through the capillary. In the simulations, a particle with diameter d = 90 µm is positioned at the center of the capillary, and the background liquid is set to have refractive index of 1.73. The transmission spectra for particles with different index values of 1, 1.73, and 2.24 are calculated and plotted in Fig. 4 for the homogeneous and the tapered domino waveguide sensors. By comparing the two plots, it can be seen that the tapered design has higher sensitivity than the homogeneous waveguide for detection of single micron sized particles, with sensitivity being improved from 0.15/RIU in the homogeneous waveguide to 0.35/RIU in the tapered design. In the homogeneous domino waveguide, the ratio of the particle volume to the whole spoof SPP-liquid interaction volume is rather small, so the change in the index of the particle does not result in a significant shift in the spoof SPP transmission spectra. In comparison, at the center of the tapered waveguide, where the cross section is narrowest (80 µm), the mode is squeezed in the lateral direction into a scale comparable with the particle's diameter, which makes the transmission more sensitive to the change of the particle index. Remarkably, at 0.2 THz, the free space wavelength is λ 0 ∼ 1500 µm, which is substantially larger than the particle's diameter, of only 90 µm. However, as we have shown, the enhanced field and wavelength shrinking provided by the spoof SPP enables sub-wavelength detection of the micro-particle with high sensitivity.
Finally, we note that in the discussion above, the liquids and domino waveguide materials have been assumed to be lossless. In practice, the materials will exhibit some loss and this will lead to the attenuation of the transmission signal, particularly near the cut-off frequency due to the tight confinement of the spoof SPP mode. For example, in our target frequency, based on the double Debye model, water has a high imaginary part of refractive index around 1-2, which will significantly attenuate the signal. A careful balance between the sensitivity and the tolerable loss needs to be optimized in this case. We suggest that one possible way to control the level of interaction between the spoof SPP waves and lossy analytes is by changing the vertical position of the capillary or operating at a frequency with a lower level of confinement. However, we stress that while materials are assumed to be lossless for simplicity, the operating principle and mechanism can be applied to realistic situations.

IV. CONCLUSION
In summary, we have numerically demonstrated a domino waveguide microfluidic sensor in the THz regime, where the spoof SPP modes can tightly confine THz waves. This greatly enhances light-matter interaction with the analytes in the capillary in a single groove. The transmission spectrum of the spoof SPP wave is highly sensitive to the refractive index of the analyte, which is on the scale of nanoliters, providing a new method for waveguide-based optofluidic sensing. Furthermore, by exploiting the insensitivity of the spoof SPP on the width of the domino blocks, the tapered domino waveguide with slowly shrinking width can greatly squeeze the spoof SPP mode into a narrower region with better field concentration, making it more promising for the detection of a sub-wavelength micro-particle. We have shown that the high sensitivity is due to the high refractive index liquid or particle decreasing the cut-off frequency of the spoof SPP mode in one groove, resulting in a significant reduction in transmission of the spoof SPP waves through this waveguide section. These findings can be of value for integrated microfluidic THz lab-on-chip devices for biosensing applications such as the analysis of proteins and detection of circulating tumor cells.