A microring resonator based negative permeability metamaterial sensor.

Metamaterials are artificial multifunctional materials that acquire their material properties from their structure, rather than inheriting them directly from the materials they are composed of, and they may provide novel tools to significantly enhance the sensitivity and resolution of sensors. In this paper, we derive the dispersion relation of a cylindrical dielectric waveguide loaded on a negative permeability metamaterial (NPM) layer, and compute the resonant frequencies and electric field distribution of the corresponding Whispering-Gallery-Modes (WGMs). The theoretical resonant frequency and electric field distribution results are in good agreement with the full wave simulation results. We show that the NPM sensor based on a microring resonator possesses higher sensitivity than the traditional microring sensor since with the evanescent wave amplification and the increase of NPM layer thickness, the sensitivity will be greatly increased. This may open a door for designing sensors with specified sensitivity.

Introduction

Due to their intriguing electromagnetic properties, a great deal of attention has been focused recently on metamaterials. The permittivity and permeability of metamaterials can be designed to continuously change from negative to positive values. Many novel metamaterial-based applications have been proposed, such as perfect lenses, cloaks, concentrators, directive antennas, superscatterers, superabsorbers, transparent devices, etc. [1-6]. Recently, great interest has been devoted to the sensing applications of metamaterials. For example, Jakšić et al. [7] investigated some peculiarities of electromagnetic metamaterials convenient for plasmon-based chemical sensing with enhanced sensitivity, and they envisioned practical applications of metamaterial-based sensors in biosensing, chemical sensing, environmental sensing, homeland security, etc. He et al. [8], studied the resonant modes of a 2D subwavelength open resonator, and showed it was suitable for biosensing. Melik et al. [9] presented telemetric sensing of surface strains on different industrial materials using split-ring-resonator based metamaterials, and desirable properties were obtained. Lee et al. [10] demonstrated experimentally the effectiveness of a split-ring resonator (SRR) array as a biosensing device at microwave frequencies. Cubukcu et al. [11] reported a surface enhanced molecular detection technique with zeptomole sensitivity that relies on the resonant electromagnetic coupling between a split ring resonator and the infrared vibrational modes of molecules. Alù et al. [12] proposed a method of dielectric sensing using ε near-zero narrow waveguide channels. Shreiber et al. [13] developed a novel microwave nondestructive evaluation sensor using a metamaterial lens for detection of material defects small relative to a wavelength. Zheludev [14] reviewed the road ahead for metamaterials, and pointed out that sensor applications are another growth area in metamaterials research. Our team has studied the performance of metamaterial sensors, and shown that the sensitivity and resolution of sensors can be greatly enhanced by using metamaterials [15-17].

WGM is a morphology-dependent resonance, which occurs when light within a dielectric microsphere, microdisk, or microring has a higher refractive index than its surroundings. In a ring resonator, WGMs form due to the total internal reflection of the light along the curved boundary surface [18]. The WGM resonance phenomenon has attracted increasing attention due to its high potential for the realization of microcavity lasers [19], quantum computers [20], sensing applications [21-29], etc. Examples of the applications of WGM sensors include biosensing [24], nanoparticle detection [25], single-molecule detection [26], temperature measurement [27], ammonia detection [28], and TNT detection [29]. However, to the best of our knowledge, there are no reports about any NPM sensors based on microring resonators operating in WGM.

In this paper, we derive the dispersion relation of a cylindrical dielectric waveguide loaded on a NPM layer, and compute the resonant frequencies and electric field distributions of the corresponding WGMs. We perform a full wave simulation of the performance of the NPM sensor, and compared it with the theoretical results. We show that the NPM sensor possesses much higher sensitivity than a traditional microring sensor, and the mechanism behind these phenomena is verified by theoretical analysis and simulation.

Theoretical Analysis

Figure 1 shows the geometry of a cylindrical dielectric waveguide loaded with a layer of metamaterials. The inner side of the e cylindrical dielectric waveguide (ε3,μ3) is loaded on a metamaterial layer (ε2,μ2). The waveguide has a four-layer structure. The material parameters of regions 1, 2, 3, 4 are denoted as (ε1,μ1), (ε2,μ2), (ε3,μ3), (ε4,μ4), respectively. The axial fields in corresponding regions for TM mode [30] are: where A, B, C, D, and are chosen here to weight the field, but they are interdependent. The functions J, Y, and K are, respectively, the Bessel functions of the first kind, of the second kind, and the modified Bessel function of the second kind. The terms , , , . β is the propagation constant, and m is the angular order. For an infinite cylindrical dielectric waveguide with negligible absorption and no axial component of the propagation constant (β = 0), TM mode degenerates to WGM [31], and Equation (1) becomes: where , , , , , is the Hankel function of the first kind. The relation between and K is . For TM mode in an infinite cylindrical dielectric waveguide, transverse magnetic fields can be obtained as:

Figure 1.

(a) Model of the four-layer cylindrical dielectric waveguide; (b) cross section of the waveguide.

The tangential fields matching equations at the boundary surfaces r=r1, r=r2 and r=r3 are expressed as:

Satisfying these conditions gives: where:

The dispersion equation can be obtained by setting | M |=0. The resonant frequency for different modes can be calculated using the software Mathematica (Wolfram Research Inc., Champaign, IL, USA). Coefficients B, , C, and can be expressed in terms of the arbitrary coefficient A, and , , , , . Parameters , , , and may be found from Equation (4). Electric field distribution for different mode can be obtained by substituting these coefficients in to Equation (2):

Results and Discussion

Simulation models of the NPM sensor based on a microring resonator are shown in Figure 2. A layer of NPM with thickness t is located on the inner side of the microring. Permittivity and permeability of the NPM are ε2 = ε0, μ2 = –μ0. Width of the microring and the waveguide is w = 0.3 μm. The outer diameter of the microring is d = 5 μm. The distance from outer microring to the waveguide is g = 0.232 μm. The permittivity of the microring and the waveguide is ε3 = 10.24ε0. Figure 2(a) is the simulation model for homogeneous sensing. The dielectric core with permittivity ε1 = εε0 is colored in light blue. Figure 2(b) is the simulation model for surface sensing. The dielectric substance with thickness t and permittivity ε1 = εε0 is attached to the NPM layer.

Figure 2.

Simulation model of the NPM sensor: (a) homogeneous sensing; (b) surface sensing.

The frequency spectrum of the NPM sensor for homogeneous sensing is simulated by the finite element software COMSOL Multiphysics (COMSOL Inc., Burlington, MA, USA), as shown in Figure 3. In the simulation, the computational space is surrounded by a scattering boundary. The excitation is set at port A of the waveguide. The spectrum is obtained by frequency sweep. From left to right, the spectral lines represent modes 25, 26, 27, 28 and 29 of the NPM sensor. The inset shows the amplification in the 191.83–191.87 THz frequency range. Table 1 shows the comparison of the analytical and simulated resonant frequency for the microring sensor and the NPM sensor. Therefore, WGMs (m = 25, 26, 27, 28, 29) in the cross section of the waveguide correspond to the modes of the microring sensor and the NPM sensor. The analytical resonant frequency of the sensor can be obtained by setting | M | = 0 (details may be found in next Section). The maximum deviation between simulation results and analytical results is 0.011 THz. The analytical results are in good agreement with the simulation results.

Figure 3.

Frequency spectrum of the NPM sensor. Thickness of the NPM layer is t = 0.09 μm. Permittivity of the dielectric core is ε = 1.

Table 1.

Comparison of the analytical frequency and simulated frequency for the microring sensor and the NPM sensor.

Mode (m)	25	26	27	28	29
Theoretical results for t = 0 μm (THz)	186.145	192.199	198.251	204.300	210.347
Simulation results for t = 0 μm (THz)	186.156	192.208	198.257	204.304	210.351
Deviation (THz)	0.011	0.009	0.006	0.004	0.004
Theoretical results for t = 0.12 μm (THz)	180.484	186.179	191.844	197.476	203.072
Simulation results for t = 0.12 μm (THz)	180.493	186.186	191.850	197.481	203.076
Deviation (THz)	0.009	0.007	0.006	0.005	0.004

Supposed that material parameters of the waveguide in region 1, 2, 3, 4 are ε1 = εε0, μ1 = μ0, ε2 = ε0, μ2 = –μ0, ε3 = 10.24ε0, μ3 = μ0, ε4 = ε0, μ4 = μ0, respectively. The resonant frequency of WGM (m = 27) in the cross section of the waveguide can be calculated by setting | M |=0. The coefficients , , , , can be easily obtained according to Equations (6)–(10). The electric field distribution of the WGM can be calculated according to Equation (2), and are shown in Figure 4(a,c). To confirm the WGM in the cross section of the waveguide corresponds to the mode of the microring resonator, we simulate the electric field distribution of the microring resonator, as shown in Figure 4(b,d). From Figure 4, we can observe that the theoretical results are in good agreement with the simulation results. Interestingly, we find that the maximum electric field is located at the interface of the NPM layer and core medium. This implies that a microring resonator loaded on an NPM layer has higher sensitivity than a traditional microring resonator without loading on the NPM layer.

Figure 4.

Electric field distribution of the WGM operating at mode 27. (a) The cross section of the waveguide; (b) the microring resonator; (c) the cross section of the waveguide loaded on NPM layer; (d) the microring resonator loaded on NPM layer. Thickness of the NPM layer is t = 0.09 μm.

To confirm the above idea, we simulated the performance of the microring sensor and the NPM sensor for homogeneous sensing, as shown in Figure 5. Permittivity (ε) of the dielectric core varies from 1 to 1.1 with an interval of 0.02. From Figure 5(a,b), we can observe that the spectra red shift with the increase of ε. Sensitivity for the microring sensor and the NPM sensor is 5.9 nm/RIU and 64.2 nm/RIU, respectively. Here, sensitivity is defined as . Figure 6(a,b) show the analytical and simulating resonant frequency for the NPM sensor and microring sensor, respectively. Simulating frequencies are calculated from Figure 5, while the theoretical frequencies are obtained by setting | M |=0. From Figure 6, we find that the simulation results are in good agreement with the theoretical results. With an increase of 0.02 in core medium permittivity, average frequency shift for the NPM sensor is very large [Figure 6(a)], but the average frequency shift of the microring sensor is quite small [Figure 6(b)]. Therefore, the NPM sensor possesses much higher sensitivity than the traditional microring sensor.

Figure 5.

Resonant frequency spectrum of mode 27 with respect to the change of core medium permittivity ε. From left to right, the curves correspond to ε = 1, 1.02, 1.04, 1.06, 1.08 and 1.1, respectively. (a) The microring sensor; (b) the NPM sensor. Thickness of the NPM layer is t = 0.09 μm.

Figure 6.

Relation between ε and resonant frequency. (a) NPM sensor; (b) Microring sensor.

To reveal the mechanism behind these phenomena, we plotted the electric field distribution of the NPM sensor along the x axis from −3 μm to −1.5 μm for mode 27, as shown in Figure 7. Permittivity of the core medium is set to be ε = 1. It is seen that the electric field intensity increases with NPM layer thickness (t). The inset shows the electric field distribution of the NPM sensor. From Figure 7, we can clearly observe that the stronger electric field of evanescent wave penetrates into the detecting region when the thickness of NPM layer increases. Figure 8 shows the relation between core medium permittivity and wavelength shift for different NPM layer thickness. Permittivity of the core medium increases from 1 to 1.1 with an interval of 0.02. Resonant wavelength shift is calculated by Δλ = λ(ε,t) – λ(1,t). For the microring sensor (t = 0), the sensitivity is only 5.9 nm/RIU. For the NPM sensor, the sensitivity increases with NPM layer thickness. When the thickness of the NPM layer is 0.06 μm, 0.09 μm, 0.12 μm, and 0.15 μm, the corresponding sensitivity will be 28.4 nm/RIU, 64.2 nm/RIU, 136.8 nm/RIU, and 240.7 nm/RIU, respectively. Therefore, the essence for the enhancement of sensitivity is the evanescent wave amplified by the metamaterial. Interestingly, we find that the sensitivity of the NPM sensor can be up to 327.3 nm/RIU when NPM thickness is 0.174 μm. But when the thickness is larger than 0.174 μm, WGM with m = 27 will be transferred to the WGM with m = 26 or 28. Details are not shown here for brevity.

Figure 7.

Electric field distribution along x axis from −3 μm to −1.5 μm for the NPM sensor operating in mode 27. The inset shows the electric field distribution of the NPM sensor, of which the NPM layer thickness is t = 0.15 μm.

Figure 8.

Homogeneous sensing. Relation between ε and wavelength shift for a variation of NPM layer thickness.

Surface sensing performance of the NPM sensor can also be analyzed according to the above procedures, and it is not shown here for brevity. Figure 9 shows the simulation results for surface sensing. Similarly, the sensitivity increases with NPM layer thickness. When the thickness of the NPM layer is 0.06 μm, 0.09 μm, 0.12 μm, and 0.15 μm, the sensitivity of the NPM sensor will be 24.1 nm/RIU, 54.9 nm/RIU, 117.7 nm/RIU, 208.9 nm/RIU, respectively. Therefore, sensitivity of the NPM sensor can be greatly improved by increasing the thickness of the NPM layer attached to its inner side. This is a novel method for sensor design with specified sensitivity.

Figure 9.

Surface sensing. Relation between ε and wavelength shift for a variation of NPM layer thickness.

Conclusions

WGMs of a dielectric waveguide with a layer of negative permeability metamaterial are theoretically analyzed, and the dispersion relation is derived. Analytical results of the resonant frequency shift and electric field distribution of the sensor are in good agreement with the simulation results. We show that the NPM sensor possesses a higher sensitivity than the traditional microring sensor, due to the amplification of the evanescent wave. Moreover, the sensitivity will be further improved by increasing the thickness of the metamaterial layer, opening a door for the design of novel sensors with desired sensitivity.