1. INTRODUCTION
High-index dielectric nanoparticles supporting distinct Mie resonances have aroused widespread research interest due to their excellent light manipulation performance at the subwavelength scale [1]. In 2010, resonance-enhanced light scattering from silicon (Si) nanowires was first observed in the visible and near-infrared spectral range [2,3]. Since then, the research on strong electric and magnetic responses and their interactions in high-index nanoparticles has revealed rich resonance-induced phenomena including structural color [4–6], directional scattering [7–10], scattering suppression by an anapole mode [11], significant enhancement of nonlinear optical response [12,13], and so on. Compared with metallic nanoparticles, dielectric nanoparticles support magnetic resonances and possess much lower optical loss. These advantages lead to the strongly localized electric field inside dielectric nanoparticles at optical resonances, greatly improving light-matter interaction.
Increasing quality factors ( factors) of resonant modes is vital for generating remarkable near-field enhancement [14–16] and achieving strong light-matter interaction [17]. In 2012, it was found that a dielectric nanoparticle placed on a gold (Au) substrate can interact strongly with its mirror image [18]. A further study revealed that the cooperation of the electric dipole (ED) excited in the nanoparticle and its virtual counterpart can be interpreted as a mirror-image-induced magnetic dipole (MMD), leading to enhanced scattering [19]. In addition, the factor of the magnetic dipole (MD) resonance supported by the Si/Au nanocavity can be effectively improved by optimizing the polarization and angle of the incident light [20]. For a nanocavity composed of a Si nanoparticle, a silica () spacer, and a silver (Ag) substrate, the resonant strength and linewidth of the ED and MD modes can be manipulated by varying the thickness of the spacer layer. It was demonstrated that the thin spacer greatly enhanced the MD resonance in the Si nanoparticle and induced phonon-assisted hot luminescence [21]. Besides, when metals are introduced into a system, it is important to reduce Ohmic loss [22]. The spacer separates the Si nanoparticle from the Ag substrate, which effectively suppresses the effect of Ohmic loss on the optical modes in the Si nanoparticle. Therefore, a nanocavity composed of a high-index nanoparticle and a dielectric-metal heterostructure acts as a good platform for realizing the enhancement of light-matter interaction at the subwavelength scale.
As the core material for optoelectric devices, Si has excellent characteristics of high refractive index and low material loss in the visible and near-infrared bands. Light emitters based on Si nanostructures have always attracted particular interest as they provide unique opportunities for highly integrated light sources for photonic circuits [23]. An on-chip Si light source could potentially achieve a higher integration density with a compact size and display better performance in terms of energy efficiency [24]. However, the quantum efficiency of bulk Si is limited by indirect bandgap structure [25]. Hence, resonance-enhanced light-matter interaction is essential to improve the luminescence performance of Si nanoparticles or nanostructures. In 2018, white light emission from Si nanoparticles placed on a substrate was demonstrated by resonantly exciting their MD resonances with femtosecond laser pulses [26]. The luminescence is caused by the interband transition of hot carriers enhanced by the Auger effect and Purcell effect, which is different from that caused by the intraband transition of hot carriers in GaAs [27]. The efficiency of the white light emission was further enhanced by improving the factor of the resonant mode at the excitation wavelength, for example, adopting the MMD resonance of a Si nanoparticle placed on a metal substrate [28] or quasi-bound states in the continuum of a Si nano cuboid [29] and a Si metasurface [30]. These strategies can substantially reduce the threshold for the luminescence burst. In addition, it was reported that the MD resonance of a Si nanoparticle can be effectively enhanced by adding a spacer in between the nanoparticle and an Ag substrate [21]. Although it has been known that the optical modes supported by a Si nanoparticle can be modified by a dielectric-metal heterostructure, a systematic study on the linear and nonlinear optical properties of a hybrid nanocavity composed of a Si nanoparticle and a dielectric-metal heterostructure is still lacking.
Sign up for Photonics Research TOC Get the latest issue of Advanced Photonics delivered right to you!Sign up now
Mie resonances with high factors can not only enhance the interaction between the optical mode and the cavity itself but also enhance the interaction between light and two-dimensional (2D) materials. A dielectric-metal heterostructure can support transverse-electric-polarized (TE-polarized) surface-bound waves [31] under the physical mechanism of optical magnetism [32]. The TE wave, which is dominated by the in-plane electric field, can fully interact with 2D material embedded into the heterostructure. For example, a narrow linewidth TE wave supported by a silicon nitride/silver () heterostructure can interact with a tungsten disulfide () monolayer to achieve strong exciton-photon coupling [33]. However, being a surface-bound wave, the confinement of the electromagnetic field in the propagation direction of the TE wave is still lacking. Therefore, it is necessary to introduce a high-index dielectric nanoparticle on the surface of the heterostructure to improve local field enhancement for further enhancing light-matter interaction.
In this paper, we investigate systematically the enhanced light-matter interaction based on hybrid nanocavities. Both low-threshold luminescence burst and strong exciton-photon coupling are demonstrated. First, the enhanced MD resonance in the hybrid nanocavity substantially improves the nonlinear-optical-absorption-induced interband transition of hot electrons. With the help of the Auger effect [34] and Purcell effect, luminescence burst from the nanocavity is stimulated by femtosecond laser pulses with low energy of . Compared with the case of a Si nanoparticle placed on a substrate, not only the threshold pulse energy is reduced by 30.8 pJ, but also the luminescence intensity is increased by one order of magnitude. Second, a monolayer is placed in between the Si nanoparticle and heterostructure to achieve strong exciton-photon coupling with the excitation of s-polarized light. By changing the incident angle of s-polarized light, the resonant wavelength of the optical mode can be adjusted covering the two exciton resonances of the monolayer. Large Rabi splitting energies of and are obtained at the A- and B-exciton ( and ) resonances, respectively.
2. RESULTS AND DISCUSSION
We first examine the effect of a heterostructure on the Mie resonances of a Si nanoparticle based on numerical simulation. As schematically shown in Fig. 1(a), a Si nanoparticle with a diameter of is placed on a substrate. The thicknesses of the Ag film and spacer are noted as and , respectively. The Si nanoparticle is illuminated by a broadband light source propagating along the direction and the backward scattering is collected by an objective. We performed the simulation with the commercial software Ansys Lumerical based on the finite-difference time-domain (FDTD) method. In Fig. 1(b), we show the evolution of the scattering spectrum of a nanocavity () when the thickness of the Ag film is fixed () and that of the spacer layer is increased. It is noticed that the intensity of the optical resonance located at is enhanced significantly when is in the range from 70 to 95 nm. In addition, it is found that the intensity of the optical resonance located around 600 nm begins to increase when rises beyond 95 nm. These behaviors imply that the optical resonances supported by the nanocavity can be effectively manipulated by simply varying the thickness of the spacer layer. Based on multipolar decomposition (see Appendix A), we reveal that the optical resonances appearing around 600 nm and 765 nm are dominated by ED and MD, respectively. The scattering spectra of the nanocavity with different are discussed in Appendix B, which indicates that a 50 nm Ag layer is sufficient to be a good reflector. We examined the effects of a heterostructure (, ) on the scattering properties of a Si nanoparticle. As shown in Figs. 1(c)–1(e), we compared the scattering spectra of Si nanoparticles with different diameters on a substrate with and without the heterostructure. The results reveal that increasing the diameter of the Si nanoparticle primarily increases the resonant wavelengths of both ED and MD resonances. In addition, it is noticed that the factors of the ED and MD resonances are improved by using the heterostructure. The factors of the ED and MD resonant modes can be calculated by the formula , where is the peak wavelength in the scattering spectrum of each mode and is the full width at half maximum (FWHM) of the scattering spectrum. Gaussian fitting functions were used to extract of each resonant mode from the scattering spectrum. Without loss of generality, taking the result in Fig. 1(d) as an example, the factors of ED and MD increase from 7.2 and 8.2 to 18.6 and 11.4, respectively, after the introduction of the heterostructure. The improvement of factor results from the destructive interference of the in-plane ED or MD with its mirror mode reflected from the Ag film. Here the heterostructure acts as a Fabry–Perot cavity. The degree of interference cancellation of the out-of-plane radiation can be tuned by varying the thickness of the spacer [16,35].

Figure 1.Optical resonances supported by a nanocavity. (a) Schematic configuration of a hybrid nanocavity. (b) Calculated scattering spectra (normalized to the maximum value) of nanocavities (, ) with variable spacer layer thickness. (c)–(e) Comparison of the scattering spectra (each spectrum is normalized to its maximum value) of Si nanoparticles with different diameters placed on a substrate (, ) (red curves) and a substrate (blue curves).
The energy band structure of Si and the carrier dynamics in a Si nanoparticle excited by femtosecond laser pulses are depicted schematically in Fig. 2(a). Since the energy gap between the conduction and valence bands at the point is 3.4 eV, the interband transition of electrons is initiated by two- or three-photon-induced absorption (2PA or 3PA) of femtosecond laser pulses. The enhancement factors of 2PA and 3PA ( and ) are proportional to the fourth and sixth power of the electric field inside the Si nanoparticle, respectively, i.e., , , where is the amplitude of the incident electric field and is the volume of the Si nanoparticle. Therefore, it is expected that the 2PA and 3PA of a Si nanoparticle will be effectively enhanced at the ED and MD resonances, generating high-density carriers in the Si nanoparticle. Owing to the Auger effect, the non-radiative recombination lifetime () of hot electrons can be significantly increased (from 0.1–1.0 ps to 10–100 ps). Theoretically, the fluorescence quantum efficiency () of a Si nanoparticle can be described by the formula , where represents the radiative recombination lifetime of hot electrons at the emission wavelength. Relying on the Purcell effect, the radiative recombination lifetime can be reduced by exploiting the high-order Mie resonances of the Si nanoparticle, such as the electric quadrupole (EQ) and magnetic quadrupole (MQ) resonances. Therefore, the combination of the Auger effect and the Purcell effect can dramatically enhance the nonlinear photoluminescence (PL) of a Si nanoparticle.

Figure 2.Physical mechanism for the enhanced nonlinear optical emission from a nanocavity. (a) Energy band diagram of Si and carrier dynamics in a Si nanoparticle excited by femtosecond laser pulses. (b) Calculated scattering and 2PA/3PA spectra for a Si nanoparticle () placed on a substrate. The scattering spectrum is normalized to its maximum value. (c), (d) Electric field amplitude distributions in the hybrid nanocavity calculated at the ED and MD resonances.
Figure 2(b) illustrates the scattering and 2PA/3PA spectra calculated for the hybrid nanocavity (, , ). One can see that the enhancement factors of 2PA and 3PA reach peak values at the ED () and MD () resonances of the nanocavity. The electric field distributions at the ED and MD resonances are shown in Figs. 2(c) and 2(d), respectively. Remarkably, the electric field is strongly localized in the Si nanoparticle at the MD resonance, leading to larger enhancement factors of 2PA and 3PA than those at the ED resonance. Compared with the Si nanoparticle placed on a substrate (see Appendix C), the enhancement factors of 2PA and 3PA are improved by and times, respectively.
We studied the nonlinear optical emission of a hybrid nanocavity (, , ) under the excitation of femtosecond laser pulses. The detailed experimental setup is described in Appendix D. In Fig. 3(a), we show the evolution of the PL spectrum of the nanocavity with increasing pumping energy. The excitation wavelength was set at the MD resonance of the nanocavity () based on the scattering spectrum shown in Fig. 3(b). In each case, the emission from the nanocavity was recorded by using a charge-coupled device (CCD), as shown in the insets of Fig. 3(a). When the pulse energy was increased from 3.9 to 15.8 pJ, the hot electron luminescence from the nanocavity increased only slightly. Notably, a burst of white light emission from the nanocavity was observed when the pulse energy exceeded 20.5 pJ. The cooperation of the intrinsic excitation and the Auger effect leads to a large number of hot electrons in the high-energy states (around the point), which significantly enhances the nonlinear optical emission. In Fig. 3(b), we compared the scattering spectra of the nanocavity before and after the luminescence burst. It is noticed that the scattering intensities of the nanocavity at the MD and ED resonances are greatly enhanced after the luminescence burst. This phenomenon can be ascribed to the annealing of the Si nanoparticle by the high temperature induced by the laser pulses, improving the quality of the crystalline structure. In addition, the scattering peaks of MD and ED are blue-shifted and red-shifted, respectively. It is caused by the giant optical force applied by femtosecond laser pulses on the Si nanoparticle. Under the optical force, the nanoparticle is squeezed in the out-of-plane direction and stretched in the plane along the polarization direction of the incident light [36]. In addition, the MQ resonance at is excited by the large-angle components of the broadband incident light focused by using an objective with a numerical aperture (NA) of 0.8 (see Appendix A). In Fig. 3(c), we present the dependence of the integrated PL intensity on the excitation pulse energy observed for the nanocavity. It can be seen that the PL intensity increases exponentially when the excitation pulse energy exceeds 20.5 pJ. For comparison, we also examined the nonlinear optical emission of a Si nanoparticle placed on a substrate. Based on the scattering spectra shown in Fig. 3(e), the excitation wavelength was chosen at 745 nm so that the MD resonance was resonantly excited. As shown in Fig. 3(d), the luminescence burst could be observed only when the excitation pulse energy was raised to 51.3 pJ. This behavior is also shown in Fig. 3(f) where the dependence of the integrated PL intensity on the pumping pulse energy is presented. By comparing the results shown in Figs. 3(c) and 3(f), it can be seen that the threshold for the luminescence burst is greatly reduced in the nanocavity when the heterostructure is used as the substrate. More importantly, it is found that the PL intensity of the nanocavity is one order of magnitude larger than that of the Si nanoparticle on the substrate.

Figure 3.Nonlinear optical emission from a nanocavity under the excitation of femtosecond laser pulses. (a) PL spectra of the nanocavity measured at different pumping pulse energies. The corresponding CCD images of the emission from the nanocavity are shown in the insets. (b) Backward scattering spectra (normalized to the maximum value) measured for the nanocavity before and after the luminescence burst. The scanning electron microscope (SEM) image of the Si nanoparticle before the excitation is shown in the inset. (c) Dependence of the integrated PL intensity on the pumping pulse energy observed for the nanocavity. (d) PL spectra of a Si nanoparticle placed on a substrate measured at different pumping pulse energies. (e) Forward scattering spectra (normalized to the maximum value) measured for the Si nanoparticle placed on the substrate before and after the luminescence burst. (f) Dependence of the integrated PL intensity on the pumping pulse energy observed for the Si nanoparticle placed on the substrate.
The above results indicate clearly that the interaction between the Si nanoparticle and the femtosecond laser pulses can be effectively enhanced by the modified MD resonance supported by the nanocavity, which possessed a larger factor. As a result, a low pulse energy threshold () for the luminescence burst is achieved under the excitation of 750 nm femtosecond laser pulses. It should be emphasized that the generation of hot electrons in the conduction band is dominated by the 3PA process when the excitation wavelength is chosen at 750 nm. In this case, the threshold for the luminescence burst is larger than that reported in the previous study where the use of 720 nm femtosecond laser pulses implies that the 2PA and 3PA processes contribute simultaneously to the generation of hot electrons [28].
Recently, it was found that a heterostructure can support TE waves with significantly reduced damping rates as compared with conventional surface plasmon polaritons (TM wave). More interestingly, the electric field of a TE wave is localized on the surface of the layer [31], which is beneficial for the interaction with a 2D material attached to the heterostructure [33]. A further enhancement between the TE wave and the 2D material is expected if a Si nanoparticle is placed on the heterostructure, forming a nanocavity. In this case, the lateral confinement of the electric field in the nanocavity leads to the electric field enhancement at the MD, ED, and MQ resonances supported by the nanocavity.
As shown schematically in Fig. 4(a), a nanocavity is created by placing a Si nanoparticle on a heterostructure. The nanocavity can be excited by the TE wave propagating on the surface of the heterostructure stimulated by s-polarized incident light in the Kretschmann-Raether configuration [37]. In the experiment, we chose a Si nanoparticle with a diameter of . The SEM image of the Si nanoparticle is shown in the inset of Fig. 4(b). We first examined the scattering spectrum of the nanocavity under normal illumination, as shown in Fig. 4(b). We could identify the MQ, ED, and MD resonances supported by the nanocavity, which appear at the wavelengths of , , and , respectively. Then, we examined the scattering of the nanocavity excited by s-polarized light with different incident angles, which can generate TE waves on the surface of the heterostructure. The scattering spectra of the nanocavity measured at different incident angles are shown in Fig. 4(c). The actual incident angle at the interface is calculated considering the refraction caused by the prism () and the () substrate. In each case, one can identify two scattering peaks in the scattering spectrum. One is fixed at while the other one has a blue shift with the increase of incident angle. It has been known that the wavelength of the TE wave decreases when the incident angle is increased. On the other hand, the Si nanoparticle supports three optical resonances located at . , and , corresponding to the MQ, ED, and MD resonances [denoted by shadows in Fig. 4(c)]. These two features are reflected in the scattering spectra of the nanocavity. The intensity of the TE wave is enhanced when its wavelength coincides with the optical resonances supported by the Si nanoparticle. The MD resonance can be found in all scattering spectra. However, its intensity is dramatically reduced at large incident angles when the wavelength of the TE wave is shifted to shorter wavelengths. It is expected that the electric field will be enhanced in the nanocavity due to the coupling between the TE wave and the Mie resonances.

Figure 4.Coupling between the TE wave supported by a heterostructure and the Mie resonances supported by a Si nanoparticle in a nanocavity. (a) Schematic showing the excitation of a nanocavity by using the TE wave propagating on the surface of the heterostructure stimulated by s-polarized incident light. The electric component () of TE wave is concentrated on the surface of the spacer, while the corresponding magnetic component () mainly exists at the interface. (b) Backward scattering spectrum of the nanocavity. The SEM image of the Si nanoparticle is shown in the inset. (c) Measured scattering spectra of the nanocavity excited by s-polarized light with different incident angles. The corresponding CCD images of the scattering are shown in the insets. The wavelength bands where MQ, ED, and MD exist are denoted by green, yellow, and red shadows, respectively. The angular dispersion of the TE wave is represented by a dashed curve, schematically. (d) Simulated scattering spectra (in a relative scale) of the nanocavity excited by s-polarized light with different incident angles.
In order to gain a deep insight into the scattering properties of the nanocavity, we calculated the scattering spectra of the nanocavity at different incident angles by using the FDTD method, as shown in Fig. 4(d). As the incident angle increases from 44° to 62°, the resonant wavelength of the TE wave has a blue shift from to . Remarkably, when the wavelength of TE wave overlaps with the intrinsic wavelength of the Mie resonances of the Si nanoparticle, scattering is enhanced and vice versa. The strongest scattering intensity is observed at , where the wavelength of the TE wave coincides with the ED resonance. It implies that the TE wave can be effectively coupled into the ED resonance. In the previous work, a polystyrene (PS) nanoparticle was used as the probe for the TE wave [33]; however, the intensity of the TE wave is only slightly modified due to the low refractive index of PS (see Appendix E). In contrast, the use of a high-index Si nanoparticle in this work can not only act as a probe but also achieve strong near-field enhancement. Therefore, the interaction between the TE wave and a 2D material will be enhanced in the nanocavity.
To verify the above expectation, we inserted a monolayer into the nanocavity and examined the exciton-photon interaction, as schematically shown in Fig. 5(a). The SEM and forward scattering images of the nanocavities with an embedded monolayer are shown in the insets. In the experiment, the region in which the spectrometer collects the scattering spectrum can be limited to surrounding the Si nanoparticle, thus eliminating the influence of other scattered signals. We chose a Si nanoparticle with and measured the angle-resolved scattering spectra of the nanocavity with an embedded monolayer, as shown in Figs. 5(b) and 5(c). The CCD images of the scattered light are shown in the insets. By comparing the spectra with the results in Fig. 4(c), it can be found that the refractive index and material loss introduced by the monolayer make the position of the resonant peak shift towards the longer wavelength at the same excitation angle. When the incident angle increased from 62.9° to 65.8°, a split of the scattering peak is observed, implying the coupling between the TE wave and of the monolayer. With the decreasing incident angle, the wavelength of the TE wave is shifted to longer wavelengths. Similarly, one can see a split of the scattering peak when the TE wave is moved to . It suggests the coupling between the TE wave and the of the monolayer.

Figure 5.Exciton-photon coupling in a nanocavity with an embedded monolayer. (a) Schematic showing the excitation of a nanocavity with an embedded monolayer and the detection of the scattered light. The SEM and forward scattering images of nanocavities with an embedded monolayer are shown in the insets. (b), (c) Angle-resolved scattering spectra of the nanocavity () with an embedded monolayer. Each of the spectra is normalized to its maximum value. The CCD images of the scattered light are shown in the insets. (d), (e) Dispersion relations of the upper and lower polariton branches (purple dots) extracted from the angle-resolved scattering spectra of the nanocavity with an embedded monolayer around the resonant wavelengths of the and . The red and blue curves are the fittings of the dispersion relations based on the coupled harmonic oscillator model. The energies of the excitons and photons (TE wave) are indicated by dashed lines.
Based on the angle-resolved scattering spectra measured for the nanocavity with an embedded monolayer, we could extract the dispersion relation of the lower and upper polaritons formed by the coupling between the TE wave and the two excitons of the monolayer. The results for the and are shown in Figs. 5(d) and 5(e), respectively. The criterion for the strong coupling of optical mode and exciton is , where is the Rabi splitting energy, and and are the linewidths of the TE wave and exciton resonance, respectively [38,39]. Theoretically, the linewidths of A- and B-excitons can be derived from the imaginary part of the relative permittivity of the monolayer [40]. In the analytical discussion, we use the theoretical linewidths of A- and B-excitons, which are 33 and 81 meV, respectively [33]. The linewidths of the TE wave obtained from Fig. 4(d) are 91 and 65 meV, respectively, at the wavelengths of A- and B-excitons. As shown in Fig. 5(d), the Rabi splitting derived from the lower and upper polariton branches at the was estimated to be , satisfying the strong coupling criterion. In addition, this value is larger than those observed in Si/Au nanocavities () [41] and nanocavities () [33]. As discussed above, the intensity of the TE wave is greatly enhanced at by the ED resonance supported by the Si nanoparticle. Consequently, the in-plane electric field is enhanced due to the strong interaction between the TE wave and the ED resonance, leading to the enhanced exciton-photon coupling observed at the resonant wavelength of the . Similarly, shown in Fig. 5(e), the Rabi splitting at the resonant wavelength of the was derived to be from the fitting of the lower and upper polariton branches. The strong coupling condition is also fulfilled in this case. The corresponding simulation results are shown in Fig. 12 of Appendix F. The calculated Rabi splitting energy () for A-exciton agrees well with the measured value (). However, the Rabi splitting energy obtained by simulation is , which is twice the measured one (). By reviewing the results in Fig. 4(c), we can see that the scattering is very weak near the wavelength of the B-exciton. At the same time, the radiation intensity of the B-exciton is much weaker than that of the A-exciton. These two reasons cause the scattered signal to be easily covered by the noise in the experiment. The signal noise caused by the divergence of the light source and the irregular shape of the Si nanoparticles also lead to the deviation of the experiment compared with the simulation results. In addition, a recent study revealed that the dissipation of coupled subsystems, such as optical modes and exciton, will cause the deviation between measured spectral splitting and theoretical eigenlevel splitting [42]. It enriches the understanding of strong light-matter interactions at room temperature.
3. CONCLUSION
In conclusion, based on nanocavities, we investigated numerically and experimentally the enhanced interaction between femtosecond laser pulses and Si nanoparticles and that between a TE wave and a monolayer. First, it was revealed that a nanocavity can support an MD resonance with a larger factor than that supported by a Si nanoparticle placed on a substrate. By resonantly exciting the MD resonance of the nanocavity with femtosecond laser pulses, we observed the luminescence burst from the Si nanoparticle with a threshold of . Compared with a Si nanoparticle placed on a substrate, the threshold for the luminescence burst achieved in the nanocavity is reduced by more than a half while the PL intensity is increased by one order of magnitude. Second, the coupling between the TE waves generated on the surface of the heterostructure and the Mie resonances supported by the Si nanoparticle can effectively enhance the in-plane electric field, which leads to an enhanced exciton-photon coupling in a nanocavity with an embedded monolayer. Relying on the electric field enhancement achieved at the ED and MQ resonances of the nanocavity, Rabi splitting energies as large as and were observed at the and of the monolayer. It indicates an enhanced exciton-photon coupling at the as compared with nanocavities. The hybrid nanocavity proposed in this work provides a useful platform for realizing the strong interaction of light with nanoparticles and 2D materials, which is beneficial for achieving nanoscale light sources.
APPENDIX A: MULTIPOLAR DECOMPOSITION
Multipolar decomposition is a key tool for analyzing the properties of resonant modes. Here, we adopt exact expressions for ED, MD, EQ, and MQ moments of dynamic sources. First, the electromagnetic field distributions of resonant modes in the nanocavity excited by plane waves are calculated by using Ansys Lumerical based on the FDTD method. Then, the scattering cross section of multipolar moments is calculated by using the formulas proposed in Ref. [43], and the results are scaled by the surface area of the Si nanoparticle. Figure 6(a) shows the result of multipolar decomposition for the resonances excited by a plane wave with normal incidence. We can see that the resonant modes that exist around the wavelength of 600 nm and 765 nm are dominated by ED and MD, respectively. The scattering cross section of EQ and MQ remains at very low levels throughout the whole spectrum. In addition, the large-angle components of the incident light produced by an objective () need to be considered. Figure 6(b) shows the result for the simultaneous incidence of s-polarized and p-polarized plane waves at the incident angle of 53°. Compared with the case of normal incidence, the spectra of ED and MD are broadened accompanied by decreased amplitude. Notably, the obliquely incident light can stimulate the MQ resonance at the wavelength of ; this is consistent with the experimental result in Fig. 3(b).

Figure 6.Multipolar decomposition for resonances of the nanocavity (, ). The dashed curve represents the sum of the scattering cross section of the four multipolar components. (a) The nanocavity is excited by a plane wave with normal incidence. (b) The nanocavity is excited by the simultaneous incidence of s-polarized and p-polarized plane waves with an incident angle of 53°.
APPENDIX B: SCATTERING SPECTRA OF HYBRID NANOCAVITIES COMPOSED OF SILVER FILMS WITH DIFFERENT THICKNESSES
Here, we numerically analyze the effect of the Ag layer on the resonance of the nanocavity. Figure 7 shows the scattering spectra with different thicknesses () of the Ag layer. It can be seen that the scattering cross section at the wavelength of rises significantly when increases from 30 to 50 nm. When the thickness exceeds 50 nm, further increase in thickness has little effect on the scattering spectrum. It confirms that the 50 nm Ag layer is sufficient for the nanocavity to support strong resonances.

Figure 7.Simulated scattering spectra (in a relative scale) of nanocavities (, ) with different thicknesses of the Ag layer. The excitation light is a plane wave with normal incidence.
APPENDIX C: SCATTERING AND FIELD ENHANCEMENT OF THE Si NANOPARTICLE PLACED ON A SILICA SUBSTRATE
The simulated scattering and 2PA/3PA spectra of the Si nanoparticle placed on a substrate are shown in Fig. 8(a). It can be seen from the scattering spectrum that the strongest ED and MD resonances exist at the wavelengths of and , respectively. Although, the ED and MD resonances improve the enhancement factor of 2PA and 3PA, the peak values of the enhancement factor are far less than those of the nanocavity. The main reason is that the ED and MD resonances with low factors supported by the Si nanoparticle placed on a substrate cannot generate strong field enhancements inside the Si nanoparticle. From Figs. 8(b) and 8(c), we can see that the strongest electric field intensity is on the surface of the nanoparticle rather than inside the nanoparticle. By comparing the results with the field distributions in Figs. 2(c) and 2(d), it can be found that the nanocavity can significantly enhance the light intensity inside the Si nanoparticle.

Figure 8.Simulated results of the Si nanoparticle () placed on a substrate. (a) Scattering and 2PA/3PA spectra. The scattering spectrum is normalized to its maximum value. (b), (c) Electric field amplitude distributions at the ED and MD resonances.
APPENDIX D: NUMERICAL MODELING, SAMPLE FABRICATION, AND OPTICAL CHARACTERIZATION
1. Numerical ModelingIn this work, the scattering spectra and field intensity distributions of hybrid nanocavities were calculated numerically by using Ansys Lumerical based on the FDTD method. In the simulation model, the Si nanoparticle was suspended 2 nm above the substrate surface to avoid point contact between the two geometries. The whole model adopted the four-level mesh accuracy built into Lumerical software. Within a cube region containing the Si nanoparticle, the mesh was further refined in three dimensions to a size of 2 nm. A total-field scattered-field (TFSF) source was used to solve scattering models. And, a perfectly matched layer boundary condition was employed to terminate the finite simulation region.
2. Sample FabricationSi nanoparticles with different diameters were fabricated by using femtosecond laser ablation. The 800 nm femtosecond laser pulses (Legend Elite, Coherent) with a duration of 90 fs and a repetition rate of 1 kHz were focused on the surface of a Si wafer immersed in deionized water using a lens with a focal length of 150 mm. The aqueous solution of Si nanoparticles was dropped onto the substrate and substrate and allowed to dry naturally, obtaining the desired hybrid nanocavities.
The heterostructure used in this work was fabricated by the following procedure. First, an Ag film with a thickness of 50 nm was plated on a substrate by using electron beam evaporation. Then, a layer with a thickness of 95 nm was deposited on the Ag film through high-frequency plasma-enhanced chemical vapor deposition (HF-PECVD) to form a heterostructure. The monolayers purchased from Sixcarbon Tech Shenzhen were synthesized on a Si substrate via chemical vapor deposition and then transferred onto the heterostructure, forming a heterostructure.
3. Characterization of WS2 MonolayerFigure 9(a) shows the PL spectrum of the sample under the excitation of a 488 nm laser with the power of 0.5 mW. A distinct single peak exists at the wavelength of corresponding to the direct excitonic transition for A-exciton; it rules out the possibility of multilayer [44]. In addition, we also measured the Raman spectrum of the sample under the excitation of the 633 nm laser shown in Fig. 9(b). According to a previous work [45], spectral splitting exists at mode in the Raman spectra of two-layer and three-layer samples. However, the measured Raman spectrum of our sample shows a clear single peak of mode at . This result further proves that the sample used in our experiment is a monolayer.

Figure 9.(a) PL spectrum of the sample placed on the heterostructure using a 488 nm excitation laser with a power of 0.5 mW. (b) Raman spectrum of using a 633 nm excitation laser.

Figure 10.Schematic diagram of experimental setups for measuring the nonlinear optical emission from the nanocavity under the excitation of femtosecond laser pulses.
The TE waves supported by the heterostructure were excited in the Kretschmann-Raether configuration, as shown in Fig. 4(a). The sample was mounted on a prism made of (K9 glass) with silicone oil applied, which has a refractive index close to that of . The TE wave was excited by utilizing the evanescent waves generated through total internal reflection at the surface of the prism. The scattering spectra of the hybrid nanocavities were characterized by using a dark-field optical microscope (Observer A1, Zeiss), with the choice of either broadband white light or s-polarized light for illumination.
APPENDIX E: SCATTERING OF A POLYSTYRENE NANOPARTICLE PLACED ON THE HETEROSTRUCTURE
The refractive index of the PS nanoparticle () is lower than that of the layer (). Thus, a PS nanoparticle placed on the heterostructure cannot support resonances with high factors. When the TE wave supported by the heterostructure meets the PS nanoparticle, it is only scattered by the PS nanoparticle without producing local field enhancement. Figure 11 shows the scattering spectra excited by s-polarized light with different incident angles. We can see that the peak wavelength of the scattered light has a blue shift with the increase of the incident angle. It is consistent with the dispersion characteristic of TE waves. In addition, the scattering peaks form a broad spectrum envelope, which indicates that the resonant modes supported by the PS nanoparticle have ultrahigh losses. Therefore, the PS nanoparticle can only act as a probe instead of being a nanocavity to achieve strong near-field enhancement.

Figure 11.Simulated scattering spectra (in a relative scale) of a PS nanoparticle () placed on the heterostructure excited by s-polarized light in the Kretschmann-Raether configuration with different incident angles.
APPENDIX F: SIMULATED SCATTERING SPECTRA OF THE HYBRID NANOCAVITY WITH AN EMBEDDED WS2 MONOLAYER
Figure 12 shows the angle resolved scattering spectra of the nanocavity with an embedded monolayer excited by s-polarized light when . When the incident angle is , the resonant wavelength of the TE wave overlaps with the wavelength of A-exciton, which is consistent with the experimental result of . In addition, the measured Rabi splitting energy at the A-exciton () is very close to that obtained by simulation (). For the coupling between the TE wave and B-exciton, the required incident angle in the experiment is around . However, this value is a bit larger than the angle () obtained in the simulation. The deviation in incident angle can be attributed to the divergence of an incident light beam. Especially, in the condition of large-angle oblique incidence, the TE wave interacting with the Si nanoparticle is mainly excited by the edge of the light spot, which makes the deviation more obvious. Besides, since the resonance of B-exciton is much weaker than that of A-exciton, its measured scattering spectrum is easily affected by the non-collimated components of the incident light and the irregular structure of the silicon nanoparticle. Thus, the experimental Rabi splitting energy differs from the simulated results.

Figure 12.Simulated angle resolved scattering spectra of the nanocavity with an embedded monolayer under the excitation of s-polarized plane wave.