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 () 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 with an electric-field intensity enhancement of 4300, which far exceeds that of a single rod ( times). Through an overall evaluation of a combined enhancement factor, which includes the factor of the system, the sensitivity of the trimer system was improved versus 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 () factors and their ability to be integrated on chip [1–9]. These sensors when combined with frequency-locking, balanced detection, and data processing techniques are able to detect wavelength shifts () 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 factor of the resonance, that is, the broadening of the linewidth. Dark mode plasmonic resonances can help solve this problem [12–16], but simulation is necessary to understand the specific impact of a given plasmonic structure on the system.
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 []. (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 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).
Sign up for Photonics Research TOC. Get the latest issue of Photonics Research delivered right to you!Sign up now
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 factor mentioned above [18], which determines the linewidth of the resonance peak in the output spectrum. Higher factors correspond to narrower resonances, whose central resonance frequency can be more precisely determined and tracked [1–3,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 with (possibly anisotropic) material parameters () and () relative to the background is given by [20–22] where and are the permittivity and permeability of free space, and are relative permittivities and permeabilities, is the permittivity of the background media (in this case, vacuum), and are the E and magnetic fields for an unperturbed WGM, represents conjugate-transpose, and are the electric and magnetic fields within the perturbation volume, is the maximum intensity within the unperturbed toroid, is the mode volume of the resonant mode [19], is the dielectric constant of silica [23], and is the polarizability tensor of the particle. In the numerator, we have assumed , 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 [24–26] 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 is the volume of the cylindrical nanorod, is the relative permittivity of the metal, is the dielectric constant of the surrounding environment, and the are geometrical factors that can be, respectively, written as where is the eccentricity of the prolate spheroid, and is a parameter that describes the particle shape. Here, we use 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 () 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).
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) 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 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 factor of the system, even with strong local E-field enhancement, so there is a trade-off between and the intensity enhancement factor, as shown in Fig. 2(c). At a wavelength of , the maximum 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 factor of , which drops to when coupled to the nanorod with an aspect ratio 5.4.
The lower 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 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: , for scalar . Then, the volume of the smallest detectable particle follows, , assuming that the smallest experimentally measurable shift, , is inversely proportional to the factor because high- 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 -factor effects, which we call the combined enhancement factor (): where and represent the factors of the coupled and bare toroid cavities, respectively. Higher values correspond to better sensitivity and/or more well-defined resonances. Unlike the peak in , the highest 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 field enhancement factor and a factor of . Using 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,31–33], the equations of several coupled modes can be written as where and represent the mode amplitude of the WGM mode and LSPR mode of the th nanorod, and correspond to their resonance frequencies, and and 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 . 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 , where , and The intrinsic factor of the coupled system can be calculated from the eigenvalues of , corresponding to the steady state of the system. Through the simulation data, the coupling coefficient () can be obtained from a matrix generated by a system in which a cavity is coupled to a single rod (see Appendix B). Further, the factor of a matrix corresponding to a multiple () rod coupled system can be calculated numerically. We compare these values with the 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.
Figure 3.(a) 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 , Eq. (9)].
4. TRIMER-INDUCED LATERAL DARK MODE
The dipole modes induced in isolated antennas lead to large energy radiation, which adversely limit the 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,36–38]. Coupling of the dark mode and WGM standing wave can significantly reduce degradation of the 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.
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 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 () is much larger than that of a single rod (), as mentioned above. At the same time, the factor of is similar to that of the single rod case () at on-resonance coupling and is twice as large as the for the three rod case () 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 of the system is , which is achieved using the trimer with a rod aspect ratio of 4.5. This is also significantly larger than that of a single rod ().
To test the robustness of our system, we studied the effect of tilted angle perturbation, length perturbation, and inter-rod spacing on . 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 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 factor or enhancement factor of the system, the overall 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 value is almost unchanged, and the observed variation is primarily due to E-field enhancement.
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 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 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 . When the rotation angle is greater than 45 deg, the decreases as the rotation angle increases [Fig. 5(d)] due to mismatch with the polarization direction of the WGM. The system has only 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 at a 30 deg angle is , which is only 11.5% lower than the maximum enhancement of . 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 combined enhancement factor versus 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 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 , 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 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 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 [42–45]. The factor is calculated from the real and imaginary parts of the eigenfrequency obtained by COMSOL. In the matrix given in Eq. (5), the parameters () and () are determined from simulations of a bare toroid and the interaction of a single rod with a cavity using the eigenfrequency solver in COMSOL. () and () 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 value. The 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 and can be written as where is an integer corresponding to the number of antinodes of the standing wave in Fig. 6(b). satisfies the resonance condition , where represents the effective refractive index of the mode. The 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 case is shown in Figs. 6(b) and 6(c).
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 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 . 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
and 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 , which can be solved from the simulation of the bare toroid coupled to a single nanorod, then from the analytical eigenvalue of the matrix, we can calculate using