• Photonics Research
  • Vol. 7, Issue 8, 939 (2019)
Cheng Li1, Lei Chen1、2, Euan McLeod1, and Judith Su1、3、*
Author Affiliations
  • 1College of Optical Sciences, University of Arizona, Tucson, Arizona 85721, USA
  • 2State Key Laboratory of Information Photonics and Optical Communications, Beijing University of Posts and Telecommunications, Beijing 100876, China
  • 3Department of Biomedical Engineering, University of Arizona, Tucson, Arizona 85721, USA
  • show less
    DOI: 10.1364/PRJ.7.000939 Cite this Article Set citation alerts
    Cheng Li, Lei Chen, Euan McLeod, Judith Su, "Dark mode plasmonic optical microcavity biochemical sensor," Photonics Res. 7, 939 (2019) Copy Citation Text show less

    Abstract

    Whispering gallery mode (WGM) microtoroid optical resonators have been effectively used to sense low concentrations of biomolecules down to the single molecule limit. Optical WGM biochemical sensors such as the microtoroid operate by tracking changes in resonant frequency as particles enter the evanescent near field of the resonator. Previously, gold nanoparticles have been coupled to WGM resonators to increase the magnitude of resonance shifts via plasmonic enhancement of the electric field. However, this approach results in increased scattering from the WGM, which degrades its quality (Q) factor, making it less sensitive to extremely small frequency shifts caused by small molecules or protein conformational changes. Here, we show using simulation that precisely positioned trimer gold nanostructures generate dark modes that suppress radiation loss and can achieve high (>106) Q with an electric-field intensity enhancement of 4300, which far exceeds that of a single rod (2500 times). Through an overall evaluation of a combined enhancement factor, which includes the Q factor of the system, the sensitivity of the trimer system was improved 105× versus 84× for a single rod. Further simulations demonstrate that unlike a single rod system, the trimer is robust to orientation changes and has increased capture area. We also conduct stability tests to show that small positioning errors do not greatly impact the result.

    1. INTRODUCTION

    Whispering gallery mode (WGM) microtoroids [Fig. 1(a)] are excellent sensors due to their ultra-high-quality (Q) factors and their ability to be integrated on chip [19]. These sensors when combined with frequency-locking, balanced detection, and data processing techniques are able to detect wavelength shifts (<0.005  fm) that correspond to the detection of a biomolecule a radius of 2 nm binding to the microtoroid [1,2,4]. Coupling plasmonic particles to WGM cavities is one way to improve the sensitivity of the systems [10]. These particles increase the frequency shift upon molecular binding by providing enhanced interaction between the optical field and the analyte due to localized surface plasmon resonance (LSPR) [10,11]. The scattering and absorption losses of the metal, however, cause degradation of the Q factor of the resonance, that is, the broadening of the linewidth. Dark mode plasmonic resonances can help solve this problem [1216], but simulation is necessary to understand the specific impact of a given plasmonic structure on the system.

    Schematic of a microtoroid cavity. (a) The E-field is normalized by the amplitude of the maximum field in the evanescent zone of the bare WGM toroid [Eo(r)]. (b) A rendering of a gold nanorod placed parallel to the TE polarization of the WGM cavity mode. The resonance frequency of a single rod is tuned by adjusting its aspect ratio, which is defined as the ratio of the length to the width (diameter) of the rod. (c) Field distribution of the excited dipole mode around a nanorod.

    Figure 1.Schematic of a microtoroid cavity. (a) The E-field is normalized by the amplitude of the maximum field in the evanescent zone of the bare WGM toroid [Eo(r)]. (b) A rendering of a gold nanorod placed parallel to the TE polarization of the WGM cavity mode. The resonance frequency of a single rod is tuned by adjusting its aspect ratio, which is defined as the ratio of the length to the width (diameter) of the rod. (c) Field distribution of the excited dipole mode around a nanorod.

    Currently, due to limitations in computing performance and time consumption, a wedge-shaped model with perfect electric conductor (PEC) boundary conditions is widely used in plasmonically enhanced WGM toroid simulations [17]. However, such boundary conditions act as mirrors, effectively replicating the plasmonic particle multiple times, which results in an inaccurate calculation of the coupled system Q factor for a single particle (see Appendix A for more details).

    Here, we use a three-dimensional (3D) eigenmode simulation model of a whole microtoroid with a major diameter of 10 μm to explore the interaction of a single cavity with one or several nanostructures. The model is implemented in COMSOL. Due to simulation constraints in both time and memory, we cannot accurately simulate larger whole toroids in 3D. The schematic of our 3D model and the field distribution of the fundamental TE mode are shown in Fig. 1(a). The polarization of the TE mode is perpendicular to the equatorial plane of the toroid cavity. Nanorods are placed perpendicularly to the equatorial plane for maximum coupling and excitation, and their near-field enhancement is shown in Fig. 1(c).

    We define a figure of merit known as a combined enhancement factor (fC) to quantify the limit of detection and resolution of our sensing system. With this figure of merit, we can predict that slightly off-resonant coupling between the plasmon and WGM provides better performance than a direct match of the WGM and plasmon resonances. Further, we discuss the improvements that a lateral dark mode supported by a plasmonic trimer structure brings to the system and verify its stability.

    Two critical parameters affect the limit of detection of a WGM sensing system. One is the Q factor mentioned above [18], which determines the linewidth of the resonance peak in the output spectrum. Higher Q factors correspond to narrower resonances, whose central resonance frequency can be more precisely determined and tracked [13,19]. The other parameter is the electric (E)-field enhancement factor. The magnitude of the shift in resonance frequency of the WGM upon interaction with an analyte is proportional to this enhancement factor. According to perturbation theory [20], the frequency shift caused by a particle of volume ΔV with (possibly anisotropic) material parameters Δε (r) and Δμ (r) relative to the background is given by [2022] where ε0 and μ0 are the permittivity and permeability of free space, ε and μ are relative permittivities and permeabilities, εb is the permittivity of the background media (in this case, vacuum), E0 and H0 are the E and magnetic fields for an unperturbed WGM, represents conjugate-transpose, E and H are the electric and magnetic fields within the perturbation volume, |E0|max2 is the maximum intensity within the unperturbed toroid, Vm=Vεr|E0|2dV/(εr|E0|max2) is the mode volume of the resonant mode [19], εr is the dielectric constant of silica [23], and α(ω) is the polarizability tensor of the particle. In the numerator, we have assumed Δμ=0, and in the denominator, we have assumed that the magnetic field contributes an equal amount of energy as the electric field and that the perturbation has a negligible impact on the total energy stored in the cavity. The polarizability tensor for a metallic nanorod can be expressed as [2426] where α and α correspond to the longitudinal and transverse modes of the nanorod, respectively. The nanorods are small enough that they can be accurately modeled based on the LSPR of prolate spheroids [25]: where ΔV=43π(d2)3+π(d2)2(Ld) is the volume of the cylindrical nanorod, εm is the relative permittivity of the metal, εb is the dielectric constant of the surrounding environment, and the G, are geometrical factors that can be, respectively, written as where e=1(dL)2 is the eccentricity of the prolate spheroid, and Rs is a parameter that describes the particle shape. Here, we use Rs=0.88 for a spherocylinder [27]. In the simulation, we put the rod parallel to the polarization of the fundamental TE mode and study only the influence of the longitudinal plasmon mode and its corresponding α. It is apparent that α is greatly influenced by the nanorod aspect ratio. This in turn affects the shift of the unperturbed WGM.

    2. COMBINED ENHANCEMENT FACTOR

    Figure 2(a) illustrates the frequency shift of the cavity mode around on-resonance (λ774  nm) coupling obtained in the simulation, which is consistent with both experimental data and perturbation theory [28]. The data for the perturbation theory is obtained by numerical calculation of Eqs. (1)–(4), where only the longitudinal polarizability α is necessary, and the unperturbed E fields and mode volumes are determined from a simulation of the bare cavity. In our calculations for different aspect ratios, the nanorod diameter is held fixed, while the length is varied, so at large aspect ratios, the assumption is that the nanorod size λ starts to fail, and the theory deviates from the numerical simulation results [29]. At the same time, the linewidth of the coupled system also varies with the resonance frequency of a single rod, as shown in Fig. 2(b).

    (a) Blue shift or red shift of the cavity mode around on-resonance coupling. (b) The relationship between the linewidth corresponding to system loss and the resonance of a single rod. (c) Q factor and enhancement factor as functions of the resonance for a single nanorod. The trend for the enhancement factor is similar to the linewidth change in (b). Extremely strong enhancements are shown for on-resonance coupling. Due to the light–matter interaction, a very strong hot spot is generated between the plasmonic nanorod and the biomolecule. (d) The relationship between the combined enhancement factor (fC) and the resonance frequency of the rod.

    Figure 2.(a) Blue shift or red shift of the cavity mode around on-resonance coupling. (b) The relationship between the linewidth corresponding to system loss and the resonance of a single rod. (c) Q factor and enhancement factor as functions of the resonance for a single nanorod. The trend for the enhancement factor is similar to the linewidth change in (b). Extremely strong enhancements are shown for on-resonance coupling. Due to the light–matter interaction, a very strong hot spot is generated between the plasmonic nanorod and the biomolecule. (d) The relationship between the combined enhancement factor (fC) and the resonance frequency of the rod.

    The inset of Fig. 2(c) shows how we calculate the intensity enhancement at the particle, which is done by adding a sphere in the near field of the nanorod to represent an analyte molecule. Here, we use the volume integration of the E-field intensity over the volume of an analyte particle as the enhancement factor instead of the maximum E-field intensity at the hot spot of the nanoparticle (NP) [10,30]. This is more accurate because the maximum E-field hot spot intensity is sensitive to the mesh of the simulation model and does not take into account the frequency shift caused by the size of the measured particle. The ratio of the integrated field involving a coupled plasmonic particle to that of the evanescent field of the bare toroid can be used to define a field enhancement factor as follows: where |Ec(r)| is the E-field magnitude of the near field of the coupled system. The enhancement factor corresponding to the E-field intensity directly affects the magnitude of the frequency shift when a particle binds to the cavity. This enhancement factor closely follows the same trend as the linewidth broadening caused by the nanorod [Fig. 2(b)]. It is evident that introducing a particle with a lossy material such as gold will significantly decrease the Q factor of the system, even with strong local E-field enhancement, so there is a trade-off between Q and the intensity enhancement factor, as shown in Fig. 2(c). At a wavelength of λ774  nm, the maximum 2500× intensity enhancement is achieved when the rod aspect ratio equals 5.4. In this case, the peak LSPR wavelength coincides with the particular WGM resonance being simulated. The unperturbed, bare toroid exhibits a Q factor of 7×106, which drops to 1.87×105 when coupled to the nanorod with an aspect ratio 5.4.

    The lower Q of the coupled system involving resonant nanorods may make it difficult to precisely lock to resonance peaks in experiments involving tracking the wavelength shift of the mode. Therefore, a balance is desired between the Q factor and the intensity enhancement factor.

    From Eqs. (1) and (3), it can be shown that the frequency shift is directly proportional to the E-field intensity times the volume of the detected particle: ΔωE2ΔV, for scalar α. Then, the volume of the smallest detectable particle follows, ΔVminΔωmin/|E|2(Q|E|2)1, assuming that the smallest experimentally measurable shift, Δωmin, is inversely proportional to the Q factor because high-Q resonances correspond to narrow peaks that can be tracked with greater precision than broad resonances. So, the volume of the smallest detectable particle is inversely proportional to a figure of merit that accounts for both the E-field and Q-factor effects, which we call the combined enhancement factor (fC): where Qc and Q0 represent the Q factors of the coupled and bare toroid cavities, respectively. Higher fC values correspond to better sensitivity and/or more well-defined resonances. Unlike the peak in fE, the highest fC does not occur when the LSPR and WGM resonances coincide. Instead, as shown in Fig. 2(d), the position of the peak is shifted significantly to a smaller aspect ratio than that of resonant coupling. We can state that the improvement to the system from a single rod with an aspect ratio of 5.0 is ideal, which corresponds to a 874 field enhancement factor and a Q factor of 6.8×105. Using fC as a figure of merit is helpful in selecting the best plasmonic structure for experimental use.

    3. INDIVIDUAL NANOROD COUPLED SYSTEM

    For multiple gold nanoantennas, we investigate the intrinsic properties of the coupled system. The single gold nanorod can be approximated as a nanocavity. According to coupled mode theory [21,3133], the equations of several coupled modes can be written as where bWGM and bn represent the mode amplitude of the WGM mode and LSPR mode of the nth nanorod, ΩWGM and Ωn correspond to their resonance frequencies, and ΓWGM and Γn represent their corresponding losses, respectively, which can be calculated from COMSOL using the eigenfrequency solver. The coupling coefficient is defined as a complex number κ, which takes into account the frequency shift of the mode and the degradation of the Q. Because the rods in the array are spaced relatively far apart, we only consider the coupling of each LSPR mode with the cavity and no direct interparticle coupling. Therefore, the above equations for the coupled system can be written as a matrix equation M ˙=HM, where M=[bWGM,b1,,bn]T, and The intrinsic Q factor of the coupled system can be calculated from the eigenvalues of H, corresponding to the steady state of the system. Through the simulation data, the coupling coefficient κ (1.43×1011  Hz) can be obtained from a 2×2 matrix generated by a system in which a cavity is coupled to a single rod (see Appendix B). Further, the Q factor of a (n+1)×(n+1) matrix corresponding to a multiple (n) rod coupled system can be calculated numerically. We compare these values with the Q factors obtained from COMSOL simulation of the multi-rod system [Fig. 3(a)]. The coupled mode theory can then be used to predict the interaction between the WGM and additional particles without having to run new COMSOL simulations.

    (a) Q factors obtained through both numerical simulations and coupled mode theory are consistent for the systems involving multiple nanorods and no direct inter-rod coupling. (b) Top view of multiple rods coupled to the cavity mode. (b) and Fig. 2(c) share the same color bar.

    Figure 3.(a) Q factors obtained through both numerical simulations and coupled mode theory are consistent for the systems involving multiple nanorods and no direct inter-rod coupling. (b) Top view of multiple rods coupled to the cavity mode. (b) and Fig. 2(c) share the same color bar.

    Figure 3(a) also shows the full simulation results of coupling several rods to the cavity, where each rod is placed at an antinode of the standing wave WGM [Fig. 3(b)]. Here, the nanorods are sufficiently far apart such that direct coupling between neighboring nanorods is negligible compared to nanorod-WGM coupling, as verified in Fig. 3(a) by the agreement between the numerical simulation results and coupled mode theory [eigenvalues of H, Eq. (9)].

    4. TRIMER-INDUCED LATERAL DARK MODE

    The dipole modes induced in isolated antennas lead to large energy radiation, which adversely limit the Q factors of the coupled systems. Dark modes, an example of which is the antibonding interaction of dipolar resonances, have a longer lifetime and lower loss compared to bright dipolar modes, but their excitation is not easily feasible from the far field [12,13,16]. Due to the nearly zero net dipole moment of dark plasmons, it is almost impossible for their out-of-phase state to interact with light in free space.

    Here, we designed a plasmonic trimer consisting of three isolated rods arranged in a symmetric triangle [34,35] [Fig. 4(a)] such that the microtoroid WGM excites one of the nanorods in the trimer, which in turn excites a dark mode resonance (lateral dark mode) of the other two nanorods [12,15,3638]. Coupling of the dark mode and WGM standing wave can significantly reduce degradation of the Q factor. The superposition of rod near-fields at the center of the trimer can also bring greater intensity enhancement factors. Figure 4(b) shows the spectrum of the lateral dark mode of a trimer and three individual rods in free space. When the ends of three rods exhibit the same charge distribution, we call it a breathing dark mode. Only when two antennas of the trimer generate an antibonding mode [Fig. 4(b)] do we refer to this as a lateral dark mode. Although the trimer spectrum exhibits a peak near 800 nm upon far-field or WGM excitation, we still call its mode a lateral dark mode because the laterally oriented lower two rods of the trimer are excited in an antibonding, or dark, mode. The apparent resonance peak in Fig. 4(b) is due to the excitation of the vertical rod, which remains bright. Because of this simultaneous combination of bright and dark modes, it is possible to excite a plasmon resonance with minimal scattering and absorption losses. Although the breathing dark mode has much lower loss, and its corresponding resonance peak is almost covered by the tail of the lateral dark mode, its coupling to the WGM is very weak. Therefore, this paper mainly studies the lateral dark mode instead of the breathing dark mode.

    (a) Plane wave excitation for three isolated rods and a gold trimer. (b) Spectral comparison of the total extinction cross section of the lateral dark mode versus the three isolated rods from (a). The lateral dark mode is excited at the peak wavelength of the red curve. The illustrations show the charge distribution of the breathing and lateral dark modes at different wavelengths. Because the coupling between the breathing dark mode and free-space radiation is so small, no peak is visible at its resonance around 725 nm. (c) Quality and intensity enhancement factors as functions of the aspect ratios of each individual rod. The red arrows show the current density direction obtained in COMSOL. The inset shows the field distribution of the excited lateral dark mode. The characteristic dark spot between the ends of the bottom two rods is clearly visible. (d) Plot of the fC of the trimer. The lateral dark mode exhibits a higher combined enhancement value than that obtained from the coupling of a single rod to the cavity alone.

    Figure 4.(a) Plane wave excitation for three isolated rods and a gold trimer. (b) Spectral comparison of the total extinction cross section of the lateral dark mode versus the three isolated rods from (a). The lateral dark mode is excited at the peak wavelength of the red curve. The illustrations show the charge distribution of the breathing and lateral dark modes at different wavelengths. Because the coupling between the breathing dark mode and free-space radiation is so small, no peak is visible at its resonance around 725 nm. (c) Quality and intensity enhancement factors as functions of the aspect ratios of each individual rod. The red arrows show the current density direction obtained in COMSOL. The inset shows the field distribution of the excited lateral dark mode. The characteristic dark spot between the ends of the bottom two rods is clearly visible. (d) Plot of the fC of the trimer. The lateral dark mode exhibits a higher combined enhancement value than that obtained from the coupling of a single rod to the cavity alone.

    Due to the plasmonic coupling between closely placed nanorods, the resonance position of the lateral dark mode shifts with respect to the spacing and position of the rods. In Fig. 4(b), we compare the extinction cross sections when three rods with the same aspect ratio of 5.4 are placed separately or placed as a trimer with an 8 nm central gap. The extinction cross sections are normalized by that of the three isolated rods. The maximum extinction of the trimer is about half of that obtained from three individual rods. At many wavelengths, including 774 nm, corresponding to the WGM resonance, the loss from the trimer is even lower than that of a single rod. Such a low-dissipation mode improves the WGM sensing system.

    We also adjust the aspect ratio of the three rods in the trimer to optimize the improvement of the trimer coupling system. Figure 4(c) shows how the Q factor and intensity enhancement vary with the rod aspect ratio. The biggest intensity enhancement factor here corresponds to the case where the spherical analyte particle is placed in the center of the trimer. At an aspect ratio of five, the maximum field enhancement factor (4300) is much larger than that of a single rod (2500), as mentioned above. At the same time, the Q factor of 1.31×105 is similar to that of the single rod case (1.87×105) at on-resonance coupling and is twice as large as the Q for the three rod case (5.98×104) presented in Fig. 3(a). Even the relatively dark region just below the center of the trimer in the inset of Fig. 4(c) can bring dozens of times of intensity enhancement, which also demonstrates that the trimer structure is beneficial to increasing the capture area for biomolecule detection. The maximum fC of the system is 105, which is achieved using the trimer with a rod aspect ratio of 4.5. This fC is also significantly larger than that of a single rod (84).

    To test the robustness of our system, we studied the effect of tilted angle perturbation, length perturbation, and inter-rod spacing on fC. Here, we chose the trimer with a rod aspect ratio of 4.5, mentioned above as a reference. As shown in the insets in Fig. 5, we observe the difference of fC in the system by introducing different geometric perturbations. Considering the need to bind biomolecules at the center of the trimer in experiments, the spacing between the rods should not be too small. Although different perturbations will impact the Q factor or enhancement factor of the system, the overall fC changes little from its initial value of 105. The spacing in Fig. 5(c) has little effect on the loss caused by the structure, so the corresponding Q value is almost unchanged, and the observed variation is primarily due to E-field enhancement.

    Influence of different perturbations on the trimer’s combined enhancement factor. (a) and (b) Small changes in angle and length can maintain the fC. (c) To ensure that the gap space is large enough for particles to bind, we study the effect of the spacing on the system when the spacing is greater than 5 nm. The illustrations for the trimer field distribution use the same color bar as in Fig. 3. (d) The fC of a single rod system decreases with increasing rotation angle. When the rod is rotated 90 deg, that is, perpendicular to the polarization of the TE mode, the overall system improvement fC is only 4. The white area in the inset shows where the magnitude of the E-field is 10× greater than that of a bare cavity. (e) The rotation of the trimer has little effect on fC. The white area in the inset shows where the magnitude of the E-field is 10× greater than that of a bare cavity. The enhancement area provided by the trimer in (e) is slightly more than twice that of a single rod in (d).

    Figure 5.Influence of different perturbations on the trimer’s combined enhancement factor. (a) and (b) Small changes in angle and length can maintain the fC. (c) To ensure that the gap space is large enough for particles to bind, we study the effect of the spacing on the system when the spacing is greater than 5 nm. The illustrations for the trimer field distribution use the same color bar as in Fig. 3. (d) The fC of a single rod system decreases with increasing rotation angle. When the rod is rotated 90 deg, that is, perpendicular to the polarization of the TE mode, the overall system improvement fC is only 4. The white area in the inset shows where the magnitude of the E-field is 10× greater than that of a bare cavity. (e) The rotation of the trimer has little effect on fC. The white area in the inset shows where the magnitude of the E-field is 10× greater than that of a bare cavity. The enhancement area provided by the trimer in (e) is slightly more than twice that of a single rod in (d).

    Rotations of the entire structure are evaluated for single nanorods and trimer structures in Figs. 5(d) and 5(e). The gold nanorods in current plasmonic-enhanced WGM systems are randomly attached on the surface of the cavity, and their orientation is challenging to control. Here, we use a rod with an aspect ratio of five (Fig. 2) to study the influence of the rotation angle on fC. When the rotation angle is greater than 45 deg, the fC decreases as the rotation angle increases [Fig. 5(d)] due to mismatch with the polarization direction of the WGM. The system has only 4× gain when rotated 90 deg. In comparison, the trimer is robust to orientation changes. Due to its rotational symmetry, we only rotate the whole trimer by 0–60 deg [Fig. 5(e)]. The 0 deg point corresponds to the idealized trimer composed of three rods with an aspect ratio of 4.5. The weakest fC at a 30 deg angle is 92, which is only 11.5% lower than the maximum enhancement of 104. The trimer system is much less affected by the rotation angle than the single rod system. The inconsistent signal amplification caused by the rotation of a single rod can add uncertainty to particle detection experiments.

    5. CONCLUSIONS

    In this work, we present a lateral dark mode of a symmetric trimer structure generated using WGM excitation and define a combined enhancement factor to analyze the improvements it brings to the system. The trimer exhibits a 105 combined enhancement factor versus 84 for a single rod. We also demonstrate the stability and robustness of the trimer structure. We verify the correctness of our simulation from the frequency shift and linewidth changes caused by detuning and the effect on the Q factor. Compared to the existing single rod plasmonic hybrid WGM system, the trimer system exhibits higher enhancement, stability, and greater contact area for particle detection. In future work, we plan to use high-precision optical tweezers [39,40] to position trimer microresonator systems for biological sensing experiments.

    APPENDIX A: Details of the Full 3D Toroid Simulation and its Comparison with the Wedge Model

    The data and profiles are obtained in COMSOL using finite element analysis. The material of the cavity is silica with a refractive index of 1.45+10 8i, and the background is air [18] with a refractive index of one. We add the imaginary component of the refractive index of the cavity to lower the Q factor to a value similar to what we would see in experiment. The material parameters of gold used here are those given by Johnson and Christy [41]. The diameter of each gold rod is 6 nm. The simulated biomolecule used to calculate the enhancement factor is a sphere with a refractive index of 1.5 and a radius of 1 nm. We constructed the mesh of the cavity by using a uniform sweeping method with a maximum size of one-ninth of a wavelength. The Q factor and field distribution of the unperturbed cavity surrounded by a perfectly matched layer are consistent with the analytical solution and the two-dimensional axisymmetric model [4245]. The Q factor is calculated from the real and imaginary parts of the eigenfrequency obtained by COMSOL. In the matrix given in Eq. (5), the parameters ΩWGM (3.8724×1014 Hz) and ΓWGM (1.03×109 Hz) are determined from simulations of a bare toroid and the interaction of a single rod with a cavity using the eigenfrequency solver in COMSOL. Ωn (3.8724×1014 Hz) and Γn (1.479×1012 Hz) are obtained from the spectrum of the NP in free space, and all rods are identical. The server used for the simulation is configured as a 56-core 2.3 GHz processor with 384 GB memory.

    The wedge simulation model is currently widely used in the calculation and prediction of relevant experimental data for microresonators [17]. The model uses PEC boundary conditions that can be considered as mirrors. Such boundary conditions apply only to bare cavities, but when the cavity is coupled to an NP or structure, the PEC boundary conditions replicate the NPs multiple times, such that the interaction between one NP and the WGM cannot be accurately predicted. This part compares the results from a wedge model with those of the full 3D model mentioned in the main body of the paper.

    Here, we study the effect of the wedge angle on the Q value. The Q value decreases with the increasing number of coupled NPs (Table 1). To compare the wedge model with a whole 3D model, we construct both models for a toroid with a major radius of 5 μm. Figure 6(a) shows the schematic of the wedge model. The wedge angle θ is inversely proportional to the number of azimuthal modes m and can be written as θ=Nπm,where N is an integer corresponding to the number of antinodes of the standing wave in Fig. 6(b). m satisfies the resonance condition mλ=2πneffR, where neff represents the effective refractive index of the mode. The N in the wedge structure is taken as 1, 3, 5, 7 and compared with the whole toroid model that is used in the main text. The corresponding field distribution for the N=7 case is shown in Figs. 6(b) and 6(c).

    (a) Schematic of the wedge model. (b) E-field distribution of a particle coupled to a toroid cavity (top view). The particle is placed at the antinode. (c) A zoom-in of near-field enhancement. (d) and (e) show the SM and ASM of the coupled system obtained by the whole 3D model using the eigenfrequency solver, respectively. The corresponding detuning frequency between the SM and ASM modes is 200 MHz.

    Figure 6.(a) Schematic of the wedge model. (b) E-field distribution of a particle coupled to a toroid cavity (top view). The particle is placed at the antinode. (c) A zoom-in of near-field enhancement. (d) and (e) show the SM and ASM of the coupled system obtained by the whole 3D model using the eigenfrequency solver, respectively. The corresponding detuning frequency between the SM and ASM modes is 200 MHz.

    The Q of the wedge structure is significantly lower than that of the whole toroid and one particle. This may be because the larger the wedge angle, the lower the number of coupled particles, so the higher the Q. But we do not know the interaction of a single rod and a cavity, nor can we calculate the influence of boundary conditions on the E field and its influence on the coupled system.

    Furthermore, when a nanoscatterer is coupled to the WGM cavity, the coupling of the scattered light to the clockwise mode and counter-clockwise mode of the cavity will generate both a symmetric mode (SM) and an asymmetric (ASM) mode [46]. SMs and ASMs correspond to the antinode or node of the standing wave where the particle is located, respectively. Since the PEC boundary condition sets the E field to zero at the boundary, the position of the WGM standing wave at a specific wavelength is fixed. Accordingly, the introduced NP does not interact with both the SM and ASM simultaneously. The 3D whole toroid model, however, can generate two modes, as shown in Figs. 6(d) and 6(e). This is another advantage of the whole toroid model. Therefore, the whole toroid model is used here rather than the wedge model.

    APPENDIX B: Coupling Coefficient Calculation

    ΩWGM+iΓWGM and Ω1+iΓ1 can be obtained from the eigenfrequency in the bare toroid simulation and the isolated nanorod, respectively. If we assume that the complex eigenfrequency of the coupled system is Ωcoupled+iΓcoupled, which can be solved from the simulation of the bare toroid coupled to a single nanorod, then from the analytical eigenvalue of the 2×2 matrix, we can calculate κ using κ=[(ΩWGM+iΓWGM Ωcoupled iΓcoupled) (Ω1+iΓ1 Ωcoupled iΓcoupled)]12.

    References

    [1] J. Su. Label-free single exosome detection using frequency-locked microtoroid optical resonators. ACS Photon., 2, 1241-1245(2015).

    [2] J. Su. Label-free biological and chemical sensing using whispering gallery mode optical resonators: past, present, and future. Sensors, 17, 540(2017).

    [3] J. Su, A. F. Goldberg, B. M. Stoltz. Label-free detection of single nanoparticles and biological molecules using microtoroid optical resonators. Light Sci. Appl., 5, e16001(2016).

    [4] J. D. Swaim, J. Knittel, W. P. Bowen. Detection limits in whispering gallery biosensors with plasmonic enhancement. Appl. Phys. Lett., 99, 243109(2011).

    [5] F. Vollmer, L. Yang. Review label-free detection with high-Q microcavities: a review of biosensing mechanisms for integrated devices. Nanophotonics, 1, 267-291(2012).

    [6] Y. Zhi, X.-C. Yu, Q. Gong, L. Yang, Y.-F. Xiao. Single nanoparticle detection using optical microcavities. Adv. Mater., 29, 1604920(2017).

    [7] D. K. Armani, T. J. Kippenberg, S. M. Spillane, K. J. Vahala. Ultra-high-Q toroid microcavity on a chip. Nature, 421, 925-928(2003).

    [8] M. R. Foreman, J. D. Swaim, F. Vollmer. Whispering gallery mode sensors. Adv. Opt. Photon., 7, 168-240(2015).

    [9] F. Vollmer, S. Arnold, D. Braun, I. Teraoka, A. Libchaber. Multiplexed DNA quantification by spectroscopic shift of two microsphere cavities. Biophys. J., 85, 1974-1979(2003).

    [10] M. D. Baaske, M. R. Foreman, F. Vollmer. Single-molecule nucleic acid interactions monitored on a label-free microcavity biosensor platform. Nat. Nanotechnol., 9, 933-939(2014).

    [11] S. I. Shopova, R. Rajmangal, S. Holler, S. Arnold. Plasmonic enhancement of a whispering-gallery-mode biosensor for single nanoparticle detection. Appl. Phys. Lett., 98, 243104(2011).

    [12] Y. Gao, N. Zhou, Z. Shi, X. Guo, L. Tong. Dark dimer mode excitation and strong coupling with a nanorod dipole. Photon. Res., 6, 887-892(2018).

    [13] T.-S. Deng, J. Parker, Y. Yifat, N. Shepherd, N. F. Scherer. Dark plasmon modes in symmetric gold nanoparticle dimers illuminated by focused cylindrical vector beams. J. Phys. Chem. C, 122, 27662-27672(2018).

    [14] D. E. Gómez, Z. Q. Teo, M. Altissimo, T. J. Davis, S. Earl, A. Roberts. The dark side of plasmonics. Nano Lett., 13, 3722-3728(2013).

    [15] D. E. Gómez, K. C. Vernon, T. J. Davis. Symmetry effects on the optical coupling between plasmonic nanoparticles with applications to hierarchical structures. Phys. Rev. B, 81, 075414(2010).

    [16] J.-S. Huang, J. Kern, P. Geisler, P. Weinmann, M. Kamp, A. Forchel, P. Biagioni, B. Hecht. Mode imaging and selection in strongly coupled nanoantennas. Nano Lett., 10, 2105-2110(2010).

    [17] A. Kaplan, M. Tomes, T. Carmon, M. Kozlov, O. Cohen, G. Bartal, H. G. L. Schwefel. Finite element simulation of a perturbed axial-symmetric whispering-gallery mode and its use for intensity enhancement with a nanoparticle coupled to a microtoroid. Opt. Express, 21, 14169-14180(2013).

    [18] M. L. Gorodetsky, A. A. Savchenkov, V. S. Ilchenko. Ultimate Q of optical microsphere resonators. Opt. Lett., 21, 453-455(1996).

    [19] T. J. A. Kippenberg. Nonlinear optics in ultra-high Q whispering-gallery optical microcavities(2004).

    [20] R. A. Waldron. Perturbation theory of resonant cavities. Proc. IEE, 107, 272-274(1960).

    [21] H. A. Haus. Waves and Fields in Optoelectronics(1984).

    [22] S. Arnold, M. Khoshsima, I. Teraoka, S. Holler, F. Vollmer. Shift of whispering-gallery modes in microspheres by protein adsorption. Opt. Lett., 28, 272-274(2003).

    [23] L. Novotny, B. Hecht. Principles of Nano-Optics(2012).

    [24] H. Chen, L. Shao, Q. Li, J. Wang. Gold nanorods and their plasmonic properties. Chem. Soc. Rev., 42, 2679-2724(2013).

    [25] D. P. Sprünken, H. Omi, K. Furukawa, H. Nakashima, I. Sychugov, Y. Kobayashi, K. Torimitsu. Influence of the local environment on determining aspect-ratio distributions of gold nanorods in solution using Gans theory. J. Phys. Chem. C, 111, 14299-14306(2007).

    [26] C. F. Bohren, D. R. Huffman. Absorption and Scattering of Light by Small Particles(2008).

    [27] S. W. Prescott, P. Mulvaney. Gold nanorod extinction spectra. J. Appl. Phys., 99, 123504(2006).

    [28] F. Ruesink, H. M. Doeleman, E. Verhagen, A. F. Koenderink. Controlling nanoantenna polarizability through backaction via a single cavity mode. Phys. Rev. Lett., 120, 206101(2018).

    [29] L. Novotny. Effective wavelength scaling for optical antennas. Phys. Rev. Lett., 98, 266802(2007).

    [30] V. R. Dantham, S. Holler, C. Barbre, D. Keng, V. Kolchenko, S. Arnold. Label-free detection of single protein using a nanoplasmonic-photonic hybrid microcavity. Nano Lett., 13, 3347-3351(2013).

    [31] B. E. Little, S. T. Chu, H. A. Haus. Microring resonator channel dropping filters. J. Lightwave Technol., 15, 998-1005(1997).

    [32] W. Suh, Z. Wang, S. Fan. Temporal coupled-mode theory and the presence of non-orthogonal modes in lossless multimode cavities. IEEE J. Quantum Electron., 40, 1511-1518(2004).

    [33] A. Yariv. Universal relations for coupling of optical power between microresonators and dielectric waveguides. Electron. Lett., 36, 321-322(2000).

    [34] J. Alegret, T. Rindzevicius, T. Pakizeh, Y. Alaverdyan, L. Gunnarsson, M. Käll. Plasmonic properties of silver trimers with trigonal symmetry fabricated by electron-beam lithography. J. Phys. Chem. C, 112, 14313-14317(2008).

    [35] S. L. Teo, V. K. Lin, R. Marty, N. Large, E. A. Llado, A. Arbouet, C. Girard, J. Aizpurua, S. Tripathy, A. Mlayah. Gold nanoring trimers: a versatile structure for infrared sensing. Opt. Express, 18, 22271-22282(2010).

    [36] T. J. Davis, D. E. Gómez, K. C. Vernon. Simple model for the hybridization of surface plasmon resonances in metallic nanoparticles. Nano Lett., 10, 2618-2625(2010).

    [37] A. M. Funston, T. J. Davis, C. Novo, P. Mulvaney. Coupling modes of gold trimer superstructures. Philos. Trans. R. Soc. A, 369, 3472-3482(2011).

    [38] T. J. Davis, K. C. Vernon, D. E. Gómez. Designing plasmonic systems using optical coupling between nanoparticles. Phys. Rev. B, 79, 155423(2009).

    [39] A. Ashkin, J. M. Dziedzic, J. E. Bjorkholm, S. Chu. Observation of a single-beam gradient force optical trap for dielectric particles. Opt. Lett., 11, 288-290(1986).

    [40] J. E. Melzer, E. McLeod. Fundamental limits of optical tweezer nanoparticle manipulation speeds. ACS Nano, 12, 2440-2447(2018).

    [41] P. B. Johnson, R. W. Christy. Optical constants of the noble metals. Phys. Rev. B, 6, 4370-4379(1972).

    [42] M. Oxborrow. Traceable 2-D finite-element simulation of the whispering-gallery modes of axisymmetric electromagnetic resonators. IEEE Trans. Microwave Theory Tech., 55, 1209-1218(2007).

    [43] J.-P. Berenger. A perfectly matched layer for the absorption of electromagnetic waves. J. Comput. Phys., 114, 185-200(1994).

    [44] B. Min, L. Yang, K. Vahala. Perturbative analytic theory of an ultrahigh-Q toroidal microcavity. Phys. Rev. A, 76, 013823(2007).

    [45] M. Oxborrow. Configuration of COMSOL Multiphysics for simulating axisymmetric dielectric resonators: explicit weak-form expressions axisymmetric electromagnetic resonators.

    [46] J. Zhu, S. K. Ozdemir, Y.-F. Xiao, L. Li, L. He, D.-R. Chen, L. Yang. On-chip single nanoparticle detection and sizing by mode splitting in an ultrahigh-Q microresonator. Nat. Photonics, 4, 46-49(2010).

    Cheng Li, Lei Chen, Euan McLeod, Judith Su, "Dark mode plasmonic optical microcavity biochemical sensor," Photonics Res. 7, 939 (2019)
    Download Citation